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Interferenzoptische Versuche zur Plattenbiegung. 
Von R. Landwehr und G. Grabert +. 


1. Einleitung. Die Brauchbarkeit der Benutzung lichtoptischer Interferenzerscheinungen 
zur Untersuchung von geringen Deformationen soll hier an dem Beispiel der Durchbiegung 
-von Platten kurz gezeigt werden, um einmal wieder auf diese einfache Methode hinzuweisen und 
zu ahnlichen Untersuchungen anzuregen. Das Ziel war die Bestimmung der Form der Biegeflache 
und der GréBe des Biegepfeiles bei einigen 
einfachen Belastungsfallen und ein Ver- 
gleich mit der Theorie. Bei solehen bekann- 
ten Biegefallen laBt sich weiterhin mit 
Hilfe des interferenz-optisch bestimmten 
Biegepfeiles leicht der Elastizitatsmodul 
des Plattenmaterials mit relativ groBer 
Genauigkeit bestimmen. 


2. Versuchsanordnung. Benutzt wird 
das Verfahren der klassischen lichtopti- 
schen Interferenzen, das schon 6fter zu 
Biegeuntersuchungen, insbesondere auch 
zur Untersuchung der Poissonschen Kon- 
stante herangezogen worden ist!. Die Ver- 
suchsanordnung zeigt Abb. 1. Die Prif- 
lingsplatte P, die der Deformation unter- 
worfen werden soll, wird tiber einen halb- 
durchlassigen Spiegel mit monochroma- 
tischem Licht der Wellenlange A — z. B. 
einer Spektrallampe, eventuell mit vor- 
geschaltetem Filter — von oben beleuchtet. 
Die Beobachtung erfolgt im reflektierten 
Licht. Méglichst dicht tber der Platte P 
befindet sich die Vergleichsplatte V. Be- 
trachtet man die Oberflache des Priiflings 
oder bildet sie auf einer Auffangebene, etwa 
einer photographischen Platte, ab, so tritt 
bei geniigender Parallelitat und Planitat 
der Luftschicht zwischen P und V ein Sy- 
stem von dunklen Interferenzstreifen auf. 
Dies entsteht dadurch, daB sich jedesmal 
zwei benachbarte (koharente) Lichtstrahlen, von denen der eine an P, der andere an V reflektiert 
wird, je nach dem Gangunterschied, den sie haben, ausloschen oder verstarken. Ist das Bindel 
gentigend parallel und fallt es geniigend senkrecht auf die reflektierenden Flachen, so verbinden 
die Streifen die Stellen gleicher Dicke der Luftzwischenschicht; es sind ,,Interferenzen gleicher 
Dicke’‘. Da nun die Streifen immer bei Anderung des Gangunterschiedes der beiden Strahlen 
um / auftreten, entsprechend einer Anderung der Dicke der Zwischenluftschicht um A/2, kann 
man sie auffassen als Schnittlinien von Ebenen, die parallel zur Vergleichsfliche V sind und 
immer denselben Abstand (A/2 ~3-10-4mm) voneinander haben, mit der Priflingsoberflache. 
Bei guter Ebenheit der Vergleichsflache erhalt man somit ein Hohenschichtlinienbild der Ober- 
flache der Platte P. 

Im vorliegenden Fall wurden die Versuche mit einem Interferenzgerat® ausgefiihrt, das wegen 
des groBen Arbeitsraumes unterhalb der Vergleichsplatte und der bequemen und feinen Fin- 


Abb. 1. Versuchsanordnung (Strahlengang schematisiert). 


I Viel z. Be i. Kohlrausch, Praktische Physik, 5.79, 84, 19. Aufl. Leipzig und Berlin 1944, 
2 R. Landwehr, Z. Instrumentenk. 62 (1942), S. 73, 


2 Landwehr und Grabert }: Interferenzoptische Versuche zur Plattenbiegung. 
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stellméglichkeit besonders geeignet war. Als Lichtquelle wurde eine Na-Spektrallampe ver- 
wendet (A=589 my). Doch dirfte bei Beriicksichtigung der oben erwahnten Bedingungen auch 
schon ein einfacherer Versuchsaufbau brauchbar sein. Ein Interferometer, bei dem die Ver- 


Abb. 2. Interferenzstreifensystem (Héhenschichtlinienbild der 
Biegeflache) bei einer an drei Randpunkten in 120° Abstand 
frei aufliegenden kreisférmigen Platte (Belastung 1,83 kg). 


gleichsflache sich entsprechend der Michelson- 
schen Anordnung seitwarts befindet und in 
die Ebene der Priflingsoberflache gespiegelt 
wird, ist nattrlich auch hier vorteilhaft, da 
die Vergleichsflache nicht mehr stért und ge- 
schiitzt liegt und die Interferenzen dadurch, 
da®B in diesem Fall der Gangunterschied sehr 
klein gemacht werden kann, besonders scharf 
zu erhalten sind. Eine weitere Steigerung 
der Scharfe der Interferenzstreifen 1aBt sich 
auBerdem durch geeignete Verspiegelung der 
in Frage kommenden reflektierenden Flachen 
erreichen!, 

Als Priiflinge dienten im vorliegenden Fall 
drei polierte Glasplatten. Die Reflexion an 
der Unterseite war bei der Platte I durch 
Mattieren dieser Seite? unschadlich gemacht, 
bei der Platte JJI durch Benutzung eines 
Farbglases herabgesetzt. Natirlich kénnen 
derartige Untersuchungen auch an Platten 
aus anderem Material (z. B. Stahl usw.) 
durchgefihrt werden, wenn nur eine ge- 
nigend ebene und spiegelnde bzw. ver- 
spiegelte Oberflache vorhanden ist. Auch 


wenn die Ebenheit nicht so hochgradig und das Ausgangsinterferenzsystem daher unregelmabig 
ist, laBt sich die Gestalt der Biegeflache nach einem besonderen Verfahren ermitteln3. 
Zur Aufnahme des Priiflings diente eine besondere Belastungsvorrichtung B, die ein in sich 


geschlossenes System darstellte und 
moglichst stabil gebaut war. Die P 
obere Platte G aus Stahl konnte 


ausgewechselt werden. Sie war mit kg 


einer Offmung in der Mitte und 8 
einem glatten Rand bzw. den ent- 

sprechenden Auflagespitzen 
sehen, gegen die dann die jeweilige 


ver- 


Glasplatte durch den abgerundeten 
Stahlbolzen des Belastungshebels 
angedrtickt wurde. Die Aufbringung 


der Last erfolgte als Mittel-EKinzel- 4 
last von unten tiber ein eingebau- 
tes Hebelsystem H, wobei die Hohe 
der Belastung durch Eichung er- 


mittelt wurde. Die ganze Be- 
lastungsvorrichtung stand auf dem 
Tisch T des Interferenzgerats und 


konnte mit dessen Stellschrauben 


so eingestellt werden, dal} die ein- 
ander zugekehrten Flachen von 
P und V_ parallel zueinander 


‘Vel. z. B. S. Tolansky, Phys. Bl. 4 (1948), S. 472. 


5 70 15 20 ae 


10 ?mm 


Abb. 3. Abhangigkeit der gemessenen Biegepfeile f von der Belastung P. 


* Hier sei eine Erscheinung erwahnt, die auch bei sorgfaltigster Bearbeitung wahrend des Mattierens der 
einen Seite der urspriinglich doppelseitig polierten Glasplatte stets auftrat, und zwar bei allen Platten etwa 
im gleichen Ausmaf: Es zeigte sich, daB die nicht mattierte Flache konkavy wurde, die ganze Platte sich 
also — vermutlich als Folge verminderter Oberflachenspannung — durchgebogen haben muBte. 


* R. Landwehr und A. Dose, Naturwiss. 36 (1949), S. 342. 
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standen. — Abb. 2 zeigt das System der Interferenzstreifen (Hohenschichtlinienbild) fir den 
Fall der an drei Randpunkten frei aufliegenden kreisférmigen Platte. 


3. MeBergebnisse. Gemessen wurden nur die Biegepfeile durch Auszahlen der Interferenz- 
streifen. Die gemessenen Werte, die einzeln allerhéchstens bis zu 4%, streuten, wurden graphisch 
ausgeglichen. Abb. 3 zeigt die Proportionalitat der MeBpunkte mit der Belastung. An diesen 
verbesserten Werten wurde dann noch eine Korrektur angebracht, die die noch vorhandene 
Deformation der Belastungsvorrichtung beriicksichtigt. Dieser Anteil hegt in der GréBenordnung 
von einigen 1-10-4 mm. Diese korrigierten Werte sind als Messung in die folgende Tabelle | 
eingetragen. Die Tabelle enthilt die Zusammenstellung der vier untersuchten Falle mit den 
Ergebnissen. 


Tabelle 1. Gemessene und berechnete Biegepfeile bei den untersuchten Biegefiillen. 


—— xs FS 


Be- f-10~-*mm | f- 107? mm | per es f-10 2mm} ne 
Reihe Platte Auflage lastung | gemessen theo- a oe theo- EONS 
| ‘in kg | (korr.) | retisch eBeoret retisch theoret. 
= | Wertes Wertes 
eel eR mca a A a oi ee coe 
I (K) (H) 
R (ERs et er et 0183 3,62 3.68 i 3.55 + 2,0 
2 50mm ¥ 4,505 8,92 9,05 == Ihe! 8,75 +1,9 
5 aig 8.518 16.88 17,11 a 16.55 +20 
ms 2x =43,8 mm (Sch) 
B ee amoRand fret 1.83 9,45 10,59 ul ae es 
ym | aufliegend 4,505 23,37 26,07 10 at a 
2r=48,4 mm 
I : (K) (H) 
(BK 7) Shah 1,83 4,64 4,69 =i 4,57 ais 
C gm oh ee 4,505 11,42 11,55 Seine 11,25 4215 
; Abstand frei 
ees aufliegend 
y K (H) 
a 3) Pee et «0266 0,87 aes O184e lone og 
D 79,6 ana OS 7 Nor 9952 3.12 3,17 erie 3.03 Bes 6 
Tere Does el 1) 5965 TAI 7.54 7 7.22 +26 
h=4,9 mm eee 5,585 18,28 18,60 Pai 17,79 +28 


4. Vergleich mit der Theorie. Die Berechnung der theoretischen Werte erfolgte nach Féppl* 
bzw. Nddai?. Fir den Biegepfeil einer frei am Rande aufliegenden kreisformigen Platte mit dem 
Halbmesser r und der Dicke h gilt bei Einwirkung einer Mittel-Einzellast P 

| Pr? 3(1—p) (3+ yp) 
f= 4 
(E = Elastizitétsmodul, 4 = Poisonsche Zahl). tie . 

Fir den Fall CG der an drei Randpunkten aufliegenden kreisférmigen Platte wurde die Durch- 
biegung in der Mitte unter der Einwirkung einer Mitteleinzellast P in der von A. Nada angegebenen 
Weise (Auflésung der Plattengleichung in Polarkoordinaten mit Hilfe eines Reihenansatzes unter 
Berticksichtigung der bestehenden Randbedingungen) berechnet. Es ergibt sich ein Biegepfeil 

; Tigh Nek : l+p 07 
| aaaee ee (0,592 — 0.769) , | 
der noch zu dem obenstehenden Wert f, hinzuzufiigen ist. Der gesamte Biegepfeil ist also fe=fi tl ye 
Bei einer an vier Eckpunkten frei aufliegenden quadratischen Platte mit der Kantenlange 2a 
und der Dicke A erhalt man fiir die Pfeilhdhe unter dem Einfluf einer Einzellast P in der Mitte 


entsprechend der Entwicklung von Féppl* 


Pa 
i Ens 192 (1—?) 


wa —8 u 
78— 64 7? pu? 


1 4, und L. Féppl, Drang und Zwang, Bd. 1, 2. Aufl. Miinchen und Berlin 1924. 
2 4, Nddai, Die elastischen Platten, Berlin 1925. _ Seen tes ups ; 
o IDIG Formeln fiir die an drei Randpunkten aufliegende kreisformige und fiir die quadratische Platte 
wurden in dieser vollstandigen Form angegeben, um eine Berechnung fiir beliebige Poissonsche Zahlen, die 
manchmal, wie im vorliegenden Fall, erwiinscht ist, zu erméglichen. Aus den angefiihrten Literaturstellen 
> 5 > 


sind die Formeln nicht ohne weiteres zu entnehmen, 
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Der Berechnung der Biegepfeile wurden die von Kunert! und Hoesch? angegebenen Werte fir 
die elastischen Konstanten bzw. im Falle des Farbglases (B) der sich aus der Zusammen- 
setzung ergebende, nur sehr angenahert geltende Wert fiir den Elastizitatsmodul® zugrunde 
gelegt (K und H bzw. Sch in Tabelle 1). Leider konnten die Konstanten seinerzeit nicht an 
den Versuchsplatten selbst bestimmt werden. 

Tabelle 1 enthalt neben den berechneten Biegepfeilen die Abweichung der gemessenen Werte 
von den berechneten in Prozent der letzteren. In den Fallen A, C und D ist die Ubereinstimmung 
sehr gut, besonders, wenn man annimmt, daf die noch verbleibenden Abweichungen weitgehend 
reell sind. Diese Annahme liegt nahe, da hier die Abweichungen fir die beiden Wertepaare EF, 4 
verschiedene Vorzeichen haben und meist gréfser sind, als den méglichen Fehlern entspricht. 
Auch im vorliegenden Fall werden also die naherungsweise abgeleiteten Formeln, insbesondere 
auch fiir die an drei Randpunkten aufliegende kreisformige Platte und die quadratische Platte, 
erneut bestatigt. Bei der an drei Randpunkten frei aufliegenden kreisformigen Platte ist der Biege- 
pfeil um 28% gréBer als bei der am Rande frei aufliegenden Platte. Fiir das Verhaltnis der Grobe 
der Biegepfeile in diesen beiden Fallen zueinander fallen die durch die Deformation der Be- 
lastungsvorrichtung und durch die allenfalls noch vorhandene Unsicherheit des Elastizitatsmoduls 
bedingten Fehler fort, da es sich beide Male um dieselbe Platte und dieselben Belastungen handelt. 
Die mitgeteilten Versuche unterscheiden sich von den von Nddai durchgefihrten, abgesehen vor 
dem benutzten Verfahren, auch dadurch, daB® im vorliegenden Fall dickere Platten untersucht wur- 


den (das Verhaltnis h/2r ist im Mittel 1/20 gegentiber 1/120 bei Nddaz). 


5. Bestimmung des Elastizitaétsmoduls. Berechnet man umgekehrt den Elastizitatsmodul 
mit der Annahme eines mittleren w-Wertes (fir BK 7 w~=0,21, fir F3 und RG1 w=0,22), so 
erhalt man fir BK 7 und F'3 die in der Tabelle 2 angeftihrten Werte. Eine Fehlerdiskussion 
ergibt als mittleren Fehler des Mittelwertes den in der letzten Spalte angegebenen Betrag?. 


Tabelle 2. Gegeniiberstellung der gemessenen Elastizitatsmoduln (in kg/mm?) fiir zwei Glassorten. 


Beobachter | Verfahren 2s Cos Genauigkeit 
BK7 EF 3 
Keeani Ultraschall 719312 | 5579 41% 
Hoesch Ultraschall 8170 | 5893 EI 0:0% 
Landw.-Gr. Plattenbiegung 8015 5736 +1,1% 


Fir die Farbglasplatte wurde ein Elastizitatsmodul von 8101 kg/mm? gefunden (gegentiber 
einem zunachst eingesetzten Wert von 7378 kg/mm?). Allerdings ergibt sich fiir die Genauigkeit 
dieses Mittelwertes nur 2,5°% wegen der gréBeren Streuung der MeBwerte. Zum Vergleich ent- 
halt Tabelle 2 die von Kunert und Hoesch gemessenen Werte sowie die von ihnen erreichte 
MeBgenauigkeit. Die noch vorhandenen Schwankungen des Elastizitatsmoduls erklaren sich 
wohl aus dem Unterschied der verschiedenen Schmelzen. Jedenfalls ist eine Bestimmung des 
Elastizitatsmoduls unter Benutzung der Plattenbiegung mit etwa derselben Genauigkeit moglich 
wie bei Anwendung der Lichtbrechung an Ultraschallwellen bzw. bei dem Verfahren der se- 
kundaren Interferenzen. 

Auf der einen Seite laBt sich also mit Hilfe der Lichtinterferenzen der Elastizitatsmodul genau 
bestimmen. Auf der anderen Seite kann man aber auch, wie gezeigt wurde, bei beliebigen Fallen 
von Biegung oder Deformation eine Kontrolle theoretischer Ableitungen vornehmen oder, wo eine 
theoretische Berechnung nicht oder nur mit groBem Aufwand méglich ist, die Durchbiegung — auch 
tiber eine gréBere Flache hinweg — rein empirisch bestimmen bzw. auch unbekannte auf bekannte 
Belastungsfalle zuriickfihren. 


1 D. Kunert, Glastechn. Ber. 16 (1938), S. 383. 

2 K. Hoesch, Z. Phys. 109 (1938), S. 606. 

3 Nach einer freundlichen Mitteilung der Firma Schott u. Gen., Jena. 

* Bei sehr genauen Messungen ist unter Umstanden eine Korrektur wegen des iiberstehenden Randes er- 
forderlich und fiir den Fall der am Rande frei aufliegenden kreisférmigen Platte die strenge Formel (vgl. z. B 
A.E.H. Love, A Treatise on the Mathematical Theory of Elasticity. $.475, 4. Aufl. NewYork 1944) anzuwenden. 


(Eingegangen am 20. Dezember 1948.) 
Anschrift des Verfassers: Dr. Rudolf Landwehr, (21b) Menden (Kreis Iserlohn), Iserlohner StraBe 19. 
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Graphische Bestimmung der Spiegelbewegungen oe 
beim Differentialwasserschlo8& von Johnson. | 


Von H. Weirich. ity | : 


1. Einleitung. Die vorliegende Arbeit! befaBbt sich mit der Ermittlung der Spiegelbewegungen aay: 
im DifferentialwasserschloB von Johnson durch ein graphisches Integrationsverfahren, das durch 7 ae 
Erweiterung der von L. Miihlhofer? fir ungedampfte Wasserschlésser entwickelten zeichnerischen pes 
_Lésungsmethode gewonnen wurde. Dieses Miihlhofersche Verfahren, das bis zur zweiten Naherung 5 a 
genau ist, und von K. Karas® auch auf gedampfte Wasserschlésser ausgedehnt werden konnte, One: 
wurde deshalb als Ausgang fir die folgenden Untersuchungen gewahlt, weil es von den vorhan- 
denen graphischen Verfahren wohl das scharfste ist und fiir die Praxis die besten Ergebnisse liefert. 
Es gestattet die Bestimmung der Spiegelstande und der Wassergeschwindigkeit im: Stollen a 
fiir jeden Zeitpunkt des Schwingungsvorganges, und so war es naheliegend, diese Verfahren ANAT? i 
auch auf das DifferentialwasserschloB zu iibertragen. Durch wesentliche Erweiterungen ist es a 
gelungen, die Bewegungsvorgange auch in diesem Wasserschlo&typ hinreichend genau zu be- 
stimmen. Das beweisen die Ergebnisse der behandelten Beispiele, die durchweg mit einer vom 
Verfasser an anderer Stelle angegebenen analytischen Lésungsmethode oder mit Ergebnissen 
anderer Untersuchungen verglichen wurden. 


2. Allgemeines. Um beim gedampften Wasserschlof die im Durchtrittsquerschnitt im Becken- 
boden (Widerstand) auftretenden Dricke, die zu Beginn einer Belastung (Eréffnung) bzw. Ent- 
lastung (Drosselung) der Turbinen GréB8twerte annehmen und fiir Stollen und Druckrohrleitung 
plétzliche gefahrliche zusatzliche Dricke darstellen, zu verringern, ordnete Johnson im Becken 
des Wasserschlosses ein Zentralrohr an, das bei Belastungssteigerung im ersten Augenblick 
die fehlende Wassermenge liefert, bei Entlastung die iberschiissige Wassermenge aufnimmt. 


28 
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SOE, Abb. 1. DifferentialwasserschloB von Johnson. 
Langsschnitt, Bezeichnungen und Vorzeichenfestsetzungen. 
Dadurch werden die Driicke im Widerstand nicht plétzlich wirk- 
F sam, sondern erstrecken sich tber einen Zeitraum von einigen 
Sekunden und verlieren dadurch ihre gefahrliche Wirkung auf 
die Druckrohrleitung und den Stollen. 
3. Wirkungsweise. Aus Abb. 1 ist ersichtlich, daB die Be- 
Unterwasserspiege! harrungsspiegel im Becken und Zentralrohr am Anfang einer 


=== Eroffnung von Teillast um z, am Ende derselben um 2 unter dem 


1 Auszug aus der Dissertation. Deutsche Techn. Hochschule, Prag 1941. 

2 L. Miihlhofer, Zeichnerische Bestimmung der Spiegelbewegungen in Wasserschléssern von Wasser- 
kraftanlagen mit unter Druck durchflossenem Zulaufgerinne. Berlin 1924. 

3 K, Karas, Bauingenieur 22 (1941), S. 335. 
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Niveau des Stauweihers auf gleicher Héhe liegen. Hierbei sind Za und z die Verlustdruckhéhen, 
herrtihrend von der Reibung im Stollen. Werden die Turbinen belastet oder entlastet, so wird 
die fehlende bzw. iiberschiissige Wassermenge zundchst ausschlieBlich dem Zentralrohr entnom- 
men oder zugefiihrt, weil im Widerstand noch keine Druckhéhe p fiir den Wasserdurchtritt zur 


Verfiigung steht. Diese ist erst dann vorhanden, wenn der Zentralrohrspiegel etwas abgesunken ~ 


oder angestiegen ist, und ergibt sich, von Zusatzdricken abgesehen, immer aus der Differenz 
der Spiegelstande im Zentralrohr und Becken z.—%. . | 

Wegen der anfangs noch kleinen Druckhohe ist auch die durch den Widerstand gedrickte 
Wassermenge viel kleiner als jene, die von den Turbinen bei Belastung angefordert oder bei 


Entlastung abgegeben wird. Es wird deshalb dem Zentralrohr weiter Wasser entnommen bzw. 


zugefihrt und dadurch der Spiegel in ihm solange absinken oder ansteigen, bis schlieBlich jene 
Druckhéhe erreicht wird, die es erméglicht, daB die ganze Fehl- bzw. UberschuBwassermenge 


dem Becken entnommen oder zugefiihrt werden kann. Der diesem Zustand entsprechende — 


Spiegelstand im Zentralrohr wird schon nach wenigen Sekunden erreicht und beschleunigt oder 
verzogert den Stolleninhalt in kirzerer Zeit als beim gedampften WasserschloB. Dadurch wird 
aber der neue Beharrungsspiegel wesentlich schneller erreicht. 


4. Die Ausgangsgleichungen fiir das Differentialwasserschlo8. Ebenso wie J. Frank1 in seiner 
Arbeit vorgeht, werden auch hier die Gleichungen fir Be- und Entlastungsfalle getrennt ent- 
wickelt, weil bei Entlastungen verschiedene Zustande wegen der vollkommen und unvollkom- 
menen Uberfalle aber den Zentralrohrrand unterschieden werden miissen. 

a) Belastungssteigerung. Mit Vernachlassigung der Elastizitat des Wassers sowie der 
Stollen- und WasserschloBwandungen gelten bei konstanten Becken- und Zentralrohrquer- 
schnitten F, und F’, die Raumgleichungen 


dz, ibs Qr 

Wen ee @ 
dz, Qu 

dicrnn Byis (2) 


Hierin bedeuten dz,/dt und dz,/dt die Spiegelgeschwindigkeiten im Becken und Zentralrohr. 
Die durch den Widerstand strémende Wassermenge Qr in (1) bestimmt man nach J. Frank? 
mit Hilfe der Bedingung 


p=k=hy (54), (3) 


in der k bzw. ky die Druckverluste im Widerstand, die sich bei Annahme eines quadratischen 
Widerstandsgesetzes mit «, als Dampfungsbeiwert und den Stollengeschwindigkeiten v und vp 
beim DurchfluB einer Wassermenge Q und der Vollastwassermenge Q, zu av? und «,v? er- 
geben. Geht man in (3) mit P=%.—%, ky=a v5 und Qy=fvo, wobei f den Stollenquerschnitt 
bedeutet, ein, so erhalt man a 


Ore f(A, (4) 


Oy 


Die dem Zentralrohr entnommene Wassermenge Qr ergibt sich aus der Beziehung 


Qr=(q—Q) —Qr, (5) 


worin q die den Turbinen durch die Druckrohrleitung zuflieBende Wassermenge darstellt. Fahrt" 
man in (5) q=fc, wobei c jene gedachte Wassergeschwindigkeit ist, die vorhanden ware, wenn 


q allein den Stollen durchflieBen wirde, und Q=fv ein, so wird mit Beachtung von (4) 


Qiu efile v (2) : (6) 


ual 
In den Gleichungen (4) und (6) ist z,—z, immer mit dem Absolutwert einzusetzen. Wegen 
der fiir z gewahlten Zahlrichtung sind bei fallenden Spiegeln dz,/dt, dz,/dt und damit Qy und 
Qu — da (q—Q)—Qr > 0 — positiv, bei ansteigenden Spiegeln negativ. Dies erreicht man durch 


entsprechende Vorzeichenfestsetzung des Gliedes eSay wobei gilt: 
1 


positives | Aes cer ies ( Zo —2y i; fiir { fallende Spiegel 
negatives by | ansteigende Spiegel. 


' J. Frank, Dtsch. Wasserw. 27 (1932), S. 101. 
20}. Frank, -a;,a-O,, S2l02% 
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Die Stollenwassermenge wird durch die Zentralrohrspiegelhohe 2, beschleunigt. Wird kon- 
_stante Stauweiherspiegelhdhe angenommen und die Tragheit der Wassermassen im Weiher 
_ und WasserschloB vernachlassigt, so erhalt man mit Beniitzung der dynamischen Grundgleichung 
die Beschleunigungsgleichung Leds 

Face Tie ioe ames (7) 
In dieser Beziehung bedeuten L die Stollenlange, dv/dt die Beschleunigung der Stollenwasser- 
masse und h die Verlusthéhe im Stollen, die sich bei quadratischem Widerstandsgesetz mit dem 
Reibungsbeiwert « zu «v2 ergibt. 
a Strémt die Stollenwassermasse vom Stauweiher zum WasserschloB (v> 0), so wird die Druck- 
_ hohe z, um die Verlusthéhe h vermindert, bei umgekehrter Strémungsrichtung (v < 0) vergréBert. 
Es gilt somit 
<0% 
_b) Entlastung. Hier miissen verschiedene Zustande des Bewegungsvorganges im Becken 


_und Zentralrohr bei der Entwicklung der Grundgleichungen beriicksichtigt werden. J. Frank! 
-unterscheidet fiinf Zustande, die sich bei teilweisen Drosselungsfallen auf zwei reduzieren lassen. 


oberes | 


Vorzeichen bei wv? fiir v 
unteres | 


N am N N NI A 

N . |. vStauwether | N N N i N_vstouwerher 

i .: S \e V 

N N N N | A N 

\ eee eon 

N N N N | N N 

N N N N | A N 

S 
Wo : SARS ISSUES : 
RAMA AA NAAM An ) SEXY GK 
Abb. 2. Zustand I, Abb. 3. Zustand II. 


Zustand I. Dieser wahrt nach Abb. 2, solange die Spiegellage z, unter dem Zentralrohrrand 
liegt. Es gelten die Gleichungen (1) und (2). Da bei Beginn des Bewegungsvorganges die Spiegel 
ansteigen, sind dz,/dt, dz,/dt und damit Q7 und Qrr — da (q—Q) —Qi< 0 — negativ. 

Zustand IJ. Hat der Zentralrohrspiegel den Zentralrohrrand erreicht, so fallt die Wasser- 

-menge Qi7 in das Becken tber (Abb. 3). Wahrend dieses Zustandes liegt der Beckenspiegel unter 
dem Zentralrohrrand, so daf ein vollkommener Uberfall eintritt. Somit ergeben sich die Raum- 


gleichungen ee Ores ‘a 
dhe Fy, ; J “SEN 
diy a Qr1—Qi (9) 
dt LE eX, _ yStauweiher 


mit der Uberfallwassermenge Qi= > 4B | 2g hal : 
Hierin sind 1, der Uberfallbeiwert, B= D, die Uber- 
fallbreite mit dem Zentralrohrdurchmesser D, und 
ha die Uberfallhohe. Mit + 4, B|/2g—=C, folgt der fur 
die weiteren Betrachtungen einfachere Wert Qi=C, hii? 
Die GréBen Qr, Qi und Qi sind fiir diesen Zustand 
negativ. 

Zustand III. Sobald auch der Beckenspiegel den Zentralrohrrand erreicht «4, geht der 
vollkommene Uberfall in einen unvollkommenen tber (Abb. 4). Die Gleichungen fiir Becken 


Abb. 4. Zustand III. 


und Zentralrohr lauten jetzt da A. Ores (10) 
PES rae 
dz. ~ Gi Oa (11) 
de pe PUES ia es 


1\ J. Frank, a.a.O., S. 103. 
ON 
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me ‘ : , 2 ete 
wobei sich die Uberfallwassermenge mit dem weiteren Uberfallbeiwert ww. zu Qi = 3 4 B 2g hid 
+ Ul, Ba 2g hut ergibt. Mit ju, B /2s= C, wird Q,=C, hat + C,aha®. Die Gré®en Qr, Qrr und 


Q;, sind fiir diesen Zustand ebenfalls negativ. 

Zustand IV. Liegen die Spiegel im Becken und Zentralrohr schlieBlich auf gleicher Héhe, 
so wird im Falle vollstandiger Entlastung der Turbinen q=0 und in diesem Augenblick auch die 
Stollenwassermenge Q=0. Im nachsten Moment wird Q<0, d.h. es stroémt Wasser im Stollen 
vom Wasserschlo8 zum Stauweiher zuriick. Diese Wassermenge wird zunachst wieder, weil 
keine Druckhéhe p im Widerstand wegen der gleichen Spiegellagen im Becken und Zentralrohr 
zur Verfiigung steht, ganz dem Zentralrohr entnommen. Durch die entstehende Spiegelsenkung 
im Zentralrohr tritt aber ein unvollkommener Uberfall aus dem Becken in das Zentralrohr ein 
(Abb. 5). Auch fiir diesen Zustand gelten die Gleichungen (10) und (11), nur werden jetzt Qz, 


Qir und Qj, positiv. 
ti} 
| _ yStauweiher 4 
N 


CEs 


1 _yStauweiher 


hhh bh bb bh hdddeptbddddedd "\ 


CLISELELDBSEEEESESELDNESTEBDOEEDs 


KZ ZL La LLL LLL LLM ! 


NSS NSS 


Ort Gi - 
SSS 


waa 


N 
HOSS 


wreed 


Abb. 5. Zustand IV. — Abb. 6. Zustand V. 


¢ 


Zustand Vv. Ist der Zentralrohrspiegel bis zur Rohroberkante abgesunken, so geht der un- 
vollkommene Uberfall wieder in einen vollkommenen iiber (Abb. 6), und es haben erneut die 
Gleichungen (8) und (9) Giltigkeit. Die GréBen Q7, Qi und Q, sind nun jedoch positiv. 

Ist schlieBlich auch der Beckenspiegel bis zur Rohroberkante abgesunken, so herrscht wieder 
Zustand I mit den Gleichungen (1) und (2), allerdings mit positivem Q7 und Qu. 

Auch bei Entlastungsfallen gilt die Beschleunigungsgleichung (7) unter Beachtung des Vor- 
zeichens von «v* bei den einzelnen Zustanden nach den angegebenen Regeln. 


5. Darstellung der Ausgangsgleichungen fiir das zeichnerische Loésungsverfahren. In diesem 
Abschnitt werden die fiir das graphische Verfahren notwendigen Gleichungen fiir die Grundform 
des Differentialwasserschlosses entwickelt. Fir Varianten dieses Wasserschlosses, z. B. mit ein- 
gebautem Widerstand im Zentralrohr, mit Schlitzen in diesem, die Zusatzdriicke darstellen 
oder mit Uberlaufen ist das Verfahren ahnlich, nur werden die Gleichungen schwieriger. Wegen 
der Verschiedenheit des Verfahrens fiir die Schwingungsvorgange bei konstantem sekundlichen 
Wasserverbrauch und konstanter Leistung erfolgt im folgenden eine getrennte Behandlung der 
Vorgange. 


A. Schwingungsvorgange bei konstantem sekundlichen Wasserverbrauch. 
a) Belastungssteigerung. Aus den Gleichungen (1) und (2) entstehen mit Verwendung von 
(4) und (6) die Beziehungen 
Zy—2y \2 
dz, ( Oy ) 


PLR MGSCN haa hs (la) 
£ 
an hf Ba ONS 
dz, os way ( Oy ) 
dt Fy i (2a) 
f 


! —2,\8_. oa i : 
wobei (= 1) mit positivem Vorzeichen einzusetzen ist. Die Beschleunigungsgleichung (7) geht 
1 


‘ 
as WR Ar bE oP TO ye 
cel 


von (3) bei gleichzeitiger Beachtung von (1) in 
s . L dv | IPS OANA iG ea ae . - ea 
—_ ———- 2 fae wvitihy “is ad 
gow rea) (a) Ca) ae 
er. Damit gilt ebenfalls Reet 
eee BON dean tee 
4 — m= + (+) (3) = (3a) ee: 
= Cael 
_ Wegen (3) findet man in dem Glied «, (3) ey den Druckverlust k im Beckenwiderstand 1 a 
wieder, der sich beim Absinken des Beckenspiegels als Druckverlust im Stollen gegenttber dem 594 


Beckenboden darstellt und beschleunigend auf die Stollenwassermasse wirkt. Bei Spiegelanstieg 
gilt’ die umgekehrte Uberlegung. Fir die Vorzeichenfestsetzung folgt daher 


oberes : ; F,\2/d2,\2,... dz; (>0 
| | PRES Vorzeichen bei a, a (S) fiir Te Beis | 

In den Gleichungen (1a), (2a) und (7a) wurde die Zentralrohrdampfung nicht beriicksichtigt, 
_ weil sie in den meisten Fallen praktisch ohne EinfluB ist}. Sie wird auch in den weiteren Be- 
_ trachtungen nicht beachtet. 


= b) Entlastung. Zustand I. Es gelten die Gleichungen (la) und (2a), jedoch jetzt mit 


; - negativem (**)'.Gliea. 
a ; 1 


F.. Zustand II. Die Gleichungen (8) und (9) gehen wegen Os Ce uber in 
: , (Gre), oa 


4 | dz it Oy ne / eae 
‘4 1 : Ane F, i (8a) ' 
f 
as Presi ti ( Zy—2y Hy i Gq hed A 
4s dz Oy Shits oe 
2 yee a Seen Sey JP rea (Qa) — x 
a 


ey 


-wobei nach Abb. 3 die Beziehung ha=z,—h, mit der Zentralrohrhdhe h, gilt. Da Qu fiir diesen _ 


ee - : ’ Z_g— By, 5 . ; 
_ Zustand negativ ist, muB hia negativ eingesetzt werden. (fa) ist ebenfalls wieder negativ. 
a : nw? > al 


Zustand III. Dic Gleichungen (10) und (11) lauten hier Me , 
| a upiberhe ee Urea Me Ae 


dz Oy f ue (10a) 
dt F, j ? 
a C 
PENT as eas Sanat ea 
dzy c—o— ( Oy ) ek Bsn if an % (11a) 
Q dt F, 
q f 


EDie veranderliche GréGBe a ist immer absolut einzusetzen. Sie errechnet sich nach Abb. 4 zu 
a=2z,—h . Ebenfalls nach Abb. 4 findet man jetzt hi=z,—z,. Mit Hilfe der letzten Beziehung 


a 
rg 
mr. 
’ 

0 


«1 Handelt es sich um Konstruktionen, bei welchen auch im Zentralrohr ein Widerstand mit dem Dampfungs- 
od beiwert «, eingebaut ist, dann gehen die Gleichungen (1a), (2a) und (7a) durch Einfiihrung des Zentralrohr- 


x _dampfungsgliedes «, (3) (2) in die folgenden Gleichungen tiber: 


fe dt , 
F,\2 /dm\2 Fy)? ( dz)244 
: Zg—%y + X | (=) | ean eae ie ) (> ) i 
dz ae ; dz, oy 
iti Fy : dt RRs , 
Rafe f 

L dv F,\2 (da,\?__ (Sy (Sy 

rar eee faa NE) de)! 
1e\2_ cy: bere.) vc: eda sift 0 se 

Das %» (3)" (42) -Glied erhalt wieder das Sate \ Vorzeichen fiir re { a0" beter 
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gewinnt man aus (10a) und (lla) die fiir das zeichnerische Verfahren besser brauchbaren Glei- 
chungen 


Sad + (Bos : ) hat 


da _ f Da CRS (10) 
dt AGy 
ii 
3 C 1 
c—v Lee Git a hi? 
diy uf F Ya! : (11b) 


dt 


| 


3 ae . . . a . 5 ‘ . 
hi? und hy? sind negativ einzufthren, weil Qj, hier negativ ist. 


Zustand IV. Fir diesen Zustand gelten ebenfalls die Gleichungen (10b) und (11}). Da 


nun Q;, positiv ist, miissen hi? und hi? positiv genommen werden. Man kann auch hier eine Um- 
formung der Gleichungen vornehmen, wenn man nach Abb. 5 a=d—ha setzt, wobei 6 =2,—h, 
ist. Man erhalt dann 
(a rye a] hat OS CNet 
dz, if ay 
dt F, 3 


(10) 


1 a HG, 1 | 

Cpe ( COLT hae SHE hu? | 

dz f 2 1 & |) ) (11e). 
dt F, , 


f 


Der Wert 6 ist immer absolut einzusetzen. Die Uberfallhohe errechnet sich jetzt zu ha=2z,—2.. . 


Zustand V. Es kommen die Gleichungen (8a) und (9a) zur Anwendung. Zu beachten ist, 
dab (#=*\Fund hi hier positiv einzusetzen sind. hy=z,—h,. 


1 
Fur alle Zustande gilt die Beschleunigungsgleichung (7a) unter Beachtung der Vorzeichen 
Fy \2 /dz,\2 
i 2 Meal ey ees ks 
bei «& v? und GS (S$) : 

B. Schwingungsvorgange bei konstanter Leistung. Im Falle konstanter Leistungs- 
regelung muf sich die den Turbinen zuflieBende Wassermenge q mit der Nutzhéhe H, so andern, 
daB qH,=C=konst. ist. Die Konstante C errechnet sich aus der Leistungsgleichung 

1000 q Hn 
fee Raa lt 
worin 1000 das spezifische Gewicht des Wassers in kg/m*® und 7 der Turbinenwirkungsgrad 
sind, zu 15N 


(12) 


Damit wird die Geschwindigkeit 
: C 
bzw. mit H,=H,)—2z, nach Abb. 1, worin Hy das Rohgefalle ist, 
C 
SS . 14 
¢ Fie) oe 


Geht man mit (14a) in (2a), (9a), (11 b) und (11c) ein, so erhalt man sofort die entsprechen- 
den Gleichungen fiir dz,/dt beim Schwingungsvorgang konstanter Leistung. Die Gleichungen 
(1a), (8a), (10b) und (10c) gelten ebenso wie die Beziehungen fir hi, a und 6 bei Entlastungs- 
fallen auch hier ane pnlene Auch die Beschleunigungsgleichung (7a) bleibt erhalten. 

Besondere Beachtung ist wieder den Vorzeichen der einzelnen Glieder zu schenken. 


Die bisher gewonnenen Gleichungen zur Erfassung der Bewegungsvorgange dienen bereits 


fir das zeichnerische Verfahren. Nun sollen aber noch einige Gleichungen angegeben werden, | | 


die zu gewissen Zeiten des Schwingungsvorganges wesentliche Vereinfachungen des Verfahrens 
erméglichen. Diese Vereinfachungen bestehen darin, daB man wie J. Frank} in den gemeinsamen 


1 J. Frank, a.a.O.,'S. 105. 
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_ extremen Spiegellagen und deren Umgebung die sehr kleine Zentralrohrspeicherung vernach- 
i! lassigt und annimmt, daB die gesamte Fehl- bzw. UberschuBwassermenge allein dem Becken ent- 

nommen oder zugefiihrt wird. Diese Annahme ist zweckmaBig, weil in diesen Bereichen des 
Schwingungsverlaufes die kleine Druckhéhe p sich bei Beibehaltung des strengen Verfahrens 
nur ungenau bestimmen laft und dadurch auch z. B. die Zahler in (la) und (2a) ungenau werden. 
Frank fihrt dann die numerische Integration wie fir ein gediampftes WasserschloB ohne Zentral- 


_rohr durch, wobei er die Druckhéhe p bestimmt, die nun im Widerstand wirksam werden mu 
und damit die angenaherte Spiegellage z, erhalt. 


6. Gleichungen fiir das vereinfachte Verfahren. Unter Zugrundelegung der obigen Uber- 

_ legungen gewinnt man die folgenden Beziehungen, die auBer in den gemeinsamen Extremlagen 
auch im Bereich vom ersten Extrempunkt des Zentralrohrspiegels bis zur ersten gemeinsamen 
Extremlage angewendet werden kénnen, weil auch in diesem Bereich die Zentralrohrspeicherung 


klein ist. Die dabei gemachten Vernachlassigungen sind praktisch ohne Hinflu8 auf die gewon- 
nenen Ergebnisse. 


A. Schwingungsvorgange bei konstantem sekundlichen  Wasserverbrauch. 
a) Belastungssteigerung. Bei vernachlassigter Zentralrohrspeicherung werden dz,/dt und 
Qir gleich Null. Die Gleichung (5) geht ther in Qr=q—Q=f (e—v), und aus Gleichung (1) wird 
dz, c—v 
Tae F: (1b) 
J. 


Damit folgt die Druckhéhe p wegen (3) zu 
PSs 


E a (ev)? (3b) 
und der Spiegelstand im Zentralrohr . 


Zo = 24 + Hy (e—v)*. 


Das Vorzeichen von dz,/dt in (1b) ergibt sich aus der Differenz c—v mit Beachtung der Zahl- 
richtung fiir v. Fir das Vorzeichen von z, ist zu beachten: 

: positives | 

negatives | 


fallenden Beckenspiegel 
ansteigenden Beckenspiegel. 


Die Gleichungen (1b), (3b) und z.=z,-+«,(c—v)? finden auch fir den zweiten gemeinsamen 
Extrempunkt und seine Umgebung bei Entlastungen Anwendung. 


Vorzeichen von «,(c—v)? fir 


b) Entlastung. Fur Zustand I und V darf das vereinfachte Verfahren wegen der groBen 
Spiegelanderungen im Zentralrohr nicht verwendet werden. 


Zustand IJ. Die Gleichung (8a) bleibt bestehen. Der Wert (9a) gleich Null gesetzt, ergibt 
mit Beachtung von z.=h,+hi 


c—Uv 


1 i oa CG, S 

TE (hz t+hi—z2,)? 4 7 ha? . (15) 
Der Klammerausdruck ist absolut zu nehmen. Fiir die Be- 
stimmung von z, muB erst h; ermittelt werden. Dies geschieht 
mit dem bereits bekannten z, aus Gleichung (8a) nach 
Abb. 7 so, daB man in einem rechtwinkligen Koordinaten- 
‘system auf der Ordinate hy, auf der Abszisse nach links 


1 
Le und nach rechts c—v bzw. —— (hy + hi — 2)? auftragt. 
J ; ica 
Die a (h,-+-hu—2,)?- Linie kann durch Annahme von zwei 
Ot 
ha-Werten im Bereich des vermuteten hj leicht gezeichnet 


\ 
\ 
\ 
\ 
\ 
bed | 
\ 


yar (het ha-2y) ic 
aN Peas 
a 1 


werden. Die auf einem Transparentpapier aufgetragenen 


—1 ,3-Linie wird dann auf der Abszissenachse so lange ver- Abbi: eee events hi 
e ustan ° 


' schoben, bis c—v auf ihr abgeschnitten wird. Die Ordinate 
des Schnittpunktes ergibt die Uberfallhéhe hy. Mit diesem hy kann aus z,=h,-+ha die GroBe z, 
sofort bestimmt werden. 
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Zustand III und IV. Die Gleichungen (10b) und (10c) gelten unverandert. Die Werte 


(11b) und (11c), gleich Null gesetzt, ergeben 


C, (G I oe +] 
c—v = tha Be a+) hab, (0); 
cv = (Fo i —) hut , (C,—C,) hi? . | (17) 5) 


Das fir die Bestimmung von 2, notwendige hy erhalt man nach Abb. 8 und 9 ahnlich wie bei 
Zustand II. Die hierzu notwendigen Werte a=2,—h,=6 findet man mit dem aus (10b) und 
(10c) gewonnenen 2. Mit hy er- 
halt man z, fiir Zustand III aus haz 
Z=2,+hi bzw. fir Zustand IV 


aus. 2,=2,—hu 


B. Schwingungsvorgange 
bei konstanter Leistung. = 
a) Belastungssteigerung. 
Auch hier gelten die Gleichungen 
(1b), (3b) und z,=z,+-p mit Be- 
achtung von (14a). Das fir die 
Bestimmung von z, notwendige 
p kann aus Gleichung (3b) hier 
Abb. 8. Bestimmung der Uberfallhthe h,,  allerdings nicht gerechnet wer- Abb. 9. Bestimmung der Uberfall- 
bei Zustand III. den, weil ¢ nach Gleichung (14a) hohe h,, bei Zustand IV. 
eine Funktion des zunachst noch 
unbekannten z, ist. Man kann aber die Differenz z.—z, wie folgt-bestimmen: Die (2a) analoge 
Gleichung lautet fiir die Vorgange konstanter Leistung mit Beachtung von (3a) 


Cc i ( Bo 2 ip 
F,\2/dz\2 
dz» if, {Ho—[ato (=) (=) ]} “4 
asi oi 
ti 
Werden hierin die Gré8en von héherer als zweiter Ordnung der Spiegelgeschwindigkeit dz,/dt 
vernachlassigt, so erhalt man 
OA age $ fF (ree corny cara oy ea ‘ a_—2 \$ 
dz, f (Ho—41)? [Mo actos (FH) Gae) 4 ( Oy ) 2b 
dbcy Py : (2b) 


uf 
Dieser Wert, gleich Null gesetzt, ergibt 
cde AER eae + C (z.—2) Za—%, \4 
fo—4) Ff (Hy—21)? ( oy ). (18) 


v 


Aus der letzten Beziehung findet man 
nach Abb. 10 bei bekanntem z, aus (1b) den 
Wert z.—z,. Damit ist aber sofort z, und 
nach Gleichung (14a) auch c_ bestimmt. 

C (—21) 


D b Vv . h b ° BSENCS TGA Lies oe 
as obere Vorzeichen bei FH, ist fiir 


B44 )? 

val 
bereits angegebene Vorzeichenfestsetzung zu 
beachten ist. 


zu verwenden, wobei auch die fir ( 


C 
A ED) ; Am einfachsten und schnellsten bestimmt 
Abb. 10. Bestimmung der Spiegeldifferenz 2,—2;. man p aber aus der Beziehung (3a). Da 
auBer z, auch dz,/dt nach (1b) immer bekannt 
ist, kann p ermittelt werden. Damit sind aber nach z,=p+2, sofort z, und wegen (14a) c be- 
kannt. Die genaueren Ergebnisse erhalt man jedoch durch Bestimmung von p auf die erstere 
Art, besonders bei Verwendung groBerer Mafstabe. 


Spiegelsenkung, das untere fiir Spiegelanstieg - 


at 


aes Weirich: Graphische Bestimmung der Spiegelbewegungen. 


tl 


: rb, b) Entlastung. Auch hier darf fir Zustand I und V das vereinfachte Verfahren nicht ver- 
-wendet werden. 


_ Zustand ITI bis IV. Die Gleichungen (8a), (10b) und (10c) fiir die Bestimmung von z, 
haben unverandert Giltigkeit. Die Gleichungen (15), (16) und (17) gehen iiber in 


Cc 1 C ; 

ssiag Sf (Ho—hz—hi) Va (hz + ha — 24)? — Ff hi? , (15a) 
ea C C 1 Cau? 

Vv Fhe) ( fF a + Ve. ) hi? — 6 hi? 5 (16a) 

"= Faimaehy ~ (GO + pa) ht + Z(G C ya. (at) 


Die Bestimmung von hi geschieht wie bei konstantem Wasserverbrauch, nur ist z. B. bei Zu- 


stand IT in Abb. 7 die = (he--hu—z,)#-Linie durch die ky 


gu ersetzen. Mit der auf einem Transparentpapier aufgezeichneten — hy?-Linie wird jetzt v auf der 


hi ie (het hi —z,)*-Linie 
Oy 


Abszisse abgeschnitten und dann hy wieder als Ordinate des Schnittpunktes gefunden. Genau 
so verfahrt man bei Zustand ITI und IV nach den Abb. 8 und 9. Mit Ay erhalt man aber z, nach 
den frither angegebenen Beziehungen. 

Die Beschleunigungsgleichung (7a) gilt fiir das vereinfachte Verfahren unverandert. 

Mit den in Abschnitt 5 und 6 sowohl fir das strenge als auch fiir das vereinfachte Verfahren 
gefundenen Gleichungen fiir die einzelnen Bewegungsvorgange kann nun an die zeichnerische 
Lésungsmethode herangegangen werden. 


7. Vorbereitung des zeichnerischen Verfahrens. Entsprechend den Ausfihrungen L. Mihl- 
hofers1, die als bekannt vorausgesetzt werden, stellen auch hier die entwickelten Gleichungen fiir 
dz,/dt, dz./dt und (7) bzw. (7a) den Tangentenverlauf der Kurven (z,,t), (z),¢) und (v,t) dar. 
Es gilt somit, z. B. unter Verwendung der Gleichungen (la), (2a) und (7) 


Z—%, \$ 
: ee Meee aes ; 
eae Lie Doe 
of f 
By—2 \$ 
ay eee ae fads (19) 
Porieg at Bao ie. he! 
f f 
Sd ke OP ee 3 
UE deg eee i ‘ats 
g g 


Fir die Durchfihrung des Verfahrens ist die Bestimmung zusammengehorender Werte, 
hier also (qi, Yi, V1), (Pi+4> Yi+y Vitg)>--- der Winkel (p, y, ») fir die Punkte (Pi; Pi» Pi )s 
(Pi+3> Pisg> Pics)» --- erforderlich. Die Punkte (pj, pj’, pi’), (Pisa: Piva Pixg) --- sind Kurven- 
punkte mit den Koordinaten (z;, t;), (Zi) t;) umd (v;, t;) baw. (i+) tisg), (Z0i+4> t;,3) und 
(visa, tizg): Da die Punkte (pj, pj’, pi’), (Pix4s Pivk» Pisg)s ++» jeweils bekannt sind bzw. be- 
stimmt werden kénnen, so sind damit auch die Werte von z,, 2, und v, also (24, 2:5 Vi), (21443 
22445 Vi+4)> ... bekannt. Die Winkel q, y und y lassen sich leicht bestimmen, wenn die Groen 


2; ee F,/f, Zy, Zo, Fo/f, 2g «v? und L/g in den Gleichungen bekannt sind. Geht man zu endlichen 
Zeitdifferenzen At tiber, so kénnen die zugehérigen Anderungen Az,, 42, und Av in diesem 
Intervall nach L. Miihlhofer in zwei Schritten bestimmt werden. Als Ergebnis erhalt man die 
Z,-, Z- und v-Linien der Diagramme XIII und XIV in den spater folgenden Abb. 12 und 13. 


1 L. Mihlhofer, a. a.O., S. 28—30. 

2 Der Winkel p wurde hier abweichend von der Mihlhoferschen Darstellung gewablt, um Verwechslungen 
mit dem Beiwert « zu vermeiden. 

3 Zur Vereinfachung der Darstellung in der Zeichnung werden die Zahler der Gleichungen mit Z,, Z, 
bzw. Z bezeichnet. Bei Entlastungsfallen kommen auBerdem noch die Bezeichnungen Z, und Z, fiir die Zabler 
der Spiegelbewegungsgleichungen bei den Zustanden II bis V vor. 
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Da bei der allgemeinsten Form des Wasserschlosses auBer L/g, das immer konstant ist, alle 
GroBen der Bewegungsgleichungen veranderlich sein kénnen, so erfordert die Ermittlung der 
GroBen Z,, Z,, Fy|f, Zo, Zo, Fo[f und « v? acht weitere Teildiagramme I bis VIII, die anschlieBend 
behandelt werden. Zur Bestimmung der Druckhéhe p bei den vereinfachten Verfahren ist auBer-| 
dem das Teildiagramm IX notwendig. | 


Teildiagramm I; F,/f-Linie. Bei bekannter WasserschloBform kann die Funktion F,(z)) 
stets zeichnerisch dargestellt werden. Damit ist auch die F,/f-Linie bestimmt, aus der man bet 
zweckmaBiger Anordnung links neben dem (z, t)-Diagramm (Abb. 12 und 13) die zu den jer} 
weiligen Punkten der z,-Linie gehérenden Werte von F/f abgreifen kann. Hat das Becken eine; 
zylindrische Form, so ist die F’,/f-Linie eine zur z-Achse parallele Gerade im Abstande F/f., 


TeildiagrammII; F,/f-Linie. Es gilt das fiir die F',/f-Linie Gesagte, jedoch auf das Zentral-. 

rohr bezogen, auch hier. / 

Teildiagramm III; (44 \4 — Z,-Linie, Es ist zweckmaBig, diese Parabel wegen des; 

i 

doppelten Vorzeichens von (ae zu beiden Seiten der Achse zu zeichnen. Aus ihr kann bei) 

1 

bekannten z,- und z,-Ordinaten zusammengehorender Punkte der z,- und z,-Linien das zugehorige | 

Z, ermittelt werden. Der positive Parabelast dient zur Bestimmung der Z,-Werte fir Spiegel- | 
senkung, der negative fiir Spiegelanstieg. 


Teildiagramm IV; ? hi?- und GS ha2-Linie. Bei Darstellung beider Kurven in: 


einem Diagramm ist es aus Grinden der Deutlichkeit angebracht, jeweils nur einen Ast der-/ 
A 3 eae ‘ : F PEG, | 
selben zu zeichnen. Man findet aus diesen Linienziigen mit hy die Werte fir Fz hit bzw. 


C,—G. , F 
—2—_“1 hy? mit Beachtung der friiheren Vorzeichenfestsetzung. 


Teildiagramm V; hi?-Linie. Dieses Schanbild stellt hi? in Abhangigkeit von hi dar.) 
Weil hu? mit positivem und negativem Vorzeichen vorkommt, ist es zweckmabig, beide Aste der} 


Parabel zu verwenden. 


C, it G 1 
Teildi Vii (2 —— }- d {26 ——}-Linie. Diese Gerade liefert die Werte 
; el es oe! un & +7) inie lese Gerade liefer le Werte 
a+ — bzw. — 6 + —— fir die Gleichungen (10b) und (11b) bzw. (10c) und (Ile). 
Prean ie f a 


Teildiagramm VII; «v?-Linie. Die «v?-Parabel dient zur raschen Ermittlung der Verlust- 
héhe im Stollen in den Gleichungen (7) und (7a). Bei Belastungsfallen treten nur positive Stollen- 
geschwindigkeiten auf, so daf nur ein Parabelast fiir negative «v?-Werte gezeichnet werden 
muB. Bei Entlastungsfallen hingegen kénnen auch negative Stollengeschwindigkeiten vor-| 
kommen, also positive «v?-Werte nétig werden. Dann sind beide Seiten der Parabel zu zeichnen. | 
Es ist zweckmafig, dieses Schaubild links neben dem (v,t)-Diagramm anzuordnen, weil dann v_ 
gleich hertiberprojiziert werden kann. Beziglich der Bestimmung des Beiwertes wird auf die. 


Arbeit von L. Mihlhofer! verwiesen. 


Teildiagramm VIII; c= -Linie fir die Falle konstanter Leistung. Die Geschwin- 


pai Cea 
pe wo Fh) 
digkeit c ist im allgemeinsten Falle eine Funktion von t und z zugleich. Hier ist c nach Gleichung 
(14a) eine Funktion von zy allein. Als Schaubild erhalt man eine gleichseitige Hyperbel, die links | 


neben dem (z,t)-Diagramm anzuordnen ist, damit z, sofort heriberprojiziert werden kann. 


ae Z By Neate : 
Teildiagramm IX; «,(c—v)?- und «, Gat z}-Linie. Diese beiden Parabeln kommen nur 


beim vereinfachten Verfahren zur Anwendung und dienen zur Bestimmung der Druckhohe Pp 
im Widerstand. Bei den Fallen konstanten sekundlichen Wasserverbrauchs ist es zweckmabig, 
die «,(c—v)*-Parabel links neben dem (v,t)-Diagramm so anzuordnen, da® die Achse in die. 

; F,\2. | 
e-Gerade fallt. Die «, 2;-Parabel wird hingegen, wie Abb. 13 zeigt, gleich in das Rechen- | 


bild XI tbertragen, weil dort die Neigung der Polstrahlen auf einer Vertikalen in der Enid 


1 L. Mihlhofer, a. a.O., S. 5—7. 
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fernung eine Einheit vom Pol 0, sofort die z,-Werte liefert. Legt man die «, (=)" z?-Parabel so, 


da8B die Vertikale als Scheiteltangente und die Horizontale durch den Pol 0, als Achse dient, 
dann kann man fiir jedes 2, gleich das zugehérige «, (a3 =p ablesen. Wird, wie bei Bei- 


spiel 2 (Abb. 13), die Vertikale in einem gréBeren Abstand als eine Einheit vom Pol 0, gewahlt, 
so wirkt sich das natirlich auf den MaBstab von 2, entsprechend aus}. 

Wegen des doppelten Vorzeichens von «, sind diese Parabeln zu beiden Seiten ihrer Achsen 
zu zeichnen. Die Bestimmung des Beiwertes «, gestaltet sich fir Be- und Entlastungsfalle ver- 
schieden und erfolgt mit Hilfe der Formeln von F. Vogt?, J. Schiiller? und J. Frank’. Auf die 
richtige Ermittlung von «, ist groBer Wert zu legen, weil sonst der Verlauf des Zentralrohr- 
spiegels wesentlich von dem in den Abb. 12 und 13 abweicht. 


Teildiagramm X; p-Linie. Diese Linie gibt deutlich die verschiedenen Druckhéhen_ p 
im Widerstand wahrend des Schwingungsvorganges wieder. Zur Konstruktion der (z,t)- und 
(v,t)-Diagramme wird dieses Schaubild jedoch nicht bendtigt. . 


_ 8. Festsetzung der Vorzeichen fiir die Achsenrichtungen und der Mafstibe. Die Vorzeichen- 
festsetzung fiir die Achsenrichtungen in den Teildiagrammen geschieht nach den Angaben L. 
Mihlhofers®. Fir die hier hinzugekommenen Diagramme ist die Vorzeichenfestsetzung aus den 
spater folgenden Abb. 12 und 13 zu ersehen. 

Die Wahl der Mafistabe erfolgt ahnlich wie bei L. Miihlhofer®, jedoch ist wegen der raschen 
Spiegelbewegungen zu Beginn des Schwingungsvorganges besonders auf die Wahl des Zeit- 
maBstabes zu achten. Die Winkel p und y erhalt man nach (19), wenn die Zahler Z,, Z,, Za und UA 
auf der Ordinate des Teildiagrammes XI so aufgetragen werden, da8 1 mi/sek durch ED/T mm 
wiedergegeben wird. Der Winkel y folgt aus (19), wobei auf der Abszisse des Teildiagrammes XII 
1 sek? der GroBe L/g durch TD/G mm dargestellt werden muB. AuBer den gewahlten Haupt- 
mafstaben fiir den gesamten Schwingungsverlauf ist es angebracht, fiir gewisse Zeitintervalle 
und Teildiagramme abweichende Mafstabe zur Erhéhung der Ablesegenauigkeit zu verwenden 


(Abb. 12 und 13). 


9. Durchfiihrung des zeichnerischen Verfahrens. Bei der Bestimmung der Spiegelbewegungen 
und des Stollengeschwindigkeitsverlaufes geht man so vor, daB zunachst mit Hilfe der Ausfiih- 
rungs- und BetriebsgréBen der Widerstandsbeiwert a, ermittelt wird. Dann werden die An- 
fangslagen p’, p” und p’” der Spiegel und der Stollengeschwindigkeit festgelegt. Nach der Auf- 
zeichnung der jeweils fiir den Be- bzw. Entlastungsfall notwendigen Teildiagramme und der 
Festlegung der Mafistabe kann an die Konstruktion der 2,-, 2,- und v-Linien in Verbindung mit 
der Niederschrift herangegangen werden. Die Durchfiihrung des ersten sowie der weiteren Schritte 
ist nach den Ausfiithrungen von L. Miihlhofer’ leicht verstandiich. 

Abweichend davon ist hier zu beachten, daB wegen der raschen Spiegelbewegungen, be- 
sonders zu Beginn des Schwingungsvorganges und da wieder hauptsachlich im Zentralrohr, 
zunachst mit kleineren Zeitintervallen gerechnet werden mufi. Es ist zweckmaBig, mit 14/; bis 
1/,, des normalen, im weiteren Verlauf des Bewegungsvorganges beniitzten Zeitintervalles At 
zu beginnen. Auch in den gemeinsamen Extrempunkten soll wegen der kleinen Druckhéhe p 
die Schrittlange verkleinert werden, weil sonst 4z,/At und Az,/At ungenau werden. Durch die 
Wahl gréf8erer Mafstabe in diesen Bereichen kann wohl diesen Ungenauigkeiten noch eine 
Zeitlang begegnet werden, dann muf man aber zum vereinfachten Verfahren tbergehen. Nach 
dem gemeinsamen Extrempunkt kann bei geniigend grofem p das strenge Verfahren wieder 


1 Wird au®er der Beckendampfung auch die Zentralrohrdampfung beriicksichtigt, dann kommt als wei- 
teres Teildiagramm die «, (42)’2}-Parabel hinzu. Beziiglich Verwendung und Anordnung derselben gilt das 


af 

iiber die (4) #-Linie Gesagte. 

2 F. Vogt, Berechnung und Konstruktion des Wasserschlosses. Stuttgart 1923. 

8 J. Frank und J. Schiiller, Schwingungen in den Zuleitungs- und Ableitungskanalen yon Wasserkraft- 
anlagen. Berlin 1938. 
4 Siehe FuBnote 1 auf Seite 6. 

.5 L, Miihlhofer, a. a. O., S. 39. 
6 L. Miihlhofer, a.a.O., S. 36. 


7 L. Mihlhofer, a.a.O., S. 39—43. 


Zusammen- t c—v—( iets. Ny 
; ite spans . (2a) 
gehorende tgy Falf Falf 
P kt / Mt “st 
unkte p pp ; Rea | Zz. Zs F,/f F,/f 
sek m/sek Einh. 
1 2 =) 4 5) 6 7 8 9 | 
Sea AGE AERO ORT ONY Sateen) Meet Date eo ere Pe Dray Tio eee Menon Ee Eee TNS e ye Sea re 
0 0 +4,240 | —4,240 0 — 4,240 
0+4 4 +4,240 | —4,240 | —1,541 —2,699 11,76 0,49 
1 4 e=0 | 44930 | —4,230 | —1,728 | —3,502 | Einh, | Einh 
1+4 4 +4,222 | —4,222 | —2,090 | —2,132 
Zahlentafel 2. Niederschrift } 
GaGa ye | 
rhea ad OBEN CUM a .. . (10b) 
Ff Ff | 
heh Cip 8 (Fa —) hit ibe . 
ii if Vey Z, 
tg y= Ff = Fp 7 tb) | 
Zusammen- ae Co 3 1 4 f / 
gehérende | * Fe oa cenz (10) 
t = = o nellls 
Liss ex F,/f F,/f 
Pp’ p yee eh Fetietar eee icy 2 16 
c—v = fF hi ae G a+ =) a ( ) 
ai 5 ms 
c » c—v hud C a+r) Bi Z Ff | F./f 
hud 
sek m/sek Einh. 
1 2 3 4 5 ned Vpn 8 9 10's od 
27+4 25,5 +1,292 | -1,292 | -0,200 -1,079 -1,279 | —0,013 
28 26 0 +1,244 | -1,244 | -0,167 -1,064. | —1,231 | -0,013 | 11,76 0,49 
28+4 27 a +1,135 | -1,135 | -0,113 —1,022 -1,135 0 Einh, | Einh. | | 
29 28 +1,033 | -1,033 | -0,075 —0,958 —1,033 0 


angewendet werden. Wahrend der Ubergang vom strengen zum vereinfachten Verfahren immer __ 
leicht méglich ist, gestaltet sich der umgekehrte Vorgang nicht so einfach. Dieser soll deshalb | 
fir einen Belastungsfall nach dem ersten gemeinsamen Extrempunkt an Hand der Abb. 11 


gezeigt werden. 


Bis zu den Punkten p;,.,, und p;,,3 wird das vereinfachte Verfahren beniitzt, fir das die 
Beziehungen (1b) und z,=z,+p gelten. Von diesen Punkten ab soll das strenge Verfahren | 
wieder angewendet werden. Da bei diesem die Spiegellage im Zentralrohr z. nach (2a) bestimmt _ 
wird, so muB zunachst der Zahler Z, des Bruches gefunden werden. Dies geschieht nach Abb. 11 | 


so, daB ae im Teildiagramm XI die parallelen Strahlen 0,—i+2 zu p;’,,—p;,, und 0,—i+2+5 ' 
ZU Piio—Piio,3 Zieht. Halbiert man den Abstand i+2-i4242, so erhalt man den Punkt | 
i424 und damit einen angenaherten Wert von Z, fir den Zeitpunkt bei p;.,. Damit ist aber, 
dacund v bekannt sind, auch der Zahler Z, in (la) bestimmt. Daraus laBt sich eine neue Spiegel- — 


lage z, fiir das strenge Verfahren, dargestellt durch den Punkt p;,, rechnen!. Mit den Zahlern 


* Die Bestimmung der neuen Spiegellage z; beim Ubergang vom vereinfachten zum strengen Verfahren _ 
beim Entlastungsfall nach dem ersten Extrempunkt geschieht wegen des Uberfalles aus Zy. 
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zu Beispiel 1, Zustand I. 


tgy =—2 = = <S (7) 
Anmerkungen 
Zy | ao p=25;—2, Fav? Z L/g | 
m sek? 
toe el Sh 1) 12 13 14 15 16 


41,640 | +1,640 —1,640 0 
+1,640 | 40,218 | —1,422 | —1,640 | —1,422 42,84 aon 

p+-1,597 | —0,190 | —1,787 | —1,636 | —1,826 sek? Begins dee. Bevenahesy pe: 
erste) -=1,044 |. —9,616 | 1,630. | 22,674 Banges (Cin (La), (aa) .(1): 


zu Beispiel 1, Zustand IIT. 
ee a Ne LN 


Anmerkungen 


—4,003|}—4,327| -0,324 | 0,324 |-0,569| 0,333 1,895 |—0,153)—4,480 Duster Tg \ 


Ubergang vom strengen zum 


—4,053|—4,340| -0,287 | 0,287 |-0,536| 0,383 1,986 |—0,141|—0,481| 42,84 | vereinfachten Verfahren 


mm = = ‘ 2: 2 2 | Gin. (10b). (11b) u. (7) str. Verf. 
4,157|—4,378 0,221 | 0,221 |—0,470) 0,487 2,174 0,118|—4,496) sek Tae {rob}, Oye, Ly Ome Vert 


—4,243)—-4,412| —0,169 | 0,169 |—-0,411| 0,573 2,330 |—0,098)—4,510 Bi weh eat avait aides 


Z, und Z, kann nun von den Punk- 
ten p;,, und p;,. ab das strenge 
Verfahren wieder angewendet wer- 
den. Der Ubergang vom verein- 
fachten zum strengen Verfahren 
mu} immer dann erfolgen, so lange 
der von p;,. auf p;,, auftretende 
Sprung praktisch noch bedeutungs- 
los ist. 

Die genaue Festlegung der Ex- 
tremwerte der z- und _ v-Linien 


erfolgt nach den Angaben L. Miihl- Abb. 11. Ubergang vom angenaherten zum genauen Verfahren nach dem ersten 
h fer og HLS Reena gemeinsamen Extrempunkt bei einem Belastungsfall. 
0 i} . 


10. Beispiele fiir die Spiegelbewegungen im Differentialwasserschlo8. Aus einer groBen Anzahl 
Jurchgezeichneter Beispiele wurden zwei ausgewahlt, welche die Anwendung des Verfahrens 


1 L, Mihlhofer, a. a.O., S. 50, 


2 Fe basi | TaRectoreroe Joe yse]JOA OA Zunysepu e707 SYONVIqIaAIOsse\ UOYOI[puNyes urequVysuoy rq SuvS10assunsurmyog SOT : se Pesut 
: Gt 
at 
: W 
; ne 
A if IX 
ct 


CHE 7-2 


Sb Ale + 


Se en ee Blu] 7-©2 er 
Ste 5 bE A | = Ss 
eee $0 0 50 pS iN 49S $2/ Pop 
on pours ae nase an 0 
Sea faa $4D-% : 
WUT-*Z 
L- 
B= 
AT pee 


OS UW 2 67L = 
é [Pee] 
0 


“BunieS1aissFunyse[agq estaM[193 
{Sunqstey Jqueisu0y 13q 
SuesioAssunsuimySg "ET “qqVy 
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Zahlentafel 3. Niederschrift 


49-21 i: 
ZL, 
RN NA aI eee aval 
0 ap RF (t2) 
C za—% \§ 
v 
FS (Ho—22) ( ay ) eee | Z, 
t — 9 seus t3 2a 
by F,/f F,/f | we 
Zusammen- C 

gehérende t __ f(Ho—#2) es on A lb 
Punkte oe. F,/f F,/f ee 

TT LT. 25 
P’ pp pax (2)*4" eee) 


Cc G 
AS RE al v— Z. Zz z F. F. 
°F (Ho-#2) . Rar as (Ho-#1) f (Ho-1) a ‘ : lf olf 
m/sek 
0 0 | +2570 | +0,790 | +1,780| — us Ot 1780) 
0+4 4 +2,612 | +0,790 | +1,822 — — +0,739 | +1,083 — 
i 4 +2,623 | +0,796 | +1,827 — — +0,774 | +1,053 — 
1+4 4 +2,647 | +0,800 | +1,847 == — +0,943 | +0,904) — 2 g 
o a 
36+4 57 +2,435 | +2,512 | —0,077 | +2,437 +0,075 —0,077 0 —0,0118 3 F 
37 58 +2,435 | +2,438 | —0,003 | +2,435 +0,003 —0,003 0 0 a= 5 
37+4 59 +2,435 | +2,364 | +0,071 | +2,435 —0,071 | +0,071 0 +0,011 Z en 
vs a 
38-44 61 +2,447 | +2,228 | +0,219 | +2,443 | -0,215 | +0,219 0 +0,034 S S 
39 62 +2,456 | +2,164 | +0,292 | +2,447 —0,283 +0,292 0 +0,045 
39 62 +2,456 | +2,164 | +0,292 — — +0,233 | +0,059 — 
3944 63 | +2,462 | +2,104] +0,358| — SO) 0,905.1! 0:63) er a 
Zahlen- 
a nc a a a a 
AusfiihrungsgréBen BetriebsgréBen 
msek—} 
Belastungsfall 48 Fm Fea Maes cet sek-1 2m os sck1 |» = Qa i 
Cm4sek-1 


2 2 sek2 m1 
ee F, se h, ie ual cml sek-1 z 
2 


e 


ne 
m oz’ sek-1 msek— Q 
i 


2 : 
Totale Entlastung | 499 | 332,55 2 1 9.0914 |, 30.62 1,64 120 4.97 | 
b. konst. sekund]. 2 pale 
Wasserverbrauch 28,27 | 13,85 —3,67 0,599 ol,2 0 0 0 } 
: 3 
efiee eis 200 24,83 20 0,0752 Vis 0,047 3 0,79 Pre 
b. konst. Leistung 3,8 3,47 _ 0,599 — 0,52 10 2,63 : 


j—S— Nd 
1 In der ersten Zeile dieser Spalten sind die zeichnerischen Ergebnisse, in der zweiten die Vergleichswerte | 


durch eigene Rechnung und in der dritten die Vergleichswerte anderer Autoren wiedergegeben. 


® Die Angaben fiir diese Anlage wurden dem angefiihrten Buche von J. Frank und J. Schiiller, S. 189 ent- 
nommen. Dort findet man auch den in der Zahlentafel angegebenen Vergleichswert. Denselben Wert erhalt 
man aus der genannten Arbeit von F. Vogt, S. 93, Gleichung (19). 


8 Fiir dieses Beispiel wurden die Angaben von J. Frank, Dtsch. Wasserw., S. 105, beniitzt, der die an- 
gegebenen Vergleichswerte durch numerische Integration gefunden hat. 


ay 


zu Beispiel 2. 


e 
— 
: 
i 
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Anmerkungen 
Dl maa ‘ 
| 4 | tetya| Fee | z | ue 
m sek? 
aes ieee 15 16 17 18 19 
+0,047 +0,047 0 —0,047 0 : 
+0,047 +0,374 +0,327 —0,047 +0,327 Beginn des Bewegungsvorganges 
+0,086 +0,445 +0,359 —0,048 +0,397 Gl. (1a), (2a) und (7). 
+0,106 +0,638 +0,532 -0,048 +0,590 
~1,024 —1,028 —0,004 —0,475 —1,503 2. Ext kt 
~1,030 ~1,030 0 ~0,448 ~1,478 oy Gi aoe fe 
~1,030 —1,027 +0,003 0,421 1,448 C LE) ane): 
—0,984 —0,955 +0,029 —0,374 —1,329 ie Ubergang vom vereinfachten 
—0,939 —0,888 +0,051 —0,352 —1,240 zum strengen Verfahren 
nach dem 2, Extrempunkt 
—0,9206 —0,888 +0,0326 —0,352 1,240 * Gln. (1b), (3a) u.(7) vereinf. Verf. 
—0,878 0,826 +0,052 — 0,334 -1,160 Gln. (1a), (2a) u.(7) str. Verf. 
tafel 4. 
_Angenommene | Abgeleitete Ergebnisse fiir die extremen Ergebnisse fiir die 
MaBstabe Malistabe Spiegellagen! Stollengeschwindigkeit 
ae Imsek1 
| | Einheit ] sek ops Gs ee ae pmsek- ae aa oe 
erage sire Te 1 max= sek 7] max= psek | ,'m {sek 1 aus pe aus a i 
2 2max,| 72 . ige ere eigen 
ee pana pha Doak : Pocaas - tose eit Rechng. ane Rechng. 
20 6 100 fee a ar ae = 4086) 3 11 ,76:,| fin £20 oeks | Fite 175 sek 
30 2 a eee ance Cale area Wel tae, ae othe OBIS | Sk BWIeeT Bi) also aro eats 
; —4,62 — — == = a 
] —1,03 58 3,01) 10,5 
12 6 100 oN ca ne +3,92 | fiir t=10 sek | fiir t=38 sek 
Be Metal Se hel Ora peeia Aaa. area Taare sai aces 
50 50 6 | +303 | 19 | —1,09 | 57,5 |+3,03| 10 , 


»zeigen und seine Brauchbarkeit nachweisen sollen. Im ersten Beispiel, das Abb. 12 wiedergibt, 
wird ein Fall totaler Entlastung von Vollast (Qa=Q,) auf Betriebsstillstand bei konstantem 
'sekundlichen Wasserverbrauch behandelt. Das zweite Beispiel in Abb. 13 zeigt einen Schwin- 
-gungsvorgang konstanter Leistung hei einer Belastungssteigerung von einer kleineren Anfangs- 
ileistung N.(Q.) auf eine gréBere Endleistung IN (02). 

In den Abb. 12 und 13 wurden die Konstruktionslinien fiir den Beginn. zweier Intervalle in 
Jallen Teildiagrammen eingetragen. Die gleichzeitig mit den Zeichnungen entwickelten Nieder- 
ischriften sind auszugsweise in den Zahlentafeln 1 bis 3 wiedergegeben. Dabei zeigen die Zahlen- 


| 3 
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tafeln 1 und 2 einige Intervalle wahrend der Zustande I und III fiir das erste Beispiel, die Zahlen- 
tafel 3 den Beste und die Umgebung der zweiten gemeinsamen Spiegellagen ey zweiten Bei- 
spiel. Sie dirften ohne weitere Erklarung verstandlich sein. 

Die Ausfiithrungs- und Betriebsgréen beider Anlagen kénnen ebenso wie die angenommenen 
und abgeleiteten Mafstabe der Zahlentafel 4 entnommen werden. Dort sind auch die Extrem- 
werte der ersten gemeinsamen Spiegellagen max = omax ZUM Leit. tj, der zweiten gemeinsamen 
Spiegellagen 2 Z1max = Zomax Zur Zeit t., sowie die erste extreme Spiegellage im Zentralrolir 2 max 
zur Zeit t, zusammengestellt. Desgleichen findet man dort den Extremwert der Stollengeschwin- 
digkeit v, zur Zeit t;, und zwei Werte fir die Stollengeschwindigkeit zu beliebigen Zeitpunkten, 
die, ebenso wie alle anderen Ergebnisse mit Vergleichswerten belegt sind. 


11. Zusammenfassung. Die Ergebnisse der Zahlentafel 4 zeigen fiir zwei in der Praxis vor- 
kommende Falle, daB nach dem behandelten zeichnerischen Verfahren die Spiegellagen und die 
Stollengeschwindigkeit fiir jeden Zeitpunkt des Schwingungsvorganges hinreichend genau be-— 
stimmt werden kénnen. Dabei wurde immer eine plétzliche Belastungsanderung der Turbinen 
angenommen. Fiir den praktischen Betrieb ist jedoch auch eine Erfassung der Falle allmablicher | 
Belastungssteigerung bzw. Entlastung erforderlich. Diese Falle konnen mit diesem Verfahren | 
nach einigen Uberlegungen ebenfalls leicht bewaltigt werden. 

Auch bei Annahme eines anderen Widerstandsgesetzes als des quadratischen fur die Stollen- 
reibung und die Dampfung im Durchtrittsquerschnitt behalt das Verfahren seine Giltigkeit, | 
was als weiterer Vorteil zu werten ist. | 

Die Frage der Stabilitat wird durch geniigend weites zeichnerisches Verfolgen des Schwin- | 
gungsvorganges beantwortet. | 


(Eingegangen am 21. Dezember 1948.) 
Anschrift des Verfassers: Dr.-Ing. Hugo Weirich, (14a) Stuttgart W, SilcherstraBe 3/II. 
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Wéolbkrafte in diinnwandigen Profilstiben. 
Von W. Fliigge und K. Marguerre. 


, 1. Problemstellung. Wenn sich in einem zylindrischen Stab die zu den St-Venantschen 
_ Torsionsspannungen gehérige Verwélbung der Querschnitte infolge einer Endeinspannung oder 
einer Veranderlichkeit des Torsionsmomentes Mr langs der Stabachse x nicht frei ausbilden 
kann, so entsteht eine Stérung des Torsionsspannungszustandes. Sie ist gekennzeichnet durch 
das Auftreten einer ,,Woélbkraftgruppe, d.h. eines Systems von Langsspannungen ohne Re- 
_ sultierende und ohne resultierendes Moment, die, im Stabquerschnitt angreifend, Verwélbungen 
hervorrufen oder schon vorhandene Verwélbungen verandern (z. B. an einer Einspannstelle 
zu Null machen). Ihre Untersuchung ist die Aufgabe der Wolbkraft-Torsionstheorie, mit der 
- wir uns hier befassen wollen. 


In einem zylindrischen Stab, dessen Querschnitt in seiner Ebene undeformierbar ist, gibt es, 
wie unsere Rechnung zeigen wird, bestimmte Wélbkraftgruppen, die sich beim Hineinwandern 
in den Stab nur der Amplitude, nicht aber der Form nach verandern. Da sie sich als die Eigen- 
funktionen einer Differentialgleichung erweisen werden, nennen wir sie Eigenkraftgruppen. 
Jede solche Eigenkraftgruppe hat eine ganz bestimmte Abkling-,,Geschwindigkeit‘‘, und man 
erhalt den zu einer beliebigen Storung gehorigen Spannungszustand im ganzen Stab, wenn 
man die Stérung nach den Eigenkraftgruppen entwickelt, d.h. die Eigenkraftgruppen an der 
Stérstelle mit solchen Amplituden zusammenfiigt, daB die vorgegebene Stérung entsteht. (Grund- 
problem: an den Endflachen x=0, x=I eines Stabes greifen beliebig vorgeschriebene Normal- 
spannungen an.) Durch Uberlagerung eines solchen Stérzustandes tiber die aus der Theorie der 
ungestérten Torsion (oder auch Biegung) bestimmten Spannungen und Formanderungen erhalt 
man dann im Einzelfall die Lésung eines unmittelbar gegebenen Problems. 


Wir beschaftigen uns hauptsachlich in dieser Arbeit mit dem so definierten Eigenwert- 
problem. Wieweit das lésbar ist, hangt ab von der Gestalt des Querschnitts: Fir Voll- 
querschnitte (wo das Stérungsproblem aber gliicklicherweise so gut wie bedeutungslos ist) 
ist eine Lésung nicht bekannt. Der Rechnung zuganglich sind bisher nur die dinnwandigen 
Profile (wo die Querschnittspunkte nicht durch zwei gleichberechtigte Koordinaten festgelegt 
werden, sondern im wesentlichen nur durch eine), und bei diesen besteht wieder ein Unterschied 
zwischen den geschlossenen und den offenen. Bei den geschlossenen ist, wenn sie nicht sehr 
kraftige Gurte haben, die numerische Auswertung der allgemeinen Theorie recht mihsam, bei 
den offenen so einfach, daB man nach Einfiithrung gewisser, die Torsionsverwélbung enthaltender 
Querschnittsfunktionen zu ubersichtlichen geschlossenen Formeln gelangt. Diese Formeln sind 
seit langem bekannt!. Wir begniigen uns hier damit, zwischen diesen Formeln und unseren 
allgemeinen Ansitzen den Zusammenhang herzustellen; denn das Ziel unserer Arbeit ist es, 
die Frage der Wélbkrafttorsion der Rohbre zu klaren und die Wélbkrafttorsion der offenen Profile 
nach derselben Methode so darzustellen, daB erstens der Zusammenhang zwischen beiden Pro- 
blemen deutlich wird, und zweitens die erwahnten Formeln fiir die offenen Profile eine neue Be- 
griindung erfahren?. 


2. Die Grundgleichungen. Wir untersuchen das Endstérungsproblem unseres prismatischen 
diinnwandigen Stabes als ein Stabproblem, d.h. unter der Voraussetzung, daB ungeachtet der 
sonstigen Verformungen die Gestalt des Querschnitts erhalten bleibe. Wie weit diese Voraus- 
setzung zutrifft, untersuchen wir hier nicht. Bei Rohren insbesondere kann, wenn die Wandung 
sehr din ist, eine groBere Anzahl starrer Spante no6tig sein. 


1 Vel. z. B. R. Kappus, Luftfahrtforschung 1937, S. 444. 

2 Die vorliegende Untersuchung ist eine in gemeinsamer Arbeit entstandene Erweiterung eines friiheren 
Aufsatzes von W. Fliigge, Torsion mit behinderter Verwélbung, der im wesentlichen die Lésung des Wélb- 
kraftproblems fiir das Rohr enthielt. Es war zur Veréffentlichung im Prandtl-Festheft der Deutschen Akademie 
- der Luftfahrtforschung (1945) vorgesehen und ist, soweit uns bekannt, nicht mehr in die Offentlichkeit ge- 

kommen. — Fiir gelegentliche Diskussionen und kritische Bemerkungen méchten die Verfasser Herrn R. Kappus 
ihren Dank aussprechen. 
3* 
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Die Linie, die die Dicke t der Profilwandung tberall halbiert, nennen wir die Skelettlinie des 
Profils. Wenn wir die Wélbspannungen o und die Schubspannungen 7 in Dickenrichtung zu | 
Resultierenden zusammenfassen, so entstehen Wo6lbkrafte (Abb. 1). 

41/2 
S= foat. 
—1/2 
und Schubfliisse 
+t/2 
Hee 1B dG 
—1/2 
die Wélbkrafte in Richtung der Stabachse, die Schub- 
fliisse in Richtung der Skelettlinie wirkend. Von den 
Wolbspannungen o werden wir ohne Schaden an- Sport 
nehmen konnen, daB sie iber die Dicke gleichmaBbig Abb. 1. Element tdxds der Profilwandung. 
verteilt sind, so daB wir einfach S—ot schreiben 
kénnen; wir kénnen sogar abschnittsweise gleichférmig verteilte (Langs-) Stringer zulassen, — 
deren Querschnitte f wir uns auf den Zwischenstrecken As ,,umgelegt“’ denken, so da8 wir fir S_ 
etwas allgemeiner schreiben kénnen?: 
Sasol, mitht,=t+f/As. 
Die Schubspannungen t dagegen setzen sich aus zwei Anteilen zusammen, die sich beide als | 
wesentlich erweisen werden, einem Mittelwert tT = T/t und einer Abweichung T von diesem Mittel- 
wert, die von der Normalkoordinate € abhangt und bei den offenen Profilen eine wichtige Rolle 
spielen wird. 

Von den drei Komponenten der Verschiebung interessieren uns nur zwei: die Achsialverschie- 
bung u und die Verschiebung v in Richtung der Tangente der Skelettlinie. Fir die vier unbekann- 
ten Funktionen S, T, u und v von x und s haben wir vier Aus- 
sagen: die Gleichgewichtsbedingung fiir die Krafte in der 
x-Richtung am Wandelement dx ds 


aS aT : : 
FE aE eee (1) 
die beiden Elastizitatsgleichungen 
du du dv 
S=En—, T=¢:(—2 42), (2a, b) 


und die kinematische Aussage der Unveranderlichkeit der Quer- 
schnittsgestalt (Abb. 2) 

<*- =a, cos B-+a,sinB + On, (3) | 
welche die tangentiale Relativverschiebung dv zweier Skelett- 
punkte vom Abstand dx ausdriickt durch die beiden Ver- 
schiebungs- und die Drehungskomponente der Relativbewegung 
zweier starrer Querschnitte; %,, %), # sind Funktionen nur von x, | 
zu deren Bestimmung spater die drei Gleichgewichtsbedingungen | 
fir die Bewegung eines Stabelementes dx in der Querschnitts- | 
Abb. 2. Entstehung der tangentialen chee EIN: pend eu: LFS | 
ae Raat Gee eng Wir eliminieren S, T und v, indem wir (2) und (3) in (1) ein- | 
g z 4 : : : é : 

0 setzen. Man erhalt dann fir u die partielle Differentialgleichung | 

zweiter Ordnung 


3 du au a | 
G ne (« aa) L Et, ie = —G —— t (a, cos B + a sin B+ On). (4). 


Da ihre Koeffizienten von x unabhangig sind, kann man sie durch einen Exponentialansatz mit 


reellem oder imagindrem Argument in eine gewohnliche Differentialgleichung tberfihren. Wir 
wahlen 


O° edz 


aus oj, % und J, 


u(S;<)'=,(s) e-#n* , 


On (S)\ = Oy, en” 4. Oa(A) Onn ean oe, (5) | 
O(x) = On ey, * 
1 Vel. z. B. H. Koller, Jahrb. d. dtsch. Luftfahrtforschung 1940, S. I 852. 
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3 rae rel Uae a ) FAB un=— © t(cun cos 8 + ayn sin B+ Pani). (6) 
E mit E ‘ 
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2 Die Wahl (5) ist dabei bestimmt durch das Ziel unserer Untersuchung: festzustellen, ob es Span- 
-nungszustande gibt, die, hervorgerufen durch eine Belastung des Stabendes «—0, mit wachsen- 
dem x abklingen. 
; Den technisch kaum interessierenden Fall einer mit s beliebig veranderlichen Wandstarke 
wollen wir hier allgemein nicht weiter verfolgen, sondern t und t, hinfort als stiickweise konstant 
ansehen. Dann 1aBt sich fiir (6) die Lésung in geschlossener Form angeben: * 
Un=Ap Cos Ans+ Bn sin Ans— f (on 008 B’ + og, sin B’ + Onri ) cos An(s’—s) ds’, (7a) 
0 
_wobei wir f’,r, statt B.1:, geschrieben haben, wenn diese GréBen von der mit s’ bezeichneten 
Integrationsvariabeln abhangen. Mit 
L(ssx)=—T5($)e45 2: (5’) 
entsprechend (5), liefert das Elastizitatsgesetz (2b) 


Tn= —Gtdn [An sin dn s—Bncosdnst J (a4, 008 8’ + a, 8in B’+ Bar, )sinAn(s’—s)ds’]. (7b) 
0 


Die durch (7a) und (7b) dargestellte Lésung der Gleichung (6) enthalt auBer dem Parameter 
An finf Konstanten 
An, Bn, Kins Xon> On . : (8) 


Ist die Wandstarke nicht langs der ganzen Skelettlinie, sondern nur abschnittsweise konstant, 
so gilt fir jeden solchen Abschnitt eine Lésung der Form (7) mit anderen Werten An, Bn, aber 
mit denselben «,, %,» Jn. Gibt es m verschiedene Abschnitte mit konstantem t, so hat man also 
insgesamt 2m-+3 Konstante. 

Zu ihrer Bestimmung haben wir zwei Arten von Bedingungen: drei Gleichgewichtsaussagen 
fiir die am Stabelement parallel zur Querschnittsebene wirkenden Krafte und Momente und 
je zwei Rand- oder Ubergangsbedingungen fiir jeden der Skelettabschnitte mit konstanter 
Wandstarke. 

Die drei Gleichgewichtsbedingungen fiir das Stabelement sagen aus, daB die beiden Quer- 
krafte Q,, Q, und das Torsionsmoment M7, konstant sind?: 

dQy dQ, dMr 
IRR i hicdwines arn oe 


Da das Stérproblem, das wir untersuchen wollen, dadurch gekennzeichnet ist, daB in den Stab- 
endquerschnitten nur Gleichgewichtsgruppen angreifen, so folgt aus den drei Gleichgewichts- 
bedingungen durch Integration 


Oi ate, Os re. 0 (9) 
Mr=0. (10) 


Die zu (9), (10) hinzutretenden Rand- und Ubergangsbedingungen fir die Skelettabschnitte 
sind zum Teil mechanischer, zum Teil geometrischer Natur. Ist das Profil ein einfacher offener 
“Linienzug, so gilt fiir die beiden Werte des Schubflusses T an den Skelettenden ,,0° und ,,1" 


T(0)=0, | T(s,)=0. (11) 


Ist das Profil ein einfaches Rohr vom Umfang s,, so bestehen die beiden Bedingungen fiir den 


doppelten Zusammenhang: 
u(s;)=u(0), T(s,)= TO). - (12a, b) 


1 Entsprechende Aussagen fiir die drei Normalspannungsresultanten (die Langskraft und die beiden Biege- 
momente) brauchen wir am Stabelement nicht zu formulieren, weil das «-Gleichgewicht fiir das Stabelement 
durch (1) gesichert ist. In der Tat folgt 


dN as dMy — ge dM, jee 
== = = ds=0, = ds=0 
dx Hae sae dx ia dx Pee 
aus (1) zusammen mit (9) und den Rand- oder Ubergangsbedingungen fiir T; bei N unmittelbar, bei My, M, 


durch Teilintegration. Und da die Liangskraft und die Momente an den Stabenden Null sind, verschwinden 
sie iiberall. 
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Man hat also fir die fimf Unbekannten in beiden Fallen finf Gleichungen. Ist die Wandstarke t 


nur abschnittweise konstant, so besteht an jeder Abschnittsgrenze eine Stetigkeitsforderung 
fiir die Verwolbung wu und (nach Abb. 3a) fiir den SchubfluB T. Sitzt an einer StoBstelle eine 
' konzentrierte Langsversteifung, so verlauft wu ebenfalls stetig, wahrend T so springt, wie es die 
Gleichgewichtsbedingung fir das Element dx des Stringers nach Abb. 3b verlangt: 


do,&— » 
Pe) nae EGF PT. 
oder T(k-1) — T,,(*) ss Ew Frit Un(*— 3") ‘ (124 
On 
a a : 
T-Vda 
Cc (Got oe dx) F 


Abb. 3. Ubergangsbedingungen fiir den SchubfluB T: 
a) Sprung in der Wandstarke t, 
b) Langsversteifung vom Querschnitt F, 
c) Verzweigungsstelle. 


Beim offenen Profil mit m Abschnitten gibt es m—1 Abschnittsgrenzen und daher 2(m—1) 
Ubergangsbedingungen, zu denen die beiden Bedingungen (11) hinzutreten. Beim Rohr gibt 
es m Abschnittsgrenzen, also 2m Ubergangsbedingungen, die dann aber die Aussagen (12) tiber 
den doppelten Zusammenhang schon mit enthalten. 

Wenn schlieBlich die Skelettlinie Verzweigungsstellen aufweist, so kommen fiir das in einem 
solehen Punkt abzweigende Skelettstiick ebenfalls zwei Bedingungen hinzu: entweder die zwei 
Stetigkeitsforderungen fiir wu an seinen beiden Enden oder, wenn der Zweig frei endet, (z. B. 
Flanschhalfte eines I-Profils), als zweite die Bedingung T=0 ftir dies Ende. Das aus Abb. 3¢ 
abzulesende Verzweigungsgesetz fiir den SchubfluB 


Thre Umat ole (13’) 


zihlt nicht als neue Bedingung, da es nur die Stetigkeitsforderung fiir T ersetzt, die man auch | 


ohne den Abzweig an der Grenze zweier Abschnitte hatte formulieren miissen. 
Wir haben also immer zur Bestimmung der Unbekannten (8) soviel homogene Gleichungen, 
wie Unbekannte vorhanden sind; diese kénnen nur dann eine von der trivialen verschiedene 


Lésung besitzen, wenn ihre Determinante verschwindet. Das liefert eine transzendente Gleichung | 
fiir die GroBe An mit unendlich vielen Lésungen Ap, durch die wegen y* Et./Gt=22 auch die Ab- | 


klingfaktoren Un bestimmt sind. Zu jeder dieser Lésungen gehért eine besondere Wdolbkraft- 
gruppe S,(s) samt dem zugehorigen SchubfluB T,, und da man zeigen kann, daB dieses System 
von Eigenfunktionen vollstandig ist, so laBt sich jede beliebige Gleichgewichtsgruppe von Wélb- 


spannungen S(s) danach entwickeln. Auf diese Weise ist es grundsatzlich moéglich, die Span- 


nungen in einem Stab zu finden, an dessen Endquerschnitt eine beliebige Gleichgewichtsgruppe 
von Normalspannungen angreift oder dessen Endquerschnitt eine beliebige Verwélbung wu auf- 
gezwungen wird. Die Schubspannungsverteilung im Endquerschnitt des Stabes ist dabei nicht 


| 


willkiirlich vorschreibbar, weil die in beliebiger Verteilung etwa angebrachten Schubspannungen _ 


infolge unserer Voraussetzung undeformierbarer Querschnitte vom Stab (nétigenfalls von einem 
starren Endspant) deformationslos so umgelagert werden, da®B sie zu den Normalspannungen 
passen. 


oN 
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Bis hierher haben wir offene und geschlossene Profile gemeinsam behandelt, und ein Unter- 
_schied zwischen beiden Typen war nur hervorgetreten bei der Formulierung der Bedingungen 
(11) und (12). Es besteht aber noch eine zweite, sehr wesentliche Verschiedenheit. Um zu er- 
_kennen, welcher Art sie ist, gehen wir zweckmaBig von den Ergebnissen der St-Venantschen 
Torsionstheorie aus!. Beim Rohr gibt es dort einen SchubfluB T, der wegen d T/ds=0 kon- | 
stant, aber von Null verschieden ist. Auf die Veranderlichkeit der Schubspannungen tiber die 
~Wanddicke zu achten, ist nicht notwendig. Beim offenen Profil dagegen verschwindet T an den 
_ freien Randern und wegen d T/ds=0 also tiberall. Da kann das Torsionsmoment nur aufgenommen 

werden durch die mit der Dickenkoordinate veranderliche Schubspannung 7; sie ist mit ¢ linear 
und in s, bis auf eine vernachlassigbare Stérung an den Profilenden, konstant. Fiir den Torsions- 
widerstand ergibt sich daher die bekannte Formel 


Ip= = [ ds, (14) 


das Integral tber die gesamte Skelettlinie erstreckt. Mechanisch bedeutet das: Beim Rohr 
-spielt die t® proportionale Plattendrillsteifigkeit der Wandung fiir den Widerstand gegen eine 
Verwindung keine Rolle; beim offenen Profil ist sie, da ein umlaufender FluB nicht entstehen 
kann, wesentlich?. Dieser Unterschied zwischen den beiden Profilarten mu nun auch bei der 

_ Wolbkrafttorsion spiirbar bleiben. Im Rohr treten unter Torsions- wie unter Wélbkraftbelastung 
Schubfliisse auf, die beide von derselben Art sind, und auf die geringfigige Veranderlichkeit der 
Schubspannungen tiber die Wanddicke einzugehen, besteht kein Anlaf. Anders beim offenen 
Profil. Da entstehen Schubfliisse nur infolge der Wolbkrafte. D.h. die Schubfliisse sind klein, und 
man mu, wenn man die Gl. (10) aufstellt, auf den Beitrag des in Richtung der Wanddicke t ver- 
anderlichen Anteils der Schubspannung achten. Selbstverstandlich wird man dabei entsprechend 
‘dem Naherungscharakter einer von der Dinnwandigkeit ausgehenden Theorie annehmen kénnen, 
da8B 7 mit ¢ linear ist und auch konstant in s, da die durch T gebildeten ,,Platten‘‘-Drillmomente 
genau wie bei St-Venant nur in unmittelbarer Nahe der Profilenden gestért werden. Dagegen 
miissen die T bei Anwesenheit einer Wélbstérung wie alle andern Spannungen mit x verdnderlich 
sein, so daf in der Formel fiir den von ihnen getragenen Anteil des Torsionsmomentes 


Mr=GIrd (15a) 


mit I; nach (14), die Verwindung # eine Funktion von x wird. Das Schnittmoment, das eine 
auBere Torsionslast tbertragt, besteht also bei wélbbehinderten offenen Profilen aus zwei An- 
teilen, dem St-Venantschen Anteil (15a) und dem Wolbanteil 


M,=fTrds, (15b) 


die beide von der gleichen GréSenordnung sein kénnen, weil die an den groBen Hebelarmen r; 
(GréBenordnung der Profilabmessungen h, b in Abb. 4) angreifenden Schubspannungen t = T/t 
klein, die an den kleinen Hebelarmen 2 t angreifenden Schubspannungen T gro sind. 

Zu den Querkraften liefert der veranderliche Anteil der Schubspannungen natirlich keinen 
Beitrag. Hier gilt beim offenen wie beim geschlossenen Profil 


0,=0,=f/Tespds, Q=Q=/ Tsing ds, 


und die fiinf Bedingungsgleichungen fiir die Konstanten (8) lauten daher 


> 


beim offenen Profil beim geschlossenen Profil 
T(0) =, u(s,) = (0), 
T(s:) =0, T(s) = T(0). fs 
Jf T cosB ds =0, JS T cos ds =0, 
fetes ds, fT sin Bds =0, 
G1I,0+fTnds=0, [ Trids =0. 


1 Vgl. z. B. K. Marguerre, Bauing. 21 (1940), S. 317. 

2 Wir unterscheiden die aus der Plattentheorie gelaufige Drillsteifigkeit je Langeneinheit einer Platte 
(die zusammen mit den beiden Biegesteifigkeiten die elastische Widerstandsfahigkeit der Platte ausmacht) 
von der Torsionssteifigkeit des Stabes. Beim offenen Profil ist die Torsionssteifigkeit G I das langs s genom- 
mene Integral iiber die Drillsteifigkeit der Wandelemente. Beim Rohr fallt. die Drillsteifigkeit der Wandung 
fir die Torsionssteifigkeit des Stabes nicht ins Gewicht. (Der mechanische Grund ist der, daB das Rohr mit 
Endspanten eine allseits geschlossene Schale ist, die auf diese Belastung nur mit ,,Membran“*-Spannungen 
antwortet, so da& Biege- und. Drillsteifigkeit ruhig verschwindend klein sein kénnen.) 
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3. Offene Profile. In Ziff. 2 haben wir alle mechanischen Uberlegungen angestellt, die wir 
zur Durchfithrung der Wolbkrafttheorie diimnwandiger Stabe brauchen. Wir wollen jetzt in die 
Diskussion des Eigenwertproblems fiir das offene Profil eintreten, und zwar wollen wir, um die 
Gedanken an Bestimmtes zu knipfen, das in Abb. 4 gezeichnete Z-Profil betrachten, das beziig- 
lich des Punktes 0 doppelt antimetrisch und daher einfacher zu tibersehen ist als der allgemeine 
Fall mit fiinf nichtverschwindenden Konstanten (8). 

Der Punkt 0 in Abb. 4 erfahrt keine Verschiebung; es ist daher On = %,=9, und wenn wir | 
die Bogenlange s von 0 aus zahlen, so lautet die Lésung (7): 

Un= An cosdAns + Bn sin dns—On fr; cos An(s’—s) ds’ , (17a) 


0 


ene = —An sin dns + Bn 08 Ans— sal r, sin An(s’—s) ds’. (17b) 
0 


Da die beiden Querkraftbedingungen (l6c,d) aus Symmetrie- 
Sa -griinden erfillt sind, besteht das Gleichungssystem (16) nur aus 
drei Gleichungen ftir die drei Unbekannten Ap, Bn, On. Anstatt 
seine Determinante formal aufzustellen und ,aus ihrem Verschwin- 
den die Gleichung fir die Eigenwerte An zu gewinnen, benutzen 
wir einfacher die beiden Randbedingungen (16a, b), also 

T(+a) = T(—a)=0 mit a=24b 


dazu, um A,, B, durch 39, auszudriicken. Es wird 


—A,sinjna+ Br cosAna—On fri sin An (s’—a) ds’=0, 
0 


+A, sindAna +. Ba cos Ana In fr, sind, (s’+ a)ds’=0. 
0 


b 
Wegen r:(—s) =—n: (s) folgt daraus 
Abb. 4. Z-Profil, autimetrisch zu 0. 9 ro i 
Bi=05 An= ~ —*— [ nsindn(s—a)ds. (18) 
; sin J, a 


Von dem System (16) bleibt nun nur noch eine einzige Gleichung, die Momentengleichung (16e) 


d, a ; a a s 
0 = Cant (= fn sin An(s—a)ds fr sind,s ds— In fir [x sin A, (s’ —s) ds‘ ds) + 
0 0 0 0 1 
+ = GIO, (19) 
poh [eee 


3 


Sie ist in #, homogen und hat eine von Null verschiedene Lésung nur, wenn /, der transzendeten 
Gleichung 


x 1 1 s/a 
An @ Taste j , . 
sin ja J 2 sin da(s—a) d= [7 sin Ans d — — jaa ff sin n(s’—s) d= a> 4 
0 0 0 0 ; i (20) 
t 
edo. 


gentgt. Die Ausrechnung der Integrale in (20) ist in unserem Beispiel sehr einfach, da r, auf 
dem Steg verschwindet und auf den Flanschen den konstanten Wert -—- h/2 hat. Es ergibt sich 


— (35) | adage (1008 And) (cos An(a—b) —c05 dna) — 7,4 Unb—sindad)] +5 (4) =0. ey 


2a A, asiniA, a 3 


Um die kleinste Wurzel dieser Gleichung zu finden, entwickelt man zweckmaBig die tri- 
gonometrischen Funktionen nach Potenzen von A,; man erhalt so 


(Rey [Se ah haar a ay 
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oder, mit 2a=h+2b 


(4,6)? = ( t )’ 4b (h+-2b)? 


b) kh (2h) 
Es gibt also eine Wurzel A, in der Nahe von Null: 
mt 
Ajo = ser 


wobei der Zahlenfaktor m fir ,,mittlere‘* Abmessungen (h~ 26) nur wenig von | verschieden ist. 
Die zugehérige Stérung klingt mit 


enmt x/b2 m Ages | CGiEt. . 
also sehr langsam ab. Mechanisch rihrt das Abklingen her von dem Widerstand, den die Drill- 
steifigkeit der Profilwandung [t’-Glied in (19)] einer elastischen Verwindung des Stabes ent- 
gegensetzt. 

Von ganz anderem Typ sind die héheren Wurzeln. Auf deren Betrag hat das Glied + (t/a)? 
so gut wie keinen Einflu8; denn wenn man in (21) den Ausdruck (2a/h)? + (t/a)?= + (t/h)? in die 
eckige Klammer hineinnimmt, kann man ihn dort zu dem Glied b/a schlagen, das dadurch nur 
um einen verschwindend kleinen Betrag geandert wird. Das heift: die Wurzeln J, der transzen- 
denten Gleichung sind fir n >0 unmittelbar die Nullstellen der eckigen Klammer. Wo diese 
liegen, ist am einfachsten an dem Beispiel a=h=2b zu tibersehen. Da laBt sich die Nullforderung 
fiir die eckige Klammer auf die Form 

il Ab ; ( cos2Ab6 


2 tan 2 cosAb 


) =(Ab—sin db) 


bringen. Die linke Seite ist periodisch und springt bei Ab=s, Ns fs ee 3m, oF .. ZWi- 
schen +0 und —o; die rechte schwankt mit einer Abweichung von +1 um die Gerade —Ab. 
Die Treffpunkte 2,b liegen daher in der Nahe der Sprungstellen und nahern sich ihnen fir 


wachsende n immer mehr. Die erste Wurzel liegt bei 4,b=2,48, und fiir h=2b6 wird daher 


Ay _ V3.2 + 

A t2 Ab 
fiir ein dinnwandiges Profil klingen also die héheren Eigenkraftgruppen um eine GroBenordnung 
schneller ab als das zu A, gehérige Hauptglied der Stérung. (Die héheren Kraftgruppen sind, unab- 
hangig von der Wandstarke — also auch fiir t—»0 —, Wolbkraftgruppen in dem in der Einleitung 
definierten Sinne, die niedrigste ist es nur bei Vorhandensein einer Drillsteifigkeit der Profil- 
wandung.) 

Was wir hier durch die Zahlenrechnung an dem Beispiel des Z-Profils mit h=2b gefunden 
haben, gilt fir offene Profile allgemein: Die zur kleinsten Wurzel gehoérige Stérform klingt 
langsam ab, die anderen im Verhiltnis nb/t schneller, wobei n der GroBenordnung nach etwa die 
Reihe der ganzen Zahlen durchlauft. Vom technischen Standpunkt aus ist also allein die niedrigste 
Wurzel der transzendenten Eigenwertgleichung interessant; fiir diese kann man aber eine Formel 
angeben, die allgemeingiiltig und zugleich besonders tbersichtlich ist. Wir wollen das am Bei- 
spiel des punktsymmetrischen, aber sonst beliebigen Profils zeigen, indem wir in Gleichung (19), 
in der rm noch irgendeine Funktion von s sein kann, den Grenziibergang A,—>0 vollziehen. Wenn 
wir G@, herausheben und A,=A,—>0 gehen lassen, wird aus (19) 


t 
—= (Oho, b? 


al f n fnle—s ete ees ts: rvs =<. (20) 
0 0 0 é 


Die eckige Klammer kénnen wir durch Teilintegration umformen. Setzen wir zur Abkirzung 
s 


ie ds'=—q(s) , y(s’)=9", 
0 
so liefert zweimalige Teilintegration des ersten Terms 


[offre alas f n{[- p-(%—s)| + [vasa -f n fara 


[afer foun fees jore 
One 0 0 0 
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Im zweiten Term wird jeder Faktor einmal der Teilintegration unterworfen: 


a 


fig weeded : + {[-e—al + fas} [ws] + f as} 
0 0 : : 
= [vas| oe +h [vast 
fore (foay 
[ortas 


Beides zusammen ergibt also 


Dafiir kann man schreiben 


mit i 
1 
P*=P—Yo> go= | vas, (22)) 
: 0 
und aus (20’) wird daher I 
| [ee eee (22’) | 
24 fords 


Die GroBe p* ist die spezielle Torsionsverwolbung des Querschnitts, die zu einer Verwindung j=] 


; 1 é 
gehért und deren Mittelwert 1 fords verschwindet. In dem von uns betrachteten Sonderfall 
0 


; , 4 2 Pas 10 ‘ | 
des punktsymmetrischen Profils konstanter Wandstirke ist das zugleich die ,,nattrliche Ver- 


wolbung, die bekanntlich durch die drei Bedingungen 
$1 $1 $1 
[otar=o, [otzaF=0, | otyaF=0. (23’) 
0 0 0 
mit dF'=tx,ds unter. den co? Verwélbungen des tordierten Profils ausgezeichnet ist. 
Indem man schon in der Differentialgleichung (6) An=0 setzt, was damit gleichbedeutend ist, 
an Stelle von (2b) die Aussage 


Sot pee 0 (24) 


zu benutzen, d. h. die Formanderungswirkung des Schubflusses T zu vernachlassigen, kann man 


zeigen, daB sich fiir den Abklinggrad Mo=Ao | Gt /Et, der niedrigsten Eigenfunktion ganz alle 
gemein der Ausdruck 


GI 
oes (25) 
ergibt, mit ; 
Cu = [et tds, ee (25’) | 
0 ! 


wobei tibrigens t, mit s durchaus veranderlich sein kann'. Die zugehorige Verwélbung u, und. 
damit auch die Wélbkraft S=Et,du/dx, verteilt sich wie die Wélbfunktion g* der St-Venant- | 
schen Torsion, und das Ergebnis unserer Rechnung ist also, das man jede andere Wolbstorung 
zerlegen kann in einen y*-Anteil und einen Rest, der aber so schnell abklingt, da® in kleinem 
Abstand von der Stérstelle nur der g*-Anteil bleibt. J 

Damit ist die Frage beantwortet, wie ein Stab mit offenem Profil auf eine beliebige Wolb- | 
stérung reagiert. 


4, Rohre. Auch fiir die Behandlung eines durch Wolbkrafte beanspruchten Rohrprofils | 
haben wir in Ziff. 2 schon alle mechanischen Aussagen zusammengestellt. Die rechnerische Ver-| 
folgung der zu den homogenen Gleichungen (16) gehérigen fiinfreihigen Determinante, aus der| 
sich die Eigenwerte bestimmen, ist auch hier sehr mithsam, so da® wir uns fir die weitere Durch-| 


* R. Kappus, a, a. O., Formel (18) oder auch K. Marguerre, Ringbuch der Luftfahrttechnik II A 9, S. 14/15. | 
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_fiihrung der Theorie analoge Beschrankungen auferlegen wie beim offenen Profil, damit nicht 


. das mechanisch und mathematisch Wesentliche durch ein zu umfangreiches Formelwerk ver- 
: deckt wird. Wir wollen voraussetzen, daB das Rohrprofil zu den beiden Achsen y und zsymmetrisch 
sei. Dann wird wieder %,=«%.—0, die Querkraftbedingungen (16c, d) werden trivial und scheiden 
aus der Untersuchung aus. Die verbleibenden drei Gleichungen werden so einfach, da8 man 


aN 


wieder A,, By mit Hilfe von (16a, b) durch 9, ausdriicken kann und aus (16e) die Gleichung 


fir An unmittelbar erhalt. Legt man noch den Nullpunkt der Bogenlange in einen der vier Sym- 


metriepunkte des Profilskeletts (beim Rechteck z. B. in eine der Seitenmitten) und zahlt die 
\ Bogenlinge, beginnend beim gegeniiberliegenden Symmetriepunkt, von —s,/2 bis +s,/2, so 
erhalten die Zusammenhangsbedingungen (16a, b) die Form 


u(+s,/2)=u(—s,/2), T(+s,/2)=T(—s,/2) , 


| und fiir die Konstanten 4,, By des Ansatzes’ (7) ergibt sich 


$,/2 


a Sk On s 
A,=0, Bu = sayy | 11008 In (s——) ds. (26a, b) 
0 


. Natiirlich kann man A,=0 auch daraus folgern, daB u als eine doppelt antimetrische Funktion 
in den Symmetriepunkten verschwinden mu8, und der Ausdruck (26) 1aBét sich einfacher schrei- 


ben, indem man ihn aus der Forderung 


5,/4 
dun Sy 
[Gods =m (4) =0 (27) 
0 
herleitet: id 
wk On . SH ‘ 
Bn = saga | 1 0084An(s— 2) ds. (26b’) 
0 


DaB die Formeln (26b) und (26 b’) gleichbedeutend sind, erkeunt man, wenn man die aus der 


- Querschnittssymmetrie folgenden Eigenschaften der periodischen Funktion r(s) beriicksichtigt. 


Die beiden Formeln gelten natiirlich nur, wenn t nirgends springt, so daB der eine Ausdruck 
(7a) fir die ganze Skelettlinie besteht. Beim Rechteck mit Eckgurten z. B., das wir in Ziff. 6 


 diskutieren werden, mite man du,/ds in (27) nach (2) und (3) ersetzen durch T,/Gt—m@n, und 
_darin T, mit Hilfe der Sprungbedingung (12’) bestimmen. 


van Stelle von un, Tn die Funktionen 


Indem man (26) in (7b) und das in (16e) einftihrt, erhalt man (da¥#,, + 0 sein soll) die gesuchte 
transzendente Gleichung fiir die unendlich vielen Eigenwerte An: 
s,/4 34/4 4/4 s 
1 s Risto i / 
Se Ai for COs An (s—) ds fin cosdnsds— [ rt lee sinAn(s’—s)ds’ds=0. (28) 


0 0 0 


_ Die Aufgabe, fir ein gegebenes Rohr die Wélbkraftgruppen zu bestimmen, ist damit grund- 
satzlich gelést. Man braucht nur fiir den aus der Querschnittsform folgenden Verlauf von ri(s) 
die Integrale in (28) auszufiihren und die dadurch entstehende transzendente Gleichung fiir An 
zu lésen. Dann geben (7a, b) sofort die Verwélbung un und den SchubfluB T,. Die zugehérige 
Wolbkraft S, folgt aus (2a). Da (28) eine transzendente Gleichung ist und daher unendlich 
viele Lésungen /, hat, so gibt es je einen Satz von unendlich vielen Funktionen u,(s), Tn(s), 
Sn(s). Sie sind die Eigenfunktionen der Differentialgleichung (6) (oder zweier anderer, die wir 
statt dessen hatten aufstellen und der Behandlung unseres Problems zugrunde legen kénnen). 
Es ist zweckmaBig, sie von den Parametern zu befreien, die noch in ihnen enthalten sind, und 


a Un — Tn 
an she 


als Eigenfunktionen unseres Problems zu benutzen. gn und yn hangen nur noch von Form und 


- GroBe des Querschnitts ab, nicht mehr vom Material des Stabes und der Intensitat irgendeiner 


Beanspruchung. 
Es laBt sich allgemein zeigen, daB die Eigenfunktionen eines jeden Systems untereinander 


orthogonal sind, d. h., dab 


fpnlsypmls) ds=0, ——f yals) pn(s) ds = 0 (29) 
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ist, wenn gn und gm und ebenso wp und Wm zwei zu verschiedenen Werten A, und Am gehorende 
Funktionen sind. Um das zu bestatigen, braucht man nur die beiden Funktionen nach (7a) 
und (7b) einzusetzen und bei der Ausfithrung der Integrale die Gleichung (28) und die schon 


erwaihnten Symmetriebeziehungen in geeigneter Weise zu berticksichtigen. 


Wir kénnen die so gefundene Lésung benutzen, um die Spannungen und Formanderungen zu 
berechnen, die in einem Rohr entstehen, wenn man seinem Endquerschnitt x=0 eine beliebige 


Verwélbung u(s) aufzwingt oder in ihm eine beliebige Gleichgewichtsgruppe von Normalkraften | 
S(s) angreifen lat. — Man kann also z. B. einen Stab untersuchen, der ein Torsionsmoment Mp | 


ibertragt und dessen Verwélbung im Endquerschnitt durch eine Einspannung verhindert ist. 
Ohne diese Einspannung wiirde ein Spannungszustand nach der Bredtschen Theorie entstehen, 
zu dem die Verwélbung! 


ua(s) = Pols) = Gye Hs) 


gehért mit der ,,Wélbfunktion‘ (Verwélbung je Einheit der Verwindung) 


s Sy 


os) =— [rid +> [rds 


0 0 


Um diese Verwélbung im Endquerschnitt zu beseitigen, iberlagern wir eine Summe der eben 


aufgestellten Lésungen 
u(s,x) = |7(s) + me Cn Pn() etn | : 
Fur den Querschnitt «=0 ist dann 


1 (3,0) = Bo p(s) + > Cngpals)] 


und das soll identisch Null sein. Die Koeffizienten C, sind also die Entwicklungskoeffizienten | 
von —qg(s) nach den ,,Eigenverwélbungen* gy, und kénnen, da diese ein Orthogonalsystem bil- | 


den, berechnet werden als Quotienten zweier Integrale 


S As) Gnls) ds 
S nes) ds 


wobei diese Integrale grundsatzlich ther den ganzen Umfang zu erstrecken sind, aber wegen 
der hier angenommenen doppelten Symmetrie des Querschnitts auch nur itiber die Halfte oder 
ein Viertel des Umfangs erstreckt zu werden brauchen, beide natirlich in gleicher Weise. Die 
endgiiltige Lésung schreibt sich dann in der Form , 


C= 


(30) 


U(S3%) =, | vs) + ys Gree mine pals) | 5 
; 


S(s, x) ies Et, 34 eo Ch [en Cmte Pn(s) > 
1 


Wa ff T(s, x) = T+ Gye >) Cn e-H'n® nls), 


B(x) = By[1 +) Cre-tns] mit B = Mr 
1 7 


GIy’ 


Abb. 5. Rechteck a+b als Rohrquerschnitt. 


SchubfluB ist. 


5. Das Rechteckrohr. Wir wollén an einem Beispiel zeigen, wie die Eigenfunktionen aus- 
sehen, und wahlen dazu ein Rohr mit rechteckigem Querschnitt, zundchst ohne Eckgurte (Abb. 5). 
Der erste Schritt ist die Berechnung der Eigenwerte j,. Dazu dient grundsatzlich Gleichung (28); 
die Rechnung wird aber in diesem einfachen Beispiel tbersichtlicher, wenn wir T, nach (7b) 
und (26b’) zunachst explizit schreiben und erst dann die in (28) geforderte Integration zur Bil- 


1 Siehe FuBnote 1 von Seite 27. 


(31) 


wobei JT) der nach der Bredtschen Formel zu bereelnenas | 
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dung von My vornehmen. Wenn man den Null 


legt, folgt aus (26b’): 
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punkt der s-Zahlung in die Mitte der Langseite a 


An ed OB at+by]42 af. a+b \}(a+5)/2 
c= iz [sin An (s— )] ees [sina ba, “| | 
4 ain a+b | 2 2 0 2 i ( 2 / Ja/2 
4 2 eLAgb 
5 Oa fib a= ae 
ian 2 v1 2 a+b z 
sin An —— 


Die Ausfithrung des Integrals in (7b) im Bereich 0 Ss <a/2 liefert 


; Se 5 (1+ cos dns); 
insgesamt gibt das fir die Langseite 
ig . 
sin - 
b a—b 2 
Pra= Git?,, (: woes mg ona) : (32a) 
sin j,, —— 


Die Rechnung wird am iibersichtlichsten, wenn wir die Bogenlange nicht tiber die Ecke hinweg 
weiterzahlen, sondern fiir die Schmalseite 6 eine neue, von deren Mitte aus rechnende Zahlung 
_einftihren, Der Ausdruck fir T, in diesem Bereiche ergibt sich dann aus (32a) durch einfache 


Vertauschung von a und b: 


0-4 


2 


Tat =GtO, (: _ 


(32b) 


Die Bildung von M7/4 erfordert nun nur die Multiplikation der Ausdriicke (32) mit den Hebel- 


armen 6/2 und a/2 und die Integration von 0 
so als A,-Gleichung 


bis zu den Ecken s=a/2 und s=b/2. Man findet 


fn ae nb sin sin Ana 
a/by\2 ee Megan Wes 2 (3) oe Pe aes 
2 (3) ce PNT) a+b an 2 D Te aye : a+b 
< nA,—x— sin A, 
GD) ‘ 2 
oder 
Ana: And 
ab jy 2 aby ae ies) 
ah ee a+b 
py sb 
oder auch he 5 are 
Ane ee ee 33 
CL Sap poms Sirs al cal arg Won ey 


Sie hat z. B. fir a=2b folgende Lésungen: 


a+b 


An > 


= 3,022,- 


= In F- = 173,16°, 356,58°, 718,29°, 898,64°, 1259,03°, 
6,223, 


12,54 15,68, 21,97. 


Die Schubspannungen sind durch (32) gegeben; fiir die Verwolbung up ergibt sich entsprechend 


b 
; b sin An > 4 = i . 
a— x b ee a 
Un = a 21, rae sin Ans fir. —— Ss S45 (34a) 
sin A, ——— 
2 
ie 
sin == 
lee "2 i i] PDS fe <a 5b 34b 
Un 2A, 5 mere sindns fur F=ss ete a (34b) 
sin A, —~— 


34 Fliigge und Marguerre: Wélbkrafte in diinnwandigen Profilstaben. Ingenieur-Archiv 
Sa AI els SMa crac a kaw) Sans 8 i 


Die beiden ersten EKigen- 
verwolbungen un/#, sind 
in Abb. 6 aufgezeichnet. 
Vor allem fiir die groBe- 
ren n unterscheiden sie 
sich nur wenig von ein- 
fachen Sinuslinien. 


Abb, 6. Eigenverw6lbungen 
des Rechteckrohres: 


a) die erste Eigenverwélbung ¢,, 


b) die zweite Eigenyerwélbung ¢@p. 


Wenn wir die Span- 
nungen in einem Rohr 
berechnen wollen, das 
durch ein Torsionsmo- 
ment Mr belastet und 
dessen Endquerschnitt 
wolbstarr eingespannt 
ist; so mlissen wir die 
St-Venantsche Verwdl- 


bung 
ee) fir die Langseite a, 
2 (a+b) (35) 
»  a(b—a) ‘ : , 
ree s fir die Schmalseite b, 
nach den Eigenfunktionen (34) entwickeln. Wir finden fiir die Integrale in (30) 
(a+)/2 
(a—b)? mp GEA 
* re n 
YP gnds Ban? Ana feb ee aaeb } 
0 : Se -F etg 5} 
was wegen (33) gleichwertig ist mit 
(a-+b)/2 
*p ds SS ale, GOR Nar. 36 
ae 2 An? (a+b) (36a) _ 


0 


ay 
ig 


: und entsprechend 


Z 


7a 
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(a+8)/2 . 9 And Roel @ 
i, Sy ene ab(a+b) (a—b)? asin” — + 6 sin? ; 
J Pn 16 An? I 16 Ane A r a+b ° (36 b) 


‘ 2 
Nachdem man die C, gefunden hat, ergeben sich die Spannungen und Verformungen aus (31). 


Der Verlauf des Schubflusses T ist in Abb. 7 fur einige Querschnitte aufgezeichnet. Man 
sieht, wie bald er die elementare Bredtsche Verteilung T’'=const annimmt. Technisch be- 
deutet das, da®B die Wélbbehinderung im Endquerschnitt nur einen 6rtlich sehr be- 
grenzten EinfluB hat 


und fiir die Deformation 
im groBen unwesentlich 
ist, sehr im Gegensatz 
zu dem Verhalten der 


- spannungen sind natiir- 
ihre GréBe gibt die hier 


~ kunft. 


| spiel von Ziff. 5 hatte sich 


Stabe mit offenem Pro- 


fil. Die 6rtlichen W6lb- 


lich betrachtlich. Uber 


dargestellte Theorie Aus- 


6. Rechteckrohr mit 
Eckgurten. In dem Bei- 


gezeigt, daB die Eigen- 
werte /, zueinander et- 
wa im Verhdltnis der 
ganzen Zahlen standen, 
daB alsobeim Rohr — im 
Gegensatz zum offenen 
Profil =o der niedrigste Abb. 7, Schubfiu8 T (s, x) an einem wolbfest eingespannten Rechteckrohr (a/b = 2). 
EKigenwert keine Sonder- Abklingen der ,,Stérung** von links nach rechts, 

rolle spielte, und dem 

entsprach das schnelle Abklingen der ganzen Stérung. Auch beim geschlossenen Profil kann aber 


eine Auszeichnung des niedrigsten Eigenwertes (oder der niedrigsten Eigenwerte, bei mehr als 


4 Gurten) und dementsprechend ein langsames Abklingen der zugehérigen Stérung eintreten, 
dann namlich, wenn das Rohr kraftige Langsgurte aufweist. Denn der in den Gurten konzentrierte 
Teil der Wélbkrafte wird durch die Rohrwandung natirlich langsamer abgebaut als die in der 
Wand selbst wirkenden Langsspannungen, und das Abklingproblem kann damit technisch sehr 
an Bedeutung gewinnen. Genau wie beim offenen Profil bringt es aber die Auszeichnung des 
niedrigsten Eigenwertes mit sich, dafs man dann das Wesentliche des Vorgangs durch eine ele- 
mentare Theorie erfassen kann; diese ist als Theorie des Kastentragers schon lange bekannt, 
und unsere Aufgabe besteht nun darin, den Zusammenhang mit dieser Theorie herzustellen. 

Wir beschranken uns zweckmaBig wieder auf ein einfaches Beispiel und wahien den doppelt- 
symmetrischen Rechteckkasten mit vier gleichen Eckversteifungen. Die Gleichungen (7a) und 
(7b) gelten jetzt fiir die vier Seiten einzeln mit Ana, Bna fiir die Langseiten a und Anp, Bno fiir 
die Schmalseiten b. Die Rechnung wird auch hier iibersichtlicher, wenn wir fir jede Seite eine 


‘besondere Zahlung der Bogenlange s einfihren, die ihren Nullpunkt in der Seitenmitte hat. 


An den Ecken sind 4-2 Ubergangsbedingungen von den Typen (12a) und (12’) zu erfillen, 


Da es keine wesentliche Mehrarbeit macht, die Wanddicken von Lang- und Schmalseiten in 


(12’) als verschieden vorauszusetzen, wollen wir das tun. Bezeichnen wir sie mit ta und ft), den 
Gurtquerschnitt mit F, so heiBen die beiden Bedingungen fiir eine der Ecken 


tine 40/2) = una(—b/2) . (31) 
und, wegen t,=t, also Ew, =GA;, 


Tna(+a/2) ~ Tno(—B/2) = 5 CA2 F [una( + a/2) + uns(—b/2)] . (38) 
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Da die aus der doppelten Symmetrie des Querschnitts folgende Bedingung Ana=Ans=9 auch)| 
hier gilt, so sind, nach Einfiihrung des Losungsansatzes (7) (getrennt fiir jede Rechteckseite),,| 
(37) und (38) zwei Gleichungen zur Berechnung von Bna und By: 


a/2 b/2 |) 
Alay hex) O,b a pet tA Oa ee ! 
Bra sin yee e J 0s an (s— 5) ds = — By» sin + 5 J 008 an (s x) ds, | 
0 0 
a Pas Oe per maar b | 
n@ n : a n na . ie ne 
i | Bue COSia a oe sin An (s—$) ds| = iif | Bas C09 ea agi [sin Ap (s x) ds| 
0 0 
= ae [una (4-4/2) + uno (—b/2)] = = | Bre sin ee 


a/2 b/2 


2S ae J cosdn (s— +) ds + ae [05 n (s—3) as]. 
: 0 


Ausfihrung der Integrale und Umordnung ergibt 


Ana Peay Lien Aes a Anb AnF 3 An b\ 
Bed («. cos —5 3 Sin ) Bay (« cos —S—— 3 sin — | )= 
al Oo, [.tha—tab b An @ Teel il ve Axa a Anb Vee Irae gee il) 
oie [a + F (ta cos 2 5) sin 5 ) 2 (1s cos 5 2 sin 5} ) 
Sete Ce AD OVO ete eA Ue ALD 
Bnasin 5 + Bn sin 5 i | 5 sin 5 T-5 Sin 5 le 
Mit 
POs cer ARO TEM ie een ea WEA Wine AD, 
A = t, cos 5 Sin —5— + te cos Pee) — A, F sin Se tyke (39) | 
folgt daraus 
Oo, 7 6b tha—t,b Anb 
Bra tat 2A | | 
panes O, [a 4 tab—tha . dna | Se 
cia yi ie pe IES ae 
Die Schubfliisse werden daher nach (7b) 
b ted 10s) 
ae | 
Tra = COnta (5 | A sin — 08 Ans) 5 a 
1 
he a TAD tyme enG 
Tx» = COnte ($ se a sin —5 05 Ans) 
und die Momentenintegration (l6e) fuhrt auf die Eigenwertgleichung 
b 2a Oe rt a-<t Neco an nau tO. a\2 6b Det hO = ty a ae ane ieee 
ts ie IU Nance EY Pane ae eo | + |(4) baby Mae tae on Gennes |=0 
oder ‘ ( by? . or : 
a th a—t, ‘ nde: is, 
t a ‘ 
‘ (tp ata b) lee sin —)— sin 5 
oder schlieBlich 
Ana Anb 2 (tp a—t, b)* 
tac = 
ote Tart tet bicsa Va thea Gc meee bye, Ota (25) 


Fur F=0, ta=t=t kamen wir zuriick auf die alte Gleichung (33). Fir F +0 interessiert uns 
vor allem die erste Wurzel J, der Gleichung (42). Ist F hinreichend gro8, so erhalten wir sie aus 
einer Reihenentwicklung der Cotangens-Funktionen: 


Dige 2th A, 4 A,b 2(t a—t,b)? 
aot - 
aie OU ESO Ay ab (ty a-+t, b) Aye. 
Daraus folgt 
8 
Nita 43 
i a Lange (43) 
Gogh ae Rikcn ta 
TG (aa-+nd)] (2 a 4) 


Der Summand = (taa+t)b) neben F ist eine Art mittragende Querschnittsflache der Rohrwandung. 


ek EMaked f teh ec} ree ars ay : 
Aah, pare ay i SS ie es : hy wea, 
1 bkrafte in diinnwandigen Profilstaben. SG Uae, lif 


lem ten Chiedsder etg-Entwicklungen in (42) und zeigt, in welcher Richtung 
sich andert, wenn at/F von Null verschiedene Werte annimmt. Mit 


se 


a 


ee i 2 foictas ©" Fo : (taa+tb) = F, 
(effektiver Querschnitt) lautet die Formel fiir den Abklinggrad 
Betis ale tg FC. : 8G 1 j 
5 ee My ee BRA G bee’ / 43 ) 


i Bo Ve 
tS ti eine aus der Kastentragertheorie wohlbekannte Beziehung, die sich dort ergibt, wenn 
jede der vier Wande als einen Trager mit hohem Steg behandelt, der lings der Eckgurte 


den Nachbarn zusammenhingt und eine Biege- und Schubverformung erleidet, die die 
nzeltragerquerschnitte eben laBt. 


Aylarh) Ao Sy 


She Ce is Sater ae 
reas ie Si Re fo -  Asymptoten: +— n-Vie 
. a Py A 
= 3 
4 a 
e 
a 
b i 
“a 
7 
Ee (arbjt 
d | [Sana 
0 | eens 2 3 


Abb. 8. Der niedrigste Eigenwert 7, in Abhingigkeit yom Flachenverhiltnis 
Wandung:Gurt. Seitenyerhiltnis b/a als Parameter. ‘ 


. 


_ Wir verfolgen die Kastentragertheorie! im Rahmen dieses Aufsatzes nicht weiter, sondern 
geben nur noch das Diagramm Abb. 8, das zeigen soll, bis zu welchem Querschnittsverhaltnis 
“Wand zu Gurt der niedrigste Kigenwert nach (43) berechnet werden darf. Wir beschraénken uns 
dabei auf den Fall ta=ts, wo (42) mit den Abkirzungen 


E 3 =|, y= BG, cH? (44) 
“Die : Hed I OSA 
| Ci acon CAM ater ee ane we 
aw ah (45) 
Uf) og ctg nee = ae 


9 of Q f 3 of VES A aay 
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annimmt, sich also darstellen laBt als eine von dem Parameter ¢ abhangige Schar von Kurven 
n(€). Die Grenzkurven dieser Schar sind 


&? a= + tg a fiir €=1 (Quadrat) , (457) : 


und 


&?, = oe a 2 fir €=0 oder €= (sehr flache Rechtecke). (45/9) 
1 Cle 7) ‘ | 


Die Formel (43) geht mit den Abkirzungen (44) tber in 
12) , 

n= —— . (46)) 

\3+2 | 

die von €=b/a unabhangige, fiir alle Kurven (45) gemeinsame Schmiegungskurve in der Um-; 


gebung des Nullpunktes (7=2€ ist die Wulleubletseansonrel: Man erkennt, daB, solange der 
Gurtquerschnitt nicht kleiner ist als der der Wande (1), die Formel (43) eine sehr gute: 


Naherung abgibt. Fur = (a+b)t > F wird sie bald unbrauchbar, und zwar aus zwei Griinden:: 


erstens wird die Formel. zahlenmaBig schlecht, weil sich die‘Kurven (45), vor allem bei flachen! 
Rechteckkasten, von der Grenzkurve mehr und mehr entfernen; zweitens ist es bei den hohen} 
Abklinggraden, die zu Flachenverhiltnissen €? > 1 gehéren, unzulassig, die niedrigste Wurzel! 
der Gleichung (42) allein als wesentlich gelten zu lassen: Die nachst héhere liegt im Falle b=a) 
z. B. bei 


1 = No = 6,85 ie ade (45)], i 


hat also knapp den vierfachen Betrag von ny =1,72, und da kann man die Abklinggrade nicht: 
mehr als der GréBenordnung nach verschieden ansehen. Technisch verliert glicklicherweise das; 
ganze Problem gleichzeitig viel an Bedeutung: Die gesamte Wolbstérung klingt, wie Abb. 7 
an dem Zahlenbeispiel a=2b, F=0 zeigt, so schnell ab, daB ihre eee meist nicht mh 
notig ist. 

Fur Rohre, deren Querschnitt gegen eine Deformation in seiner Ebene hinreichend geschiitzt 
ist (sei es durch hinreichende Dicke der Wandung, sei es durch viele starre Spante), gilt also 
das folgende Ergebnis: Wenn Langsgurte ganz fehlen oder schwach sind, klingen Wélbstérungen 
sehr schnell ab, haben also technisch meist keine Bedeutung. Will man sie bestimmen, so kann 
das nur mit Hilfe der hier gegebenen Theorie geschehen, die auf unendlich viele Eigenkraft- 
gruppen fihrt. Sind kraftige Langsgurte vorhanden, so klingt eine bestimmte Wélbstérungs- 
gruppe langsam ab. Zu ihrer Untersuchung, die technisch von groBer Wichtigkeit sein kann. 
reichen die Kastentragerformeln aus!. 


1 Siehe FuBnote 1 von Seite 37. 


(Eingegangen am 24. Dezember 1948.) 
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Der Wirkungserad der Planetengetriebe. 
oa . Von R. Poppinga. 


x 1. Einleitung. Wenn man sich bei der Lésung einer getriebetéchnischen Aufgabe aus be- 
stimmten Grinden fiir ein Planetengetriebe entscheiden méchte, so interessiert vielfach die Frage, 
wie hoch sein voraussichtlicher Wirkungsgrad im Vergleich zu dem des in der Regel verwendeten 
Standgetriebes sein wird. In der nachfolgenden Untersuchung wird an einer Reihe von Bei- 
spielen gezelgt, wie man derartige Vergleichszahlen wherschlagig yorausberechnen kann. 


ww 
i) 


WW 


Zz Abb. 1. Schematische Darstellung der behandelten Planetengetriebe, 


st = Steg, Planetenradtriger. 


_ Es werden insgesamt 17 riickkehrende (gleichachsige) Planetengetriebe behandelt, die in 
Abb. 1 schematisch wiedergegeben sind. Die Typen I bis X stellen die grundlegenden Bauformen 
lrei- und vierradderiger Planetengetriebe dar und XI bis XVII haufig verwendete zusammen- 
yesetzte Planetengetriebe. Alle Getriebe haben je ein An- und ein Abtriebsglied, sind also reine 
Ubersetzungsgetriebe im Gegensatz zu den Verzweigungsgetrieben (Summengetriebe oder 
Jifferentiale), auf die hier nicht eingegangen wird. 


_ 2. Wirkungsgrad des Standgetriebes. Bei der Berechnung des Planetengetriebe-Wirkungs- 
rrades wird davon ausgegangen, da dieser durch den Wirkungsgrad des Standgetriebes aus- 


‘edriickt werden kann, der sich fiir jede Verzahnung nach folgender Beziehung errechnen laBt:! 2 


It 
ess 


Hi t—v = 1—\ Ewen sat (1) 


4 2% 


fe M. ten Bosch, Vorlesungen tiber Maschinenelemente. Berlin 1940. 
2. R. Poppinga, Stirnrad-Planetengetriebe. Stuttgart 1949. 
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Verlustgrad Vv 


Ubersetzung -t 
Urbs ipo AZA 4 6 8 70 


Abb. 2 und 3. Der Verlustgrad yon AuSen- und Innenverzahnungen in 
Abhiangigkeit von der Ubersetzung bei verschiedenen Ritzelzahnezahlen. 


Annahmen: Zahnform- und Teilungsfebler vernachlassigt. Zahnreibungs- 


ziffer 4=0,05. Evolventenverzahnung ohne Korrektur, 4= 20°. Kopf- 
kiirzung der Innenzahnrader. 


Verlustgrad’ Vv 


c) Bei gleichem i und z, ist der Verlustgrad einer SMUG CISL groBer als der einer 


16 7 
% 
oo 
42 
08 
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, Uberseteung 40 
Zz ¥ 6 moe, 70 
Innenverzahnung. 


nH = Nav] Nan = Verzahnungswirkungs- 
grad = Verhaltnis der Abtriebs- 
(Nutz-) leistung zur Antriebs- 
leistung, a 
ae eae Verlustgrad = Verhalt- 
der Verlustleistung zur An-} 
risbsleeenne! 
& = Zahnform- und Téilfehlerfaktor; 
yw = Zahnreibungsziffer, 


ip Uberdeckungsgrad, 


bes ck Ryley = Ubersetzungsverhaltnis, 

z, = Zahnezahl des aufenverzahnte 
Rades (kleines Rad, Ritzel), 

z. = Zahnezahl des innenverzahnte 


oder auBenverzahnten Rade 


(groBes Rad), 


das —-Zeichen gilt fir die Au®enver- 
zahnung (Gegenlaufigkeit), i 
das +-Zeichen gilt fur die Innenver- 
zahnung (Gleichlaufigkeit). 


Wirkungsgrad laBt die Lagerverlust 
auBer acht, die nach Budnick! bei heu- 
tigen Getrieben mit Walzlagern, gute 
Schmierung und geringen Zapfendurch 
messern gegentiber den Verzahnungsver 
lusten auch nicht ins Gewicht fallen und 
daher vernachlassigt werden kénnen. 


Um den Einflu8 der Ubersetzung i, 
der Zahnezahlen und der Verzahnungs- 
art (AuBen- oder Innenverzahnung) au 
den Wirkungsgrad zu zeigen, ist in den! 
Abb. 2. und 3 der Verlustgrad y als Funk-,; 
tion von i mit z, als Parameter bei ver-; 
einfachten Annahmen (&=1, bu = 0,05) 
dargestellt worden, und zwar je fir) 
AuBen- und Innenverzahnung. Die er-- 
haltenen beiden Hyperbelscharen gestat-} 
ten folgende Aussagen: | 


| 


a) Allgemein fallt der Verlustgrad 
beider Verzahnungsarten mit wachsender 
Zahnezahl ab. Eine Erhéhung der Ritzel-| 
zahnezahl im Bereich z,=10 bis 20 fahrt) 
besonders bei niederer Ubersetzung zu 
einer starken Abnahme des Verlustgrades. 

b) Bei gleichbleibendem z, fallt der? 
Verlustgrad einer AuSenverzahnung mit| 
steigender Ubersetzung ab, wahrend er| 
bei der Innenverzahnung zunimmt. 


1 A. Budnick, Zeichnerische Behandlung von Kraften und Momenten in Koppel- und Radertrieben. 


VDI.-Forsch.-Heft 388 (1938). 


/ 


ey = oN 
_ XVIII. Band 1950. 
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: d) Bei i-Werten etwas iiber 1 weisen Innenverzahnungen sehr geringe Zahnreibungsverluste 
auf. Der Verlustgrad ist hier um ein Vielfaches geringer als der von AuBenverzahnungen. 


7 Neuzeitliche Getriebe mit auBenverzahnten Radern erreichen mittlere Verlustgrade von 
9,5 bis 2% fir jeden Zahneingriffl. Uber Innenverzahnungen sind bisher keine Ergebnisse 
_veroffentlicht worden. DaB sie aber bei gleichen sonstigen Bedingungen wesentlich geringere 
_Verluste aufweisen miissen als AuBenverzahnungen, mag folgendes Beispiel zeigen. Wahlt man 
~§=1, u=0,1, z,=100, z,=110 (damit i= + 1,1), so ergibt sich nach (1) fiir die AuBenverzah- 
“nung v=0,6% und fiir die Innenverzahnung »=0,03 %. 

_ Dié hier vorwiegend behandelten Planetengetriebe sind mit Innen- und AuBenverzahnungen 

in gemischter Bauart ausgefihrt. Bei ihnen kann man nach Budnick gleichen Verlustgrad fiir 
_ jeden Zahneingriff ansetzen, ohne einen grofen Fehler zu begehen. Von dieser Vereinfachung 
_ist daher auch bei den folgenden Wirkungsgrad-Rechnungen Gebrauch gemacht worden. Es 

wurde durchweg mit einem Verlustgrad von 1% je Verzahnung gearbeitet. Legt man aber diesen 
Wert auch fiir reine Innengetriebe zugrunde, so wiirden diese beim Vergleich zu ungimstig 


_beurteilt. Man mu hier geringere Verlustgrade ansetzen, die man tiberschlagig nach (1) ermitteln 
kann. 


3. Die dreiradderigen Planetengetriebe I und II. Da nach einer Beziehung zwischen dem 
Wirkungsgrad des Planetengetriebes und dem des Standgetriebes gefragt ist, wird zu unter- 

- suchen sein, ob und in welcher Weise beide Getriebeformen miteinander zusammenhangen. Das 

-Planetengetriebe I nach Abb. 1 geht in ein Standgetriebe tiber, wenn der Steg festgestellt und 
das Hohlrad beweglich gemacht werden. Kinematisch kann man diese Umwandlung auch so 
deuten, als ob das gesamte Getriebesystem mit einer Drehzahl umlaufen wirde, die gleich gro, — 
aber entgegengesetzt der Stegdrehung ist. Der Steg ist dann relativ zum Gehduse in Ruhe, und 
das derart rotierende Getriebe kann als ein Standgetriebe mit angenommenem bekanntem Ver- 
lustgrad aufgefaBt werden. Die Uberlagerung des Planetengetriebes mit dieser Zusatzdrehzahl 
veradndert weder die Relativdrehzahlen seiner Glieder untereinander noch die Krafte, Momente, 
Leistungen und Verluste. Ventilationsverluste. kénnen wie die Lagerverluste auBer acht bleiben. 
Nimmt man bei diesem Getriebe an, dem Rad | werde die Leistung N, zugeftihrt und der Steg 
fiihre die Leistung Ns ab, so weist das durch die erwahnte Zusatzdrehung entstandene Getriebe 
als Antriebsglied ebenfalls das Rad 1 auf, wahrend der Abtrieb vom Steg auf das Hohlrad 2 
uibergegangen ist. Die Verlustleistung dieses Standgetriebes ist dann 


Nv=N,—N,- (2) 


Der Wirkungsgrad des Zahneingriffs der Rader 1; 3 wird vereinbarungsgemaf gleich dem der 
Verzahnung 3;2 gesetzt. Wird er mit 7 bezeichnet, so folgt mit N,=(N, 7) yo 


Ny=N,(1—7*) | (3) 
oder mit y=1—v7, auch angenahert, wenn man v?=—0 setzt, da y sehr klein ist, 
Ny=2yvN,. (4) 


Die GréBe N, als die Umfangsleistung des Rades | ergibt sich zu 
N, = Pye ys= Py 04s - 


Der doppelte Zeiger bei w soll kennzeichnen, daf es sich um eine relative Winkelgeschwindig- 
keit handelt, namlich um die des umlaufenden Rades | gegeniiber dem, gemaB der oben getroffe- 
nen Annahme, festgestellten, also ruhenden Steg st. Die Umfangsgeschwindigkeit v,s entspricht 
der Geschwindigkeit, mit der die Zahne in der Verzahnung 1;3 (und auch 2;3) aufeinander ab- 
rollen, der sogenannten Walz- oder Zahneingriffsgeschwindigkeit?. 

Bei dem wirklichen Planetengetriebe ist diese Walzgeschwindigkeit natirlich nicht mit der 
Umfangsgeschwindigkeit identisch. Wie aus dem als bekannt vorausgesetaten Geschwindigkeits- 
plan nach Kutzbach? * hervorgeht, setzt sich die Umfangsgeschwindigkeit eines Rades aus der 


1 Siehe FuGBnoten | auf Seite 39 und 40. 
2 H. v. Thiingen, Z. VDI. 83 (1939), Se 730) 


3 Siehe FuBnote 2 auf $. 39. 
4 R. Zajonz, Die zeichnerische Untersuchung von Stirnrad-Planetengetrieben. Diss. Techn, Hochschule 


Dresden 1938, 
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Wialz- und einer sogenannten Zahnmitnahmegeschwindigkeit zusammen (vgl. Abb. 4). Hier 
ergibt sich 
Ty Wy9 ats ry (@y2— stg) aa i Tsp sto, Wi Sct hah ioe 
(Umfangsgeschwindigkeit) (Walzgeschwindigkeit) (Mitnahmegeschwindigkeit) 


OE, 2 cory im Gehauseé 


fi pee 
N 
SI 
OCs 
ES 
Abtries y 4 ( a he Sk Ane by 
NGetriebemittellinie” | Mitnahme ap Walegeschwindigke/t 


GeSCIW. Ty W549 7(ye—Us¢z2) 


Abb. 4. Der Geschwindigkeitsplan nach Kutzbach fir das Getriebe I. 


Sonnenrand (auBenyerzahnt), 2 = Hohlrad (innenverzahnt), 
Planetenrad (auBenverzahnt), st = Steg, Planetenradtrager. 


= 
ll 


-Entsprechend laBt sich die Antriebs- baw. LE in eine Zahnwalzleistung Nw und 


eine Zahnmitnahmeleistung Ny zerlegen: 
N,=NwitNmy, ’ 
worin 
Nwy= P11 (@12— stg) 


ist. Hier bedeutet der Zeiger bei m, da es sich um relative Winkelgeschwindigkeiten des Rades | 
bzw. Steges st gegenitber dem ruhenden Rad 2 handelt. Die Walzgeschwindigkeit ist gleich der 
Relativgeschwindigkeit des Rades | gegeniiber dem jetzt umlaufenden Steg. Wie erwahnt, sind 
aber die Relativgeschwindigkeiten yon Stand- und entsprechendem Planetengetriebe immer 
gleich, so daB hier gesetzt werden kann . 
st = 2 — Wat + (5) 


Fur die Bemessung der Verzahnung ist die Walzgeschwindigkeit bzw. -leistung ma®gebend, 
bei diesem Planetengetriebe also nur ein Bruchteil der Antriebsleistung N,. Daraus kann be- 
reits gefolgert werden, daB seine Verzahnungen geringer beansprucht sind und sein Wirkungs- 
grad hoéher.sein muB als der des Standgetriebes. Die Gleitbeanspruchung der Verzahnung wird 
ee durch den umlaufenden Steg entlastet, da dieser einen Teil der Antriebsleistung durch Zahn- 
kupplung bzw. -mitnahme tbertragt. 


Wird der Wirkungsgrad des Planetengetriebes mit 7 bezeichnet, so lautet die allgemeine | 


Gleichung fiir seinen Reibungsverlust 
Ny=N, (1-7); (6) 
worin 


Ni =P Oy» 


ist. Da beim Ubergang des Planetengetriebes in das Standgetriebe sich die Verluste nicht ver- 


andert haben, kénnen (4) und (6) gleichgesetzt und daraus ae den Wirkungsgrad der Planeten- 


getriebe I oder II folgende Gleichung gewonnen werden: 


at . 
= ib 2r(1 ae (7) 
Bei Vertauschung von An- und Abtrieb, also bei treibendem Steg, liBbt sich entsprechend er- 
rechnen . 
SeevcN. 
ree | (8) 


In (7) und (8) ist 
1=MSo. Rad/NSteg = Ubersetzungsverhaltnis. 


NM danni tei 
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: Es erweist sich hier als zweckmaBig, i immer als das Verhaltnis der héheren zur kleineren Dreh- 
 zahl auszudriicken, und zwar unabhangig von der Richtung des Kraftflusses. Da dieser dann nicht 
_ mehr aus dem Drehzahlverhaltnis hervorgeht, wie es DIN 868 entspricht, wird jeweils besonders 
: angegeben, ob vom Schnellen ins Langsame oder umgekehrt angetrieben wird. 

Bei Gleichlaufigkeit von An- und Abtrieb erhalt i ein positives, bei Gegenlaufigkeit ein 
“negatives Vorzeichen. In (7) und (8) wie in allen folgenden Gleichungen ist i jedoch mit seinem 
_ Absolutwert einzusetzen, da der Drehsinn bereits beriicksichtigt ist. 

-_ Die Antriebsrichtung ist hier praktisch ohne Einflu® auf den Wirkungsgrad, denn (7) und 
_ (8) ergeben etwa die gleichen Werte. Der Wirkungsgrad ist nur vom Verlustgrad des Stand- 
_ getriebes und von der Ubersetzung abhangig, und zwar fallt er bei konstantem » mit wachsen- 
dem i leicht ab, wie die voll ausgezogenen Kurven in Abb. 5 erkennen lassen. 


_ 4. Die vierraderigen Planetengetriebe II bis X. Fiir die gleichlaufigen Getriebe III und IV 
_ gelten (7) und (8) und entsprechend die stark ausgezogenen Kurven der Abb. 5. Der Wirkungs- 
grad der gegenlaufigen Getriebe V und VI ist ebenfalls nach diesen Gleichungen zu berechnen, 
wobei i lediglich negativ einzusetzen ist. Man erhalt 


rte eat ge Soy 


ese 2y(1 | “\ (9) 
mA (1— 2p) 
sel (10) 


Auch hier ist der Wirkungsgrad praktisch unabhangig von der Richtung des Leistungsflusses. 
5 Im Gegensatz zu den Getrieben I bis IV steigt aber 7 jetzt bei konstantem y mit zunehmender 
_ Ubersetzung schwach an, wie die 
_ strichpunktierten Kurven in Abb. 5 v 
zeigen. Man erkennt weiter, daB der 3 
- Wirkungsgrad bei gleichem i und y 
etwas geringer ist als der der Gleich- 
laufgetriebe, wahrend der entspre- 
- chende Standgetriebe- Wirkungsgrad 
zwischen beiden liegt. Fur eine Uber- 
setzung von t=4 und einen Ver- 
lustgrad von v=2%, beispielsweise 


betragt der Wirkungsgrad der Pla- 


- 


- neten-Gegenlaufgetriebe V und VI ) var j 


95%, wahrend das _ gegenlaufige 
Standgetriebe unabhangig von der AV aes 
Ubersetzung 96%, und die gleich- I VA 098. 
laufigen Planetengetriebe I bis IV 1 5 


-sogar 97% aufweisen wirden. Ah 


Die Getriebe IX und X stimmen I! 7 Ae 

in ihrem grundsatzlichen Aufbau 
mit den Getrieben V und VI iiber- / 
ein. Jedoch dreht der Steg im U7 
ersten Fall schneller, im zweiten y a 
Fall langsamer als das Sonnenrad. 
GemaB der Vereinbarung i >1 lau- Abb. 5. Wirkungsgrad-Schaubild der Getriebe I bis VI. 

ten die Gleichungen fiir IX und X — Gleichlaufige Getriebe I, II, III und IV, 
dann bei treibendem Sonnenrad ig eacvermaneo, Cpimobe’ san aiy te 


n=1—29(i+1) (11) 


———— ——} 


und bei treibendem Steg 
ei y 
ag OREN, oe 
in_beiden Fallen ist 
i=Nnsteg/NSo. Rad (hier immer gréBer als 1) 


Die Gleichungen (11) und (12) ergeben im gesamten 1- und y-Bereich sehr geringe Wirkungsgrade, 
die mit wachsender Ubersetzung stark abfallen, vgl. Kurvenscharen in der linken Halfte von 
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Abb. 6. Aus ihr geht weiter hervor, daf der Wirkungsgrad nicht mehr unabhangig von der 


Antriebsrichtung im Getriebe ist, sondern bei treibendem Steg wesentlich héhere Werte aufweist 
als bei treibendem Sonnenrad. Im letzten Fall kénnen die Getriebe fir bestimmte Werte von 4 
und » sogar selbsthemmend werden!. Die Selbsthemmungsgrenze ist in dem Sehaubild durch 


of Verlustgrad V 
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Gegen/lgufige Getriebe yy Gleich/ladufige Getriebe 


Abb. 6. Wirkungsgrad-Schaubild der Getriebe VII bis X. 


Gleichliufige Getriebe VII und VIII. 
—:—:— Gegenlaufige Getriebe IX und X. 


die Kurve 7=0 gekennzeichnet. In der selbsthemmungsfreien Richtung (treibender Steg, also 
Antrieb vom Schnellen ins Langsame) betragt der Wirkungsgrad dann etwa 50%. 

Die Gleichlaufgetriebe VII und VIII entsprechen in ihrer Gleichung fiir den Wirkungsgrad 
grundsatzlich den Getrieben I bis IV, nur da sie als Besonderheit entgegengesetzten Leistungs- 
flu8 aufweisen als diese. Es liegt hier der interessante Fall vor, da trotz treibenden Rades | 
der Kraftflu8 von Rad 2 auf Rad | ausgeiibt wird, worauf bereits Brandenberger? hinwies. Da- 


1 Siehe FuBnote 2 von Seite 39. 
® H, Brandenberger, Masch. Bau u. Betrieb 8 (1929), S. 249. 
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ene die Walzleistung umgekehrt ab wie bisher, und man erhalt bei treibendem Son- 
nenra 


* Sh 2 pt ; 

P ag Ly Pesae (13) 
und bei treibendem Steg 

- er SW 

: | USAR Byki 1) (14) 


Diese Werte sind in der rechten Halfte der Abb. 6 aufgetragen worden. Der Kurvenverlauf ent- 


spricht grundsatzlich dem der Getriebe IX und X, nur ist der Wirkungsgrad unter gleichen Be- 
dingungen etwas hoher als bei diesen. 


_ 5. Zusammengesetzte Planetengetriebe. Die Zahl der méglichen Bauformen ist bei diesen 
Getrieben nahezu unbegrenzt. An dieser Stelle seien einige Getriebe behandelt, die zu zwei 
Hauptgruppen gezahlt werden kénnen. 
_ Die Getriebe der einen Gruppe sind dadurch gekennzeichnet, da® sie einfach durch Hinter- 
einanderschaltung einfacher Planetengetriebe gebildet werden. Durch Zusammenfassung ver- 
schiedener Glieder zu nur einem Glied entstehen dabei oft Bauformen, die diese Hintereinander- 
schaltung nicht mehr ohne weiteres erkennen lassen. 

Das Kennzeichen der zweiten Gruppe besteht in einer inneren Leistungsverzweigung, die 
durch ein Ausgleichgetriebe, das als Teilgetriebe auftritt, hervorgerufen wird. 


a) Hintereinanderschaltung von Getrieben. Der Wirkungsgrad eines durch Hinter- 
einanderschaltung von Teilgetrieben gebildeten Getriebesatzes ist naturgemaB gleich dem Pro- 
dukt der Teilwirkungsgrade. Die Hintereinanderschaltung der Getriebe I und VII wirde bei- 
spielsweise den Gesamtwirkungsgrad 


H=NIN VI 
ergeben. Fiir 71 war (7) und fiir yy (13) gefunden worden. Die Multiplikation beider Gleichungen 
fihrt nach einigen Umrechnungen zu dem Ausdruck 

; t7(1 — 2) + 2y 
a iy(1 — 2v) + vt : 


(15) 
Es bedeuten 


yn = Gesamtwirkungsegrad, 

itr = Teiltibersetzung des Getriebes I, 
ivi = Teiliibersetzung des Getriebes VII, 

i = iptyi = Gesamtibersetzung, 


y = Verlustgrad je Zahneingriff des Standgetriebes. 


Uber die Teilwirkungsgrade wurde ausgesagt, da® 1 mit zunehmender Ubersetzung geringfigig 
abnimmt, wahrend yy sehr stark heruntergeht. Der Gesamtwirkungsgrad wird also um so 
hoher liegen, je gréBer die Teiliibersetzung ir gegeniiber tyrr ist. 
Bei sehr hoher Gesamtiibersetzung wird an Stelle des Getriebes I besser das Getriebe III 
vorgeschaltet, weil mit diesem héhere Teiliibersetzungen bis etwa 30 zu erzielen sind. Damit 
erhalt man dann auch bei héchsten Ubersetzungen noch einigermaBen gute Wirkungsgrade, 
bei t=1000 und y=2,5% z. B. 37%, wahrend das Getriebe VII alléin hierbei nur étwa 7=2% 
erreichen wide. 

Dieser Vergleich benachteiligt allerdings das Getriebe VII insofern, als sein Verlustgrad nach 
(1) weit unter dem hier angenommenen Wert y=2,5% gehalten werden kann. Leider sind, wie 
schon erwahnt, uber den Verlustgrad von reinen Innengetrieben geringen Zahnezahlunterschiedes 
bisher keine Ergebnisse veroffentlicht worden. Man hat daher so gut wie gar keine Anhaltspunkte 
dariiber, mit welchen Wirkungsgraden man bei sehr hohen Ubersetzungen, die diese Innen- 
getriebe eben bei geringen Unterschieden in den Zahnezahlen (wobei die Ubersetzung des be- 
treffenden Standgetriebes nahe | liegen wiirde) zu erreichen gestatten, rechnen kann, Bei 
dem heutigen Stand der Schmier- und Fertigungstechnik darf man annehmen, da beispiels- 
weise die Ubersetzung i=1000 von dem Getriebe VII mit einem Wirkungsgrad von 60 bis 80% 
verwirklicht werden kann, wahrend das zusammengesetzte Getriebe dabei héchstens auf 50 
bis 60°% kommen dirfte. 

Nach diesen Betrachtungen erhebt sich die Frage, ob sich der konstruktive Mehraufwand 
gegeniiber dem Getriebe VII tiberhaupt gelohnt hat. Dessenungeachtet ist aber das obige zu- 
sammengesetzte Getriebe in vereinfachter Form haufig angewendet worden. Die Vereinfachung 
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besteht darin, daB die beiden feststehenden Hohlrader sowie die damit eingreifenden Planeten 
rider zu je einem Rad zusammengefaBt werden. Auf diese Weise entsteht das in Abb. 1 mit X 


if. 
bezeichnete Getriebe, das offenbar erstmals von Wolfrom* angegeben. wurde. Abgesehen vonil) | 


. . al hic . : | 
konstruktiven Unterschieden stimmt dies Getriebe vollig mit dem obigen aus I und VII zu 


sammengesetzten Getriebe tberein. Der Wirkungsgrad beider Getriebe kann aus der Abb. if 
ermittelt werden. In dieser ist die Gesamtiibersetzung i als Abszisse, das Verhaltnis iz/y als: 


Ordinate und 7 als Parameter dargestellt worden, wobei die Gleichung (15) zugrunde zu legen) 


ist. Die 7-konst-Kurven werden dann mit groBer Annaherung Geraden, und fur gleiche Werte 


Teildbersetzung  Verlustgrad 
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Abb. 7. Wirkungsgrad-Schaubild des Getriebes XI. 


iz/v erhalt man mit geniigender Genauigkeit denselben Wirkungsgrad. Fiir eine Gesamtiber-_ 
setzung 1=100, eine Teiltibersetzung 1;=6 und einen Verlustgrad v=2,5°% laBt sich z. B. ein 


Wirkungsgrad von étwa 55% ablesen. Durch Erhéhung der Teiliibersetzung auf beispielsweise ,, 


i;=10, die baulich durchaus tragbar ist, kann eine beachtliche Wirkungsgradsteigerung erzielt 
werden. Auf noch héhere Werte kommt man mit einem dreistufigen Planetenradsatz, d. h., wenn 


man an Stelle des Getriebes I das Getriebe III als Teilgetriebe verwendet. Bei tragbarem AuBen- | 
durchmesser erreicht dies ohne weiteres eine Teiliibersetzung von etwa 25, die zu einer Steigerung | 


des Wirkungsgrades unter gleichen Bedingungen bis auf mindestens 80° fiihren diirfte. | 


b) Planetengetriebe mit innerer Leistungsverzweigung. Diese Getriebe unter- | 
scheiden sich von den bisher behandelten einfachen Bauformen dadurch, da deren festgestelltes 
Glied durch Zwischenschaltung einer Ubersetzung vom Getriebe selbst bewegt wird. Von den | 
zahlreichen Losungsméglichkeiten seien hier die in Abb. | schematisch dargestellten Getriebe XII | 


bis XVII behandelt, die der Praxis entnommen sind. 


Man erkennt, daf} alle sechs Getriebe ein Ausgleichgetriebe mit den Gliedern 1, 2, 3, st sowie 


ein Ubersetzungsgetriebe mit den Radern 1’, 2’, 3” (bzw. I’, 2’, 3’, 4” beim Getriebe XVII) ent- b 
halten. Die Ubersetzung dieses Teilgetriebes ist im folgenden immer mit 1’ bezeichnet worden. 


1 U. Wolfrom, Werkstattstechnik 6 (1912), S. 615. 
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_. Um bei diesen Getrieben die Formel fiir den Wirkungsgrad aufstellen zu kénnen, muB zu- 
-nachst untersucht werden, in welcher Weise sich die Leistung in dem Ausgleichgetriebe verzweigt. 
Da diese Untersuchung aus Raummangel nicht fiir alle sechs Getriebe durchgefiihrt werden 
kann, sei nur die Leistungsverzweigung des Getriebes XII an Hand der Abb. 8 ermittelt und 
-anschlieBend die 7-Formel abgeleitet. 

_ Unter der Annahme, daf® die fest verbundenen Rader 1 und I’ antreiben und der Steg st 
die Leistung abfihrt, mu8 vom Planetenrad 3 im Punkt C Leistung in den Steg abflieBen. All- 
gemein gilt aber fiir ein Leistung abgebendes Glied, da® Kraft und Geschwindigkeit an ihm 
entgegengesetzte Richtung aufweisen, vgl. Vjc und Psc in Abb. 8. Die gleich groben Umfangs- 
krafte P34 und P 3g halten der doppelt so groBen Lagerkraft das Gleichgewicht, wirken also 
entgegengesetzt. 
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Abb. 8. Der Leistungsflu8 im Getriebe XII. 
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-Fahrt man in dieser Weise fort, Kraft- und Geschwindigkeitsrichtung an den Radern zu be- 
stimmen, so kommt man zu folgendem Ergebnis. Dem Planetenrad flie8t nur vom Sonnenrad 1 
Leistung zu. Diese teilt sich in die beiden Leistungen Ns und N, auf. Die Teilleistung N, ist 
von auf en nicht wahrnehmbar. Sie kreist innerhalb des Getriebes und sei mit ,,Blindleistung“ 
bezeichnet, und zwar flieBt sie in ihrem Kreislauf tiber die Rader 2’, 3’ und 1’ wieder in das 
Rad 3 zuriick. Im Rad | vereinigt sie sich mit der Antriebsleistung. Die Leistungsverzweigung 
ist demnach durch folgende Gleichung gekennzeichnet: 


Ni= Na } Ng=Nan N, . (16) 


Zur Ermittlung des Wirkungsgrades setzt man wieder die Zahnreibungsverluste gleich der 
Differenz der Walzleistungen entsprechend dem LeistungsfluB 


Ny=Nan— Nav=(Nuy— Nw) +(Nux — Nur’) - 
Wird je Verzahnung wieder mit konstantem Wirkungsgrad #» gerechnet, so gilt auch 
Nan—Nap=(1—m.2) (Nua + Now’) 
Division durch Nan fihrt auf folgenden Ausdruck fir den Wirkungsgrad: 


Nab ay (Nur Nw 17 
| . Nave e TL aea h a: Oe 
Das erste Glied in der Klammer Ny,/Nan wird aus (16) gewonnen. Darin kann 

Ny = Nyy =v" Nw =e tgt tN eh pS (18) 


gesetzt werden, da das Teilgetriebe 1’ 2’ 3’ ein Standgetriebe ist und die Rader 2 und 2’ fest 
miteinander verbunden sind. Damit geht (16) in 
Ny =Nantno? No (19) 


uber. 
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Werden die Umfangsleistungen N, und N, auf Grund der Beziehungen * 
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durch die Walzleistungen ersetzt, so erhalt man aus (19) 
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Abb. 9, Wirkungsgrad-Schaubild des Getriebes XII. Annahme; Verlustgrad je Verzahnung »=1%. 


Darin bedeuten 


und 


Das zweite Klammerglied der Gleichung (17) Nwy/Nan folgt aus (18) und (21) und ergibt 


sich zu 


Nw ny Gi—1) 
Nan t+i/—nAG—I) ¢ 
Fihrt man (22) und (23) in (17) ein und setzt naherungsweise 


(23), 


ny=(1—»)? ~ 1—2y, 
ni=(1—v)4 w 1-49, 
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by 


80 erhalt man die gesuchte Beziehung fiir den Wirkungsgrad 


+1 2vi/(1— =) 
y eRe (24) 
+14 49(— 1) 
s Diese Gleichung ist in Abb. 9 fir einen Verlustgrad y=1°% zur Darstellung gebracht worden. 
4 Wie man sieht, weisen auch hier die 7-Kurven linearen Verlauf auf, wenn als Abszisse die Teil- 
: tbersetzung VY und als Ordinate die Gesamtibersetzung i gewahlt werden. Ein Vergleich mit 
Abb. 7 zeigt, daB das Getriebe XII nur wesentlich geringere Wirkungsgrade zu verwirklichen 
_ gestattet als das Getriebe XI, wenn gleiche Bedingungen zugrunde gelegt werden. Eine weitere 
_ Uberlegenheit des Getriebes XI gegenttber dem erstgenannten besteht in dem erheblich geringeren 
_ Bauaufwand. Die praktische Anwendung des Getriebes XII, auf das Eichelkraut im Jahre 1913 
das Patent DRP. 266635 erteilt wurde, kann daher nicht empfohlen werden. 
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Abb. 10. Wirkungsgrad-Schaubild des Getriebes XVII. Annahme: Verlustgrad je Verzahnung y=1%. 


Zu der Gruppe der hochiibersetzenden, jedoch sehr unwirtschaftlich arbeitenden Getriebe 
ist auch das Getriebe XVII zu zahlen, das nach Arlt! als Aufzuggetriebe verwendet worden ist. 
Die Ableitung der Beziehung fiir den Wirkungsgrad, die entsprechend den obigen Ausfithrungen 
zu erfolgen hat, wurde an anderer Stelle durchgefihrt’. Von ihrer Wiedergabe sei deshalb hier 
abgesehen und nur kurz Abb. 10 besprochen, in der sie dargestellt ist. 

Man erkennt, daB der Wirkungsgrad dieses Getriebes noch unter dem des Getriebes XII liegt. 
Die Ursache ist auch hier die auBergewohnlich hohe Blindleistung, von der das Getriebe durch-’ 
flossen wird, das somit als sehr ungtinstig zu beurteilen ist. Folgendes Zahlenbeispiel bestatigt 
dies Werturteil. Bei einer Ubersetzung 1=200 und einem Verlustgrad »y=1% erreicht das Ge- 


1 W. Arlt, Zusammenhange und Gesetze der Drehzahlverhiltnisse ebener Raderwerke. Mitt, des Verbd. 
Alter Herren d. Staatl. Akademie fiir Technik, Chemnitz, 1936. 
2 Siehe FuBnote 2 von Seite 39. 


50 - Poppinga: Der Wirkungsgrad der Planetengetriebe. 


triebe XI einen Wirkungsgrad von rund 70 °/,, das Getriebe XII etwa 57 “/ und das Getriebe XVII 


nur 48°. Zusammenfassend kénnen die drei Getriebe dahingehend beurteilt werden, daf sie | 


zwar bei verhaltnismabig kleinem AuBendurchmesser und geringer Anzahl von Radern héchste 
Ubersetzungen ermdglichen, da® ihr Wirkungsgrad dabei jedoch meist untragbar niedrig ist. 


Allenfalls kann er noch beim Getriebe XI in zulassigen Grenzen gehalten werden, wahrend die~ 
beiden anderen Getriebe so unwirtschaftlich, arbeiten, da von ihrer Anwendung abgeraten | 


werden mub. 


giinstige Wirkungsgrade. Unter normalen Betriebsbedin- 
gungen kénnen sie im ganzen Ubersetzungsbereich leicht 


trieben allerdings nicht auf so hohe Ubersetzungen wie 


angeordnet werden. 


Alle vier Getriebe sind mit unerheblichen Abweichun- 


gen zur praktischen Anwendung gekommen. Das Getriebe 
XIII findet man als Rickwartsgang in dem bekannten 
Kraftwagen-Wechselgetriebe von Wilsont. Die Firma 
Demag baute friher in ihre Elektroziige Getriebe ein, 
die grundsatzlich dem Typ XIV entsprechen®. Die 
Albatros-Werke haben in ihrem Boco-Flaschenzug 
das Getriebe XV mit einer Ubersetzung i= 86,5 verwendet. 

Die Untersuchung dieser Getriebe auf Leistungsver- 
zweigung fiihrt zu dem Ergebnis, daB in XIII und XVI 
eine Blindleistung auftritt. Sie ist hier jedoch nur immer 
ein Bruchteil der nach aufen tibertragenen Leistung; 


RTI 
Wi 
Yl 


Slate 
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Die Getriebe XIII bis XVI erreichen durchweg sehr — 


iitber 90°% gehalten werden. Man kommt mit diesen Ge- | 


mit den drei obigen Getrieben. Die obere Grenze durfte | 
etwa bei 1600 liegen, wenn zweistufige Zwischenrader 


wahbrend sie bei den Getrieben XII und XVII ein Viel- ~ 


N,>N, oder N,> Ny 
im Sondertall N, =N, 


Abb. 11. Die drei Méglichkeiten einer Leistungs- 
verzweigung innerhalb von Getrieben, die als Teil- 
getriebe ein Ausgleichgetriebe enthalten. Verluste 
sind nicht beriicksichtigt. 

Oberes Bild; 

Getriebe XII und XVII: Ny=N 4, +2 . Blind- 
leistung N, ein Vielfaches der tibertragenen Lei- 
stung Ny, bzw. N4,- Ungiinstige Verzweigung! 
Mittleres Bild; 

Getriebe XIII und XVI: Ny=N 4,+N,._ Blind- 
leistung ein Bruchteil der iibertragenen Leistung. — 
Brauchbare Verzweigung ! 

Unteres Bild: 

Getriebe XIV und XV: Nj4,=N,+N,. Keine 


Blindleistung vorhanden. Ideale Verzweigung! 


faches derselben betrug. Eine derart geringe Blindleistung 
ist durchaus tragbar, und die Getriebe sind in jedem Fall 
verwendbar. 


Die Getriebe XIV und XV sind blindleistungsfrei. Bei 


ihnen teilt sich die zugefihrte Leistung in zwei Zweige, die. 


. sich im Abtriebsglied wieder vereinigen. Durch die einzel- 


nen Verzahnungen flieBt nur ein Bruchteil der Antriebs- 
leistung.. Die Verzahnungen werden also geringer bean- 
sprucht als bei einem gewohnlichen Getriebe, dessen 
Verzahnungen nacheinander von der Gesamtleistung 
durchflossen werden. 

Nach obigem kann man also bei einer Leistungsver- 
zweigung in einem Planetengetriebe drei Falle unterschei- 
den, die in der Abb. 11 veranschaulicht wurden. Sie gelten 
grundsatzlich fiir alle Planetengetriebe mit mehreren Lei- 
stungszweigen, worauf bereits v. Thiingen® hingewiesen hat: 


Der Wirkungsgrad der obigen vier Getriebe ist praktisch unabhangig von der Teiliiber- 


setzung t’ im Teilgetriebe 1’ 2’ 3’. Auch die Gesamtiibersetzung i hat nur im unteren Bereich 
merkbaren EinfluB. Ab i>10 kann der Wirkungsgrad als konstant gelten, und bei Werten i> 20 
ist er fiir alle Getriebe annahernd gleich groB. Er ist dann nur noch von der GréBe des Verlust- 
grades y abhangig. Aus diesen Griinden sind Schaubilder der bisherigen Art hier nicht angebracht, 


und man kann sich damit begniigen, den Wirkungsgrad in Abhangigkeit von der Ubersetzung — 


bei verschiedenen »-Werten darzustellen. In Abb. 12 sind derartige Kuryen aufgetragen worden, 
wobei der Verlustgrad zu 1% gewahlt wurde. Um einen Vergleich mit den vorher behandelten 
Getrieben zu erméglichen, sind in dieser Abbildung auch fiir sie die entsprechenden Kurven 
wiedergegeben. Die Diskussion dieser Kurven und die Bewertung der einzelnen Getriebe wird 
im folgenden Abschnitt vorgenommen. 


1 Siehe FuBnote 4 von S. 41. 
2 Siehe FuBnote 1 von S. 40 und FuBnote 2 von S. 44. 
3 Siehe FuBnote 2 von S. 41. 
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6. Bewertung aller Getriebe. Den hichsten Wirkungsgrad weisen nach Abb. 12 die Getriebe 
I, H, 11 und IV mit gleichem Drehsinn von An- und Abtrieb auf. Sie tbertreffen die entsprechen- 
den Standgetriebe um ein geringes, fiir die die Gerade 7=0,98 gilt. Die gemeinsame Kurve der 
vier Planetengetriebe liegt etwas oberhalb dieser Geraden und fallt mit steigender Ubersetzung 
sanft ab, wobei die Gerade die waagerechte Asymptote darstellt. 

: Die Wirkungsgrad-Kurve der gegenlaufigen Planetengetriebe V und VI, bei denen der Steg 
dangsamer umlauft als das Sonnenrad, nahert sich der genannten Geraden asymptotisch von 
‘unten. Hier steigt also der Wirkungsgrad mit wachsender Ubersetzung allmahlich an. 


Wirkungsgrad 7 
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Abb. 12. Wirkungsgrad. aller behandelten Getriebe in Abhangigkeit yon der Ubersetzung 
bei einem Verlustgrad je Verzahnung von y=1% = const. 


Gleichlaufige Getriebe. —-—-— Gegenliufige Getriebe. 


Wahlt man bei diesen beiden Getrieben die Zahnezahlen derart, daf der Steg schneller dreht 
als das Sonnenrad, erhalt man sehr geringe Wirkungsgrade, wie die Kurven VIII und IX erkennen 
Raaben. Dabei liegt die Kurve des gegenlaufigen Getriebes IX unter der des Gleichlaufgetriebes VIII. 


Die Hohlradgetriebe VII und X sind in Abb. 12 nicht aufgenommen worden, da ihr port 
grad bei geringem Zahnezahlunterschied der eingreifenden Rader wesentlich kleiner als Le 
gehalten werden kann, vgl. Gleichung (1). Mangels entsprechender Versuchsergebnisse konnen 
dartiber jedoch keine zuverlassigen Zahlenangaben gemacht werden. Nur bei gréBeren Zahne- 
zahlunterschieden nahern sich die 7-Kurven denen der Aubengetriebe VU und IX an. Als 
Anhaltspunkt kann dafiir etwa gelten, dafi die Planetenrad-Zahnezahl mindestens unterhalb 
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der halben Hohlrad-Zahnezahl liegen mufB, so da mehrere Planetenradpaare auf dem Umfang 
angeordnet werden kénnen. 

Im Ubersetzungsbereich bis etwa i=30 mu von der Anwendung der Getriebe VII bis xX 
abgeraten werden. Die Getriebe I bis IV sind wegen geringeren Bauaufwands und hoheren 


Wirkungsgrades dabei viel giinstiger. Dies wird jedoch anders bei hohen Ubersetzungen. Wegen 


zu groBen AuBendurchmessers sind den letzteren hier obere Grenzen bis etwa i=30 gesetzt 
(Getriebe IV), wahrend die Getriebe VII bis X bis zu héchsten i-Werten mit sehr kleinem AuBen- 
durchmesser bauen. Dies gilt besonders fiir die Innengetriebe VII und X, die bis zu Ubersetzungen 


i=30000 ausgefithrt worden sind!. Derart hohe i-Werte sind nur méglich, wenn Planetenrad | 
und Hohlrad um ein bis zwei Zahne voneinander abweichen?, was sich auf die Zahnreibungs- 
verluste sehr giinstig auswirkt und zur Folge hat, da der Wirkungsgrad selbst bei diesen auBer- | 


gewohnlich hohen Ubersetzungen nach (1) in tragbaren Grenzen gehalten werden kann. 


Die AuBengetriebe VIII und IX sind normalerweise nicht fiir so hohe Ubersetzungen ge- 
eignet, jedoch sind durch besondere Verzahnungen auch hier Ubersetzungen bis ungefahr 1=1000_ 
bei tragbarem AuBendurchmesser erreicht worden. Da hierbei aber immer unginstigere Zahn-_ 


eingriffsverhaltnisse vorliegen als bei Innenverzahnungen, ist auch der Wirkungsgrad geringer. 


Er dirfte bei t=1000 z. B. auch bei bester Ausfihrung der Verzahnung und giinstigsten Schmier- 


bedingungen nicht weit iiber 10° gebracht werden kénnen. y 

Das Getriebe XI ist aus den Getriebe I und VII zusammengesetzt und fiir sehr hohe Uber- 
setzungen geeignet. Es liegt mit seinem Wirkungsgrad, wie die Kurve XI deutlich zeigt, wesent- 
lich giinstiger als die Getriebe VII bis X. 


Auch die Getriebe XII und XVII weisen bei héchsten Ubersetzungen nur geringen Durch- | 
messer auf. Die Kurven in Abb. 12 lassen jedoch ihren schlechten Wirkungsgrad erkennen, | 


weshalb ihre praktische Anwendung abzulehnen ist. 


Bemerkenswert ist der hohe Wirkungsgrad der vier letzten Getriebe XIII bis XVI. Wie 


die Kurven zeigen, liegt er durchweg bei allen Getrieben oberhalb 95%. Entsprechend der | 
steigenden Kurvencharakteristik weist allerdings das Getriebe XIII bei kleineren Ubersetzungen | 


einen schlechteren Wirkungsgrad auf, bleibt aber immer noch tiber 90%. Die Kurven der tibrigen 


Getriebe fallen mit steigender Ubersetzung allmahlich ab. Die gemeinsame Asymptote aller | 


vier Getriebe ist die Gerade 7 =96 %. Dies ist der Wirkungsgrad, den ein entsprechendes Stand- 
getriebe mit ebenfalls sechs Radern und vier Verzahnungen aufweisen wiirde, und zwar un- 
abhangig von der Ubersetzung. 


7. Zusammenfassung. Es werden Rechenverfahren und Diagramme entwickelt, die beim 
Entwurf von Zahnrad-Planetengetrieben zur tiberschlagigen Ermittlung ihres voraussichtlichen 
Wirkungsgrades dienen kénnen. Dabei werden insgesamt 17 Planetengetriebe behandelt, die 
im Maschinenbau verbreitet Anwendung gefunden haben. 

Allgemein ist aus den Ergebnissen zu schlieBen, da die Zahnreibungsverluste bei den ein- 


zelnen Getrieben selbst unter Annahme gleicher Bedingungen, wie Bearbeitung der Zahne, | 


Schmierung usw., sehr verschieden gro und insbesondere stark von der Ubersetzung abhangig 
sein kénnen, wie man es in diesem Ausmafs von den normalen nicht umlaufenden Standgetrieben 


nicht gewohnt ist. Diese sind bestimmten Bauformen von Planetengetrieben im Wirkungsgrad 


etwas unterlegen, anderen wiederum weit itiberlegen. 

Hochtbersetzende, raumsparende Planetengetriebe sind gewéhnlich in ihrem Wirkungsgrad 
gefahrdet. Durch Verwendung von Innenverzahnungen geringer Zahnezahldifferenz von Hohl- 
und Planetenrad kann der Wirkungsgrad aber auch bei sehr hohen Ubersetzungen in ertrag- 
lichen Grenzen gehalten werden. 

Bei zusammengesetzten Getrieben hangt der Wirkungsgrad unter Umstanden stark davon 
ab, in welchem Verhaltnis man die Gesamtibersetzung auf die einzelnen Teilgetriebe aufteilt. 

AbschlieBend lassen sich die Ergebnisse dieser Arbeit dahingehend beurteilen, daB sie zwar nicht 
Versuchsergebnisse ersetzen, wohl aber eine Ubersicht dariber geben kénnen, welche Wirkungs- 
grade bei den einzelnen Bauformen von Planetengetrieben im Héchstfall erwartet werden diirfen. 


Gebaut von der Masch.-Fabrik Achar A.G., Ziirich, und der Firma Prometheus, Berlin-Reinickendorf. 
Nur bei Sonderverzahnungen méglich, 


Siemens & Halske nach DRP. 584484. 
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Zur Berechnung der Poinsotbewegung. 
Von R. Grammel. 


1. Einleitung. Das Verhalten des unsymmetrischen kraftefreien Kreisels (d. h. eines beliebig 
gestalteten Kérpers, der in seinem Schwerpunkt gestutzt wird und auber der Schwerkraft und 
dem ihr entgegengesetzt gleichen Stiitzdruck keinen weiteren Kraften ausgesetzt ist) kann 
vollkommen anschaulich mit dem Bild der Poinsothewegung beschrieben werden: Abrollen 
des korperfesten (zum Tragheitsellipsoid des Stiitzpunkts ahnlichen) Poinsotellipsoides bei festem 
Mittelpunkt (Stiitzpunkt) auf der invariablen Ebene, die auf dem raumfesten Vektor des Dreh- 
impulses (Schwungs, Dralls) senkrecht steht!. So klar dieses Bild den raumlichen Ablauf der 
Bewegung wiedergibt, sein zeitlicher Verlauf ist damit keineswegs veranschaulicht, und auch 
die Vorrichtungen, die Poinsot, Sylvester, Boltzmann und andere erdacht haben, um die Dreh- 
-geschwindigkeit richtig darzustellen, sind nicht hinreichend durchsichtig. 


Der formelmaBige Ausdruck der Poinsotbewegung spiegelt diese Schwierigkeit ebenfalls 
wieder. Die bekannte Kirchhoffsche Lésung des Problems vermag zwar die Drehkomponenten 
im kérperfesten Hauptachsenkreuz sowie zwei von den drei Eulerschen Lagewinkeln durch die 
Jacobischen elliptischen Funktionen sn, cn und dn in der Zeit t auszudriicken; der dritte Euler- 
sche Winkel erscheint aber als Integral, dessen Auswertung” bis jetzt nicht befriedigt. Die ele- 
-gante Darstellung der Klein-Caleyschen Bewegungsparameter durch Thetafunktionen hat haupt- 
sichlich zum Ziele, die raumliche Gestalt der Bahnen der Kreiselpunkte zu finden sowie die 
Pol- und Spurbahnen, die der Beriihrpunkt des Poinsotellipsoids auf ihm selbst und auf der 
invariablen Ebene beim Abrollen beschreibt?. 


- Nun ist aber die Poinsotbewegung doch die unmittelbare Verallgemeinerung der sogenannten 
regularen Prazession, die der symmetrische kraftefreie Kreisel als allgemeinste Bewegung 
-vollizieht. Wie diese, so kann auch jene gedeutet werden als Uberlagerung einer Eigendrehung 
und einer Prazessionsdrehung, nur mit dem Unterschied, daB beim unsymmetrischen Kreisel 
beide im allgemeinen (d. h. abgesehen von den permanenten Drehungen um die Hauptachsen) 
nicht mehr gleichformig sind, sondern um bestimmte Mittelwerte mit der gleichen zeitlichen 
Periode schwanken, mit der dies auch der Winkel zwischen der raumfesten Prazessionsachse 
und der kérperfesten Eigendrehachse (langste oder ktrzeste Hauptachse des Poinsotellipsoids) tut. 


Das Ziel der folgenden Uberlegung ist, dies in Formeln zu fassen, d.h. die Mittelwerte der 
drei Bewegungskomponenten anzugeben sowie die Fourierreihen ihrer Schwankungen. 


Wir werden uns dabei einer Naherungsrechnung bedienen, die jedoch fir einen groSen Teil 
‘des Bereiches aller Poinsotbewegungen keinen allzu groBen Fehler aufweist und wenigstens 
fiir den vorliegenden Zweck als ausreichend genau angesehen werden kann. 


2. Die Drehkomponenten und die Eulerschen Winkel. Zur Vorbereitung miissen wir zuerst 
die Grundformeln der Kirchhoffschen Lésung zusammenstellen. Es sei ein kérperfestes, rechts- 
handiges Achsenkreuz (x,y,z) zugrunde gelegt, das mit dem Hauptachsenkreuz des Schwer- 
punkts zusammenfallt. In ihm seien @x, Wy, W: die Komponenten des Drehvektors und A, B, C 
die zugehérigen Drehmassen (Haupttragheitsmomente), fur die die Rangordnung 


A Bes O }-oder A a B= C (1) 


gelte. Dann lauten die Eulerschen Bewegungsgleichungen 


A Le LBL C) os 0, (2) 


1 L. Poinsot, Théorie nouvelle de la rotation des corps, Paris 1834. 
2 Die Formel von E. T. Whittaker, Analytische Dynamik, S. 160, Berlin 1924, zeigt in der Form 
- Zap =. 21 mt 
div gl a, oe 
wo g(t) eine periodische komplexe Funktion ist, daB @ noch nicht die ,,mittlere“ Prizessionsgeschwindigkeit 
sein. kann, Auerdem liefert die Spezialisierung auf den symmetrischen Kreisel ein unrichtiges Ergebnis. 
3 F, Klein und A. Sommerfeld, Uber die Theorie des Kreisels, S, 454, Leipzig 1897/1910. 
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und zyklisch weiter. Ihre Lésungen sind 


Wx=% dn ot, 
Wy=f sn ot, (3), 
O:=y cn ot. 
Um die Konstanten «, 6, y, o sowie den Modul k der elliptischen Funktionen dn, sn und cn / 
zu bestimmen, denken wir uns den Kreisel zur Zeit t=0 mit einem Drehstof angetrieben, | 
dessen Achse in der (x, z)-Ebene liegt und dort die Drehkomponenten @xg und Gz erzeugt, 
wogegen wy,)=0 ist. (DaB diese Inbetriebsetzung des Kreisels die Allgemeinheit seiner Bewegung — 


in keiner Weise einschrankt, ist aus dem Poinsotschen Bild bekannt.) Dann ist 


K=DWxo > VY =Wrz9 >» ; (4) | 


und die Lésungen (3), in (2) eingesetzt, liefern die drei Gleichungen 


—Aawok? =(B—C)Pon, (5) 
BB o = (C—A) @xp Wx » (6) 
—Cowmo = (A—B) Bax (7) 


zur Bestimmung von £, o und k?. Multiplikation und Division von (6) und (7) gibt 
Li(de BYA=C)p be © 


o° BC Wx0> (8) . 
C(A—C) , 
B= Sap Oo. (9) 


Division von (5) und (6) endlich gibt wegen (9) 
G(B2.C) a, 


2 
EE AUB ose (10m 


Zufolge (1) sind o,f und k sicher reell; auf die einfache Diskussion ihrer Vorzeichen brauchen wir 
hier nicht naher einzugehen, da dies fiir das Folgende belanglos ist, ebensowenig auf den be- 
kannten Beweis dafiir, da’ man je nach der Rangordnung (1) die Bezeichnung der x- und 2-_ 
Achse stets so einrichten kann, daB 0 < k?=1 wird. Die untere Grenze k?=0 bedeutet hier — 
den symmetrischen Kreisel (B=C), die obere k?=1 den Sonderfall der trennenden Polbahn, 
der uns hier nicht beschaftigt, da die zugeh6rige sogenannte asymptotische Bewegung in jeder 
Hinsicht bekannt ist. 


Die raumfeste Richtung des Drehimpulsvektors D und sein unveranderlicher Betrag D sind 
durch seine beiden Anfangskomponenten Aq@sy und Caz) festgelegt. Von ihm aus definieren 
wir die drei Eulerschen Winkel 0, y und py folgendermaffen: 6 sei der (als spitz voraussetzbare) 
Winkel zwischen dem Vektor D und der positiven x-Achse, g der Drehwinkel des Kreisels fir 
seine Drehung um die x-Achse, betrachtet von einem Beobachter, der die Drehung der Ebene 
(+x, D) um die Achse D mitmacht; diese letzte Drehung werde durch den Winkel y gemessen. 
Zur Zeit t=0 sei der Winkel p=0 fiir die positive y- und z-Achse; ebenso nehmen wir in diesem 
Augenblick den Drehwinkel y=0 fiir die (+x, D)-Ebene. 


Mit diesen Festsetzungen werden die drei Komponenten des Drehimpulses 


Awx=D coso, 
Boy=D sind sing, (11) 
Cw: =D sind cosy. 


Aus (11) folgt mit. (3), (4) und (9) 


cosd = Aan ot, (12) 
t be —~ SD iot 
oe cnot (13) 
mit der positiven Abkirzung 
9 B(A_—C) 
C(A_B) * (14) 
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Endlich projizieren wir die Drehkomponenten noch orthogonal auf diejenige Achse in der (y,2)- 
Ebene, die auf ihrer Schnittgeraden mit der Ebene (+x, D) liegt und mit D den Winkel + —6 
bildet. Da weder wx noch dd/dt und dqg/dt solche Projektionen liefern, ist 


dy 


oe sind = wy sing + @: cos 
oder wegen (11) ; 
1 Dae D 
—* = sin? @ + rei cos? , (15) 


wenn wir den trivialen Fall sind=0 der permanenten Drehung um die x-Achse ausschlieBen. 


An diese langst bekannten Formeln kniipfen wir an, wenn wir jetzt die ,,Eigendrehung » 
und die ,,Prazessionsdrehung* y — wie wir beide auch bei der Poinsotbewegung nennen mégen — 
anschaulicher als bisher darstellen, insbesondere in ihre Mittelwerte und ihre Schwankungen 
zerlegen wollen, um einen tieferen Einblick in diese Bewegungskomponenten zu bekommen. 


3. Die Eigendrehung der Poinsotbewegung. Wir erinnern uns an den aus der Theorie der 
elliptischen Funktionen bekannten Zusammenhang der Funktionen sn und en mit den Theta- 


funktionen #, und #, 
sn ot YB) 


cnot Vk %2(6t) 


(k’ = 1k? ), | (16) 


wobei fiir die Thetafunktionen die Reihenentwicklungen gelten 


aot 
U (ot) = 2 Ca sin 3K 


Sarr are aot ici, e TOE 
qisin 3 ae + q°/4sin 5 2K Heth | 


(17) 
mot mot mot 
GAct) = 2 (as cos aK ae he cos 3 ratte q’!4cos 5 aK +... ‘) ; 
hierbei ist K’ 
Ree eee Ne ute 


und K sowie K’ sind die vollstandigen elliptischen Normalintegrale erster Gattung fiir die Moduln 


k und k’= V1—R, die man mit dem Wert k (10) aus den Tafeln! entnehmen kann. Nach (13), 
(16) und (17) wird also 


qd sins —q?sin3s+q*sin5s—+:.. aot 
tsp=|/ (s= 2%). (19) 


k’ coss+q?cos3s+q®cosos+... 


Man erhalt offenbar einen sehr viel handlicheren Ausdruck, wenn man im Zahler und Nenner 
die Reihen schon nach dem ersten Glied abbricht und also genahert 


eT 
isp |p es (20) 


setzt. Eine Abschatzung des dabei begangenen Fehlers zeigt die folgende Tabelle, in welcher 
k—sin@ ist. Wir wollen diesen Fehler zunachst einfach einmal in Kauf nehmen, kommen aber 
spater noch auf ihn zuriick. 


6) qQ Fehler < 
Se 2354 7 LO SeulOme 
10° 3,69 - 10-6 (Soo Ker 
20° 6,08 - 10-5 2° 10m 
30° Silke LOR Io tke 
40° LACK orm 
50° 3,02 - 10-8 Lets? 
60° 7,40 - 1073 2e LOn2 
70° 1 20s pyonl Ket: 
80° 4,27 « 107” en Ome 


1 Jahnke-Emde, Funktionentafeln, S. 79 oder genauer S. 86, 4. Aufi. Leipzig 1948. 
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Um @ in seine beiden Teile (Mittelwert und Schwankung) zu zerlegen, differentiieren wir (20) 
nach der Zeit und erhalten 


(14 tg 9) = SE U4 (1+ tats) 


der eee! 
sate dy Va 1+ tgs 
dt 2K b+ tats 

und dies formt man leicht um in | 

dp ano VRI (21), 

at yay K ko 
dabei ist 

I iia) 
ty : an i 22 
ae l+<ecos2s me e aa ( )} 


Weil le|<1 bleibt, so kann man R in die stets konvergente Reihe 


R=1—<« cos 2s+e? cos? 2s—e* cos? 2s+—.-- 


entwickeln und danach gliedweise integrieren. Bentitzt man dabei die goniometrischen Formeln | 


Cos 8st eg = a [eos (2m— yt ("% i cos (2n—3) «+ 
Te = es dos (2ri—5) at: ot Go cos «| ; 
cos?" x= ——— Tae : ~ [0s anxt (4 "\ cos (2n—2)x%+ 


+ en cos (2n—4)x-+---+ (A "ie cos 24% 4. '> (| : 


so erhalt man 


[ Rima ee sin + — ay sin 22" +-ay si 


mit den oe 


-SCQ" aber yay. 
Sera ee. 


Diese Koeffizienten kann man leicht berechnen. Schreibt man die Reihen @, a, und a, aus- 
fiihrlich an, so erkennt man sofort, daB 


ies 1 ee he a 
Seren a Carre neg | 
2 
ay aay) 5 “ (25) { 
\ a \ q 
oe | 


ist. Um die weiteren Koeffizienten zu ermitteln, geht man vom Pascalschen Binomialkoeffizienten- 


Satz 
He a aan i ence ) 


aus, fiihrt rechts y=v’ —1 ein und erhilt 


a) (23), 


(24) | 


| 


eo) Sines Ve to) , Phe es / ‘ rates 
EO etme a 
und daher durch Vergleich mit (24) alsbald die Rekursionsformel 
Os a (22s a a2) (n> 2), (26) 
welche alle weiteren Koeffizienten a3, a4,... vollends zu berechnen, d.h. in ¢ auszudricken er- 


laubt. 


11, Band 1950. Grammel: Zur Berechnung der Poinsotbewegung. 57 


Nunmehr folgt aus (21) durch Integration gemaB (23) und mit dem Wert ay von (25) 


p=ot +5 (0). (27) 
Dabei ist 
@ MO 
a ; We = OK (28) 
der Mittelwert der Eigendrehgeschwindigkeit und 
| nen ty [ORES _ wot 

Ft)=— 77 (a sin K — a, sin 2 + a; sin 3 oe — + +) (29) 


die Schwankung von gy um seinen Mittelwert wet, und zwar ausgedriickt als Fourierreihe in der 
Periode 2K/o der Drehung . Hiermit ist — im Rahmen unserer Naherung — die Eigendrehung 
‘der Poinsotbewegung tibersichtlich zerlegt und dargestellt. 

| Im Falle B=C des kraftefreien symmetrischen Kreisels wird k=0, also K=z/2 und also 
“®e=0 = |[(A—B)/B] @xg, wie von der reguliren Prazession dieses Kreisels wohlbekannt ist; 
-auBerdem wird k’=l=1, also e=0, und damit verschwindet die Schwankung, wie es sein muB, 


identisch, da (wie leicht nachzurechnen) a, wie &, a, wie e2 und allgemein a, wie e” verschwindet. 


4. Die Prizessionsdrehung der Poinsothewegung. Man formt (15) um in 


dy D 1— te 

“dt. IBC |(B C) (B C) cee 
oder mit der Naherung (20) 

dy D k’cos* s — I sin? s 

dt BABS? Oy |(B eC) (BG) tenterTans| 
oder etwas anders geschrieben 
S dy OD beat Akl 

SORE Dee OMENG: [(B+o) | B Os kee R)| 


-und integriert dies gemaB (23) und mit Beachtung von a, (25) zu 
| p=Opt+ pit) - (30) 
Hier bedeutet 


D 
2BC 


Op 


Poben te | 
B+0)4+(B-o)+ =| | 31 
| (B+O4 (B= 0) (31) 
den Mittelwert der Prazessionsgeschwindigkeit und 


me D(B—€) 2kh'l = K ee eitog! vg wot nye AOR 
y(t) = BC ee ae (a, sin “2* —a, sin 2" + a, sin 3 K —+-++:) (32) 


die Schwankung von y um seinen Mittelwert wpt, und zwar ebenfalls ausgedriickt als Fourier- 
reihe in der Periode 2K/o der Eigendrehung gy. Damit ist — wieder im Rahmen unserer Nahe- 
rung — auch die Prazession der Poinsotbewegung in ihre kinematischen Bestandteile zerlegt. 


Man sieht aus (31) und (32) deutlich, inwieweit mit B + C diese Prazession von der regularen 
Prazession wp=D/B des symmetrischen kraftefreien abweicht. Mit B=C geht tatsachlich qp (31) 
in diesen Wert tiber und 7(t) verschwindet identisch (wegen an ~ &"). 


5. Der Prazessionswinkel der Poinsothewegung. SchlieBlich wollen wir die entsprechende 
Entwicklung auch noch fir den Winkel 6 zwischen Drehimpuls- und x-Achse durchfihren. Weil, 
wie aus der Theorie der elliptischen Funktionen bekannt, 


qr Ps(2t) 
dn ot= |/k 9,(61) (33) 
ist, wo fiir die Thetafunktionen #, und @ die Reihen 
t 
9, (at) =1+4 2 (qe0s = +-gteos 272% 4.4% cos 3 + +++), | os 
mot not mot 
8, (01) =1—2(q.cos “2 — gt cos 272" + g? cos 33 — + +++) | 


58 Grammel: Zur Berechnung der Poinsotbewegung. Ingenieur-Archiv 
NO Es uN at 9 1 A Rs ac ae SLES il Eo 


gelten, so diirfen wir mit noch kleinerem Fehler als vorhin nach (12) genahert 


t 
142008 


K 
setzen. Dabei ist 6, der Winkel 6 beim Bewegungsbeginn. 


ae, yor : 
cos 0= Vv COS Oo mit cosd9= —, (35) 


1—2q cos 


Entwickelt man (35) in eine Reihe, wobei man folgerichtig nur bis zur Potenz q? fortschreiten 
kann, so erhalt man genahert 
: cosd=uy+ u(t) - (36) 
Hier bedeutet 


up = (1+ 492) |/ k’ cosd, | (37) 3) 


die mittlere Entfernung des Punktes x=1 der x-Achse von der Ebene durch den Schwerpunkt 
parallel zur invariablen Ebene und 


mot wot 


H(t) = 4g] R cos d9 | (1434?) cos K +4 cos 272" +g? cos 3 K tee] (38) 


die Schwankung dieser Entfernung als Fourierreihe mit der Periode von 9g. 


Auch hier sieht man deutlich die Abweichung von der regularen Prazession. Mit B=C, also 


k’=1, K’=c0, q=0 kommt uy=cosd, und u(t) = 0. 


6. Abschatzung des Giiltigkeitsbereichs. Mit Annaherung an den symmetrischen Kreisel | 


(B=C) oder an die stabile permanente Drehung um die x-Achse (@.=0) geht q und damit 
der Fehler in (20) und allen folgenden Formeln gegen Null. Der Fehler wachst dagegen mehr 
und mehr bei Annadherung an den Sonderfall k=1, der bekannt ist als sogenannter Fall der 
trennenden Polbahn, also wenn zu irgendeiner Zeit die Drehachse genau oder nahezu genau 
mit der Hauptachse der mittleren Drehmasse B zusammenfallt. Das sind aber Poinsotbewegungen 
von zweifelhafter Stabilitat, die man sowieso nicht mehr unter dem Bild einer Prazession fassen 
kann. Die Grenze fiir die Genauigkeit unserer Lésung ist also sachlich begrindet: diese Formeln 
kénnen naturgemaB die Poinsotbewegung nur solange durch eine Prazession mit tiberlagerten 
Schwankungen approximieren, als sie selbst einigermafen das Geprage einer Prazession besitzt; 
je mehr sie sich davon entfernt, um so mehr muf unsere Lésung versagen. 


An Hand der Fehlertabelle in Ziff. 3 laBt sich leicht abschatzen, fiir welche Poinsotbewegungen 
unsere Formeln bei einem Fehler von vorgeschriebenem Héchstbetrag noch Geltung beanspruchen 
diirfen. Fiir 6=0 sind sie exakt; der gréBte bei jenem Fehlerhéchstbetrag noch zulassige Off- 
nungswinkel 0, des Kegels der ,,mittleren“’ Prazession kann sofort ermittelt werden. Es gilt 
namlich zufolge (37) und (35) mit den zu gegebenem Fehlerhéchstbetrag gehérenden Tabellen- 
werten @, und q, 


a Aw 
cos 0, = (1+ 4q?) |/ cos O ey 
; ( a ay : | 2 Oxo” + C Wz9" 


oder gemaB (10) wegen k=sinO@, 
(1+ 4q,) cos O, 


SCC B) 
TAB es 


cos 6, mit m 


(39) 


LaBt man beispielsweise einen Fehler von 2% zu, so kommt (laut Tabelle) mit 0,=60°, 
qi; =7,40- 10-8 
1,03 - 0,707 


cos 6 2S es 
+ V¥+0,75 m 


und dies gibt fiir folgende Verhaltnisse A:B:C die dahinter angeschriebenen Werte von 6,: 
A:B:C=4:3:2 ~6;=38,5°, 
A:B:C=2:3:4 0, = 62:57), 
AB: C= 2:07 LO MOga Oa. 


Unsere Naherung ist also innerhalb einer Fehlergrenze von 2°%% fiir alle Poinsotbewegungen — 


zwischen 6,=0 und 69 =0, brauchbar. Man erkennt tberdies, daB unsere Naherung fir d<C 


u ee daB s mfealtodicx sogenannten epizykloidischen Pornrorhensetneen besser 
ae die periaykloidischen. 


uch den Fehler von (36) kann man cachuragtich abschatzen, indem man priift, ob die aus ere 
bis OP fiir t=0 folgende Ate feet ‘ Fe 


os Uy + u(0)= cos dy “4 


oe ag Lidge hag (1b a4 4) = “e | | | 
qiteise: Man findot an Beispielen, da® ihr Fehler regelmaBig arheblich unter'der E ehlergrenze | 
per orelc in Ziff. 3 bleibt. Der Fehler yon (27) und (30) laBbt sich ebenso an der Gleichheit 


0 HO + (oy eu 


“ 


) COS Oy = Wag 
i 


5 
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Das Spektrum isotroper Turbulenz im statistischen Gleichgewicht". 


Von J. Rotta. 


1. Einleitung. Fir die, Eigenschaften der Turbulenz ist die spektrale Verteilung der Energie 
der Schwankungsgeschwindigkeiten tiber die verschiedenen Frequenzen oder die Korrelations- 


funktionen bestimmend. Um diese kennzeichnenden Funktionen fir isotrope, zeitlich abklingende _ 


Turbulenz theoretisch zu ermitteln, sind schon eine Reihe von Versuchen unternommen worden. 


In einigen neueren Arbeiten wird die Tatsache, daf die kinetische Turbulenzenergie von den | 
groBten Turbulenzelementen in die kleineren Wirbel flieBt und dort schlieBlich durch molekulare — 


Reibung vernichtet wird, als Grundlage fiir eine statistische Theorie der Turbulenz benutzt. 


Dabei eréffnet der von W. Heisenberg! angegebene Ansatz fiir die Erfassung des turbulenten’ 
Energieaustausches neue Moglichkeiten zur Behandlung noch offener Fragen. Auf dieser Grund- — 


lage wird hier das Spektrum der Turbulenz untersucht. Eng damit verbunden ist die Behand- 
lung der Frage nach dem zeitlichen Abklingen der Energie. 


2. Physikalische Grundlagen. Das Spektrum der Turbulenz einer inkompressiblen Flissig- | 
keit wird durch eine zeitlich verdnderliche Funktion F'(k,t) definiert, wobei die Wellenzahl | 


k=27/Wellenlange ist, derart, daB die kinetische Turbulenzenergie je Masseneinheit 


EW) == = f F(kn dk (2.1) | 


betragt. Dabei ist 12 der quadratische Mittelwert des Vektors der turbulenten Bewegung; da- | 


gegen ist der lineare Mittelwert von u gleich Null. 
Der physikalische Sachverhalt ist folgender: Der in dem kleinen Spektralbereich dk einer 


bestimmten Wellenzahl k.enthaltenen Turbulenzenergie F(k,t) dk wird durch turbulenten Aus- 


tausch einerseits von gréBeren Wellen dauernd Energie zugefiihrt, andererseits wird aber auch 


auf die gleiche Weise Energie an die Bereiche kleinerer Wellen angegeben; ein dritter Teil der 


Energie wird durch molekulare Reibung in Wairme umgesetzt. Letzterer Anteil ist nur bei den 
kleinsten Turbulenzelementen wesentlich; bei den gréBeren Wellen bleibt er ohne Bedeutung. 


Aus der Summe dieser drei Anteile ergibt sich die gesamte Energieanderung ee dk des 


betrachteten Spektralbereiches. Driickt man diese Uberlegungen in Formeln aus, so ergibt sich 
eine Bestimmungsgleichung fiir das Spektrum, wenn man geeignete Beziehungen fiir die Ande- 


5 F(k, : A 
rung der Energie oe dk hinzunimmt. 


Die gesamte Dissipation betragt 


© ; | 
are kd Borers d E(k, t).. 4, 
Sian iL Soe dk. (2:2)84 
0 . 
Andererseits gilt aber auch (vgl. z. B. die zitierte Arbeit von Heisenberg) 
' © i! 
S(t)= vf 2F(k,t)k2dk (2.3) 
0 
mit » = —— = kinematische Zahigkeit. 


Q 
Die Gleichungen (2.2) und (2.3) reichen zur Festlegung von 0 F'(k,t)/dt nicht aus, so daR 


noch weitere Beziehungen erforderlich sind. A. N. Kolmogoroff? untersucht eine lokalisotrope — 
Turbulenz, die in einem gewissen Bereich isotrop, homogen und zeitlich konstant ist. In diesem | 


Falle ist die ‘von einem bestimmten Spektralbereich_an die Bereiche kleinerer Wellen abgegebene 


Energie gleich der von gréBeren Wellen zugefiihrten Energie, wenn man von der Aufzehrung 
durch molekulare Reibung absieht [0 F(k,t)/dt=0]. Dieselbe Annahme liegt den Untersuchungen | 


* Auszugsweise vorgetragen auf der Tagung der Ges. f. angew. Math. u. Mech. am 22. 9, 1948 in Gottingen. 
W. Heisenberg, Z. Physik 124 (1948) S. 628. 
G. K. Batchelor, Proc. Cambridge Philos, Soc. 43, 4 (1947) S. 533. 


= 


tw 
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on C. F. v. Weizsdcker! zugrunde. Unter Beachtung der Tatsache, da die Turbulenzenergie 
iauptsachlich in den gr6Bten Elementen (Wellenlingen) enthalten ist, stellt die Annahme 
|F’(k,t)/dt=0 auch bei zeitlich abklingender Turbulenz eine sehr gute Annaherung dar, wenn 
nan speziell den Bereich betrachtet, dessen Wellenlangen klein gegen die groBten vorkommenden 
ind. Fur die gréBeren Wellenlangen gibt diese Annahme die tatsachlich obwaltenden Vorgange 
ffenbar nicht in befriedigender Weise wieder. Der Vergleich der Rechnung von W. Heisenberg 
nit gemessenen Spektren von Windkanalturbulenz zeigt, daB die Ubereinstimmung zwischen 
technung und Messung in-einem wesentlichen Wellenbereich zu wiinschen ubrig laBt. 

Es soll deshalb versucht werden, einen anderen Ansatz fiir 0F'(k,t)/0t zu finden, der in einem 
réBeren Wellenbereich Ubereinstimmung mit Messungen ergibt. Hierzu kann uns folgende 
Uberlegung dienen. Der Abklingvorgang der isotrop vorausgesetzten Turbulenz sei quasistationar, 
o da8B sich nach und nach in allen betrachteten Wellenbereichen statistisches Gleichgewicht 
instellt. Dies ist méglich, wenn die Zeit, in der sich die Intensitaét der jeweils betrachteten 
Wellen merklich andert, kleiner als diejenige Zeit ist, in welcher die gesamte Turbulenz- 
nergie ihre Grose merklich andert. Da die zeitliche Anderung der Intensitét der kleinen 
Wellenlangen schneller erfolgt als die der grofSen, andererseits das Tempo der zeitlichen Ande- 
‘ung der Gesamtenergie im wesentlichen durch die gréBten yorkommenden Wellenlangen be- 
timmt wird, ist zu erwarten, daB sich nach Ablauf einer endlichen Zeit in dem Wellenzahlen- 
sereich k>% ein Gleichgewichtszustand eingestellt hat. Dabei ist ? eine unbekannte GréBe, die 
nit wachsender Abklingzeit abnimmt und schlieBlich gegen Null geht. 

Das Kennzeichen des ausgebildeten statistischen Gleichgewichtes ist nun, daB das Aussehen 
les Spektrums unabhangig vom Zustand der Turbulenz im Augenblick ihrer Entstehung ist. 
formal kénnen wir also fiir das Spektrum im ausgebildeten Gleichgewichtszustand schreiben 


F(k,t)=E(t) L(t) fly Re) « (2.4) 
vobei voraussetzungsgemaB f(y, Re) nur eine Funktion der dimensionslosen Verdnderlichen 
y=kL(t) (2.5) 
st und auBerdem von der Reynoldszahl der Turbulenz 
Rees aU (2.6) 


ibhangt. Die Bedeutung der Konstanten x wird im nachsten Abschnitt erlautert. L(t) ist eine 
cennzeichnende Linge der Turbulenz, die allgemein ebenso wie E(t) eine Funktion der Zeit t ist. 
Jie Lange L in (2.4) und (2.5) wollen wir als eine mittlere Wellenlange definieren: 


fF (k, 
Le) = % f = ak. (2.7) 
0 


Jurch (2.4) kommt eine Ahnlichkeitsbeziehung zum Ausdruck. Und zwar gilt die Ahnlichkeit 
tir eine Flache im k, F, t-Raum. In dieser allgemeineren Betrachtungsweise des Problems glaubt 
ler Verfasser einen Fortschritt gegeniiber bisher bekannt gewordenen Ansatzen zu sehen, bei 
velchen vermittels Ahnlichkeitshypothesen, die sich nur auf eine Kurve im zweidimensionalen 
Koordinatensystem beziehen, die Eigenschaften der Turbulenz studiert wurden. 


3. Gleichungen des statistischen Gleichgewichtszustandes. Wir kénnen nun daran gehen, eine 
F(k,t a PaaS AAS : ; 
reeignete Beziehung fiir die Energieanderung ee dk fir die im statistischen Gleichgewicht 


yefindliche Turbulenz zu entwickeln. Wir gehen dabei von der Beziehung (2.4) aus. 
Die Differentiation von (2.4) nach der Zeit t und der Wellenzahl k ergibt dann beziehungsweise 


, Re dL 
IFG9 _ 2 10 fly, Re) + EW) [fly Rd +E 7] SP + | 


at Bore (3.1) 
PEOLO Sea a 

9 F (kt) af (y, Re) e OF Re) x 3.2 

pe ae ee ae ree 


liminiert man in (3.1) die GréBe f(z, Re) mit (2.4) und d f(z, Re)/j x mit (3.2), so erhalt man 


aF(k,t)  dEjfdt F(k,t)] dL/dt | ofly, Re) dRe 1g 4 
Nee Zi F (ht) + | (ht) +h = ry FEO LOTR ar 88) 


tC. F. v. Weizsdcker, Z. Physik 124 (1948) 5. 614. 


as 
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Da sowohl L als auch E eindeutige Funktionen der Zeit t sind, kann man L(t) auch in Abhangig- 
keit von E(t) ausdriicken. Es ist also 
dL dL dE 
died Bodie, 


und 
dL/dt E dL dE/dt dlog L dE/dt 


(Re Bae Pe CORY S) dlogE E 


Wir setzen der einfacheren Schreibweise wegen — 
dlog L(t) _ 


peat Yen 3.4 
dugeey * Soe | 
und fihren nun noch 
dE/dt  S _ 
econ C(t) (3.5) 
ein, wobei € und C von k unabhangige Funktionen der Zeit sind; so ergibt sich aus (2.6) 
d Re 1 $s 
Bee 09 
und wir erhalten dann als endgiiltige Form fiir 0F'(k,t)/dt 
aF (kt) Af ACD, l urea} | 
a =C}F(hst)—C| F(t) hate (= —¢) Ret. (3.7) 


worin 0F'(k,t)/d Re = E(t) L(t) Of (y, Re)/0 Re ist. Fir die weiteren Rechnungen brauchen wir eine: 


Sik ig, 2 i : ‘ | 
Aussage iiber die GréBe und die zeitliche Abhangigkeit von ¢. Man kénnte zunachst einfach) 
¢€ konstant annehmen. Richtiger ist es jedoch, die Funktion F(k,t) im Bereich der kleinsten k- 
Werte zeitlich konstant anzusetzen, wie sich leicht an Hand der ftir die Dissipation geltenden) 


Beziehung (3.15) nachweisen Ja8t. Fir (AS) F =0 geht dann (3.7) tber in 
> 
. aF(h) IF(k)] , 1p OF) 
lim F(R) —C | F(R) +8 - + Ret | 45 Re he = 0: (3.8) 


Diese partielle Differentialgleichung ]4®t sich durch den Ansatz | 
f(z, Re) =konst y5—1q(Re) , (3.9) 

also F(k,t)=konst ELS k-1q(Re) befriedigen. Wegen der Unbestimmtheit von ¢ wird jeder Wert 
&+1 in (3.9) der Beziehung (3.8) gerecht, so daB die Funktion f(y, Re) noch von einem Para- 
meter & abhangig ist, tiber dessen GréBe die bisherigen Uberlegungen nichts aussagen. Physi-} 
kalisch sinnvolle Losungen erhalt man nur fir €>1. 
Setzt man (3.9) in (3.8) ein und eliminiert die unbekannte Funktion g(Re) mit (3.9), so: 
erhalt man eine Beziehung zwischen der GréBe € und dem Parameter &: 


1+ 4/2 | 
ae (3.10)) 


& == 
wobei 


ye hee Pees ORE): ote Re afly,Re) 
Tir Nas ) Re Thay f(y, Re) 0 Re 


ist. Der Wert 7 hat in jedem Intervall, d.h. bei jeder Reynoldszahl einen anderen Wert, wahrend & 
den gleichen Wert im ganzen Re-Zahlbereich beibehalt. Die spatere Rechnung zeigt, daB 7 nur 
sehr klein bleibt, und bei mittleren Reynoldszahlen seinen GréBtwert erreicht (GréRenordnung: 
nx —0,1); fir Re=c und Re=0 wird 7=0. Im wesentlichen ist also (x 1/E; die GréBe C ist) 
nahezu konstant. Die Beziehung (3.10) liefert mit (3.4) den zeitlichen Zusammenhang zwischen’ 
E und L. ; 
Bei der weiteren Rechnung kann man unmittelbar an die von W. Heisenberg’ hergeleiteten} 
Beziehungen anschlieBen. Betrachtet man nur den Teil des Spektrums zwischen den Wellen-} 
zahlen 0 bis k, so ist die Dissipation fiir diesen Teil analog (2.2) 


il 


k 
a OF (Kt) ay | 
Sk | Veen (3. 11)) 
0 7 


1 Siehe FuBnote 1 von S. 60. 
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Nach Einsetzen von (3.7) und (3.10) ergibt die Integration 


E+ 2(€ + ) 0 Re 


wenn E; die in dem Teil von 0 bis k des Spektrums enthaltene Turbulenzenergie 
: k 


Ex = f F(k,t) dk’ (3. 13) 
0 


k 
1 A Noah / 
Se = C [E— eT Pea oe re Re [ A, ax|, (3. 12) 
(1 he 


ist. 

Die Beziehung (3.12) stellt eine ziemlich allgemeine Bedingung fir das statistische Gleich- 
gewicht der Turbulenz dar. Voraussetzungen sind nur, daB das Spektrum einer von den An- 
fangsbedingungen unabhangigen Formen zustrebt und daB F(k,t) im Bereich k—>0 zeitlich 
konstant bleibt, deren Berechtigung dadurch offenbar wird, daB der Betrag S;, der Dissipation 
nach (3.15) fiir k—+0 von héherer Ordnung gegen Null geht als die Energie Ex nach (3. 13). Be- 
kannt ist nun der Verlauf des Spektrums in der Nahe k=0 


F(k) ~ ks-?, (3. 14) 
wobei € ein zunachst frei wahlbarer Parameter ist. Wir wissen lediglich, da8 €>1 sein muf, 
um brauchbare Lésungen zu bekommen, d.h. solche, bei denen fiir k=0 auch F(o,t)=0 ist. 
Im Abschnitt 6 werden wir noch eine weitere Aussage tiber das Verhalten von F(k,t) bei kleinen 
k-Werten und tiber die GroBe & bekommen. : 

Die Auswertung von (3.12) erfordert eine zweite Gleichung ftir Sz, die die turbulenten Aus- 
tauschvorgange in geeigneter Weise erfaBt. Diese steht uns in den von W. Heisenberg aus Ahn- 


lichkeitsbetrachtungen hergeleiteten Beziehungen zur Verfigung. Es gilt dann 
y 2 


Seam f 2 (i) KeaK, (3.15) 
0 


wenn ?% die Summe der Zahigkeiten aus molekularer Reibung und turbulentem Energieaustausch 
darstellt: a. Hit 
Ws sy 
ee atx [YS ) ak’. (3. 16) 
k 


Die GréBe x soll hierbei eine reine Zahl sein, die weder von der Wellenzahl k noch von irgend- 
welchen Stoffeigenschaften abhangt. W. Heisenberg hat versucht, die GréBe x aus den hydro- 
dynamischen Gleichungen rechnerisch zu bestimmen, ist aber zu keinem endgiiltig befriedigenden 
Ergebnis gekommen; lediglich kennt man die GréBenordnung von x. Fir unsere Rechnungen 
betrachten wir x als eine GréBe, tiber die man zunachst noch frei verfiigen kann, und die dann 
spater so festgelegt werden soll, daB sich eine méglichst gute Ubereinstimmung zwischen Rech- 
nung und Messung ergibt. Selbstverstandlich kann dieses Vorgehen nur dann einen physikalischen 
Sinn haben, wenn beim Vergleich mit Versuchsergebnissen verschiedener Herkunft den Berech- 
nungen stets der gleiche Wert von x zugrunde gelegt wird. Notigenfalls wird man % eine gewisse 
Re-Zahlabhangigkeit zubilligen kénnen. Aus Grinden der rechnerischen Behandlung wird die 
Konstante x in die Reynoldszahl nach (2.6) hineingenommen. 

Fir die Rechnung ist nun die Tatsache wichtig, da} der EinfluB der Re-Zahl auf F'(k,t) 
praktisch nur dort zu Auswirkung kommt, wo die kinematische Zahigkeit wirksam wird, das ist 
also im Gebiet der kleinsten Wellenlangen. Im ubrigen k-Bereich ist OF (k,t)/dRe nur klein. 
Fur kleine und maBig groBe k-Werte wird der dritte Summand in der eckigen Klammern von 
(3.12) daher nur ‘klein gegen die anderen Summanden sein. Aber auch dort, wo die Wirkung 
der kinematischen Ziahigkeit wesentlich ist, wird wegen des hier herrschenden raschen Abfalles 
von F(k,t) mit k der von der Re-Zahl abhangige Anteil nur ein unbedeutendes Korrekturglied 
neben dem Betrag Ex, der hier Ex ~ FE wird, bleiben. Man darf deshalb eine die wirklichen Ver- 
haltnisse sehr gut wiedergebende Naherung erwarten, wenn man in (3.12) dF'(k,t)/d Re =0 setzt 
und mit der vereinfachten Gleichung! 


se renC Ea kF(k.1)| (3.17) 


arbeitet. 

1 Pir den Spezialfall €=2 hat W. Heisenberg, Proc. Roy. Soc. A 195 (1948) S. 402, die gleiche Beziehung 

Pefanden und aerate das Spektrum der Turbulenz diskutiert. Die auf die Beziehung (3.17) fiihrenden Uber- 

esrccs hat der Verfasser in von obigen Darlegungen etwas abweichender Form durchgefiihrt, ohne von 
=) . * 

der Existenz dieser Heisenbergschen Arbeit zu wissen. 


! 


j 
; 


Fir hinreichend groBe Wellenzahlen k nahert sich E, einem konstanten Endwert,. maralie ] 
dem Wert E der Gocamienonere nach (2.1) und kF'(k,t) geht gegen Null. Damit geht (3. | 


tiber in 


lim S, =konst . (3.18 

k—+o 
Fir groBe k-Werte ist (3.17) also identisch mit der Bedingung, die den Untersuchungen vo 
v. Weizsticker und Heisenberg! zugrunde liegt. Fiir beliebig groBe Reynoldszahlen wird dariber 


i 
hinaus der Ausdruck | 
| 
| 


a ping ee BEM) gy 


beliebig klein gegen Ex, so daB also bei groBen Re-Zahlen (3.18) auch die asymptotische For 


der vollstandigen Gleichung (3.12) ist. 


4. Durchfiihrung der Rechnung. In (3.17) werden die Gleichungen (3.13), (3.15) und (3. 16)! 
eingesetzt. In Anbetracht der zu erfassenden Bereiche von F'(k,t) und k werden wie bei Heisen- 
berg! logarithmische Koordinaten eingefihrt: 

F (k, t) 
F(kyt) ° | 
wobei k, irgendeine Wellenzahl innerhalb des vorkommenden Wellenbereiches und F (Iigst) den) 
zu ky gehérigen Wert von F(k,t) bezeichnet. Die so umgeformte Gleichung zwischen x und w) 
hat die Form 


logs p= — w = loge = 40 (4)s 


x 


ale ae dx) [aac fae =e (4.1) 


—— 


Hierbei treten drei Parameter auf: 


1. die GréBe &, 


2. die Reynolds-Zahl Fk 

RK — uES (4.2) 
3. das Dissipationsverhaltnis =C 

c (4.3) 


a x |/ ed 

Die Zahl der Parameter 1aéBt sich durch Einfihren einer neuen Variablen fir w jedoch auf zwei 
reduzieren. Wiirde man z. B. e~“ =(1/c?)e~” setzen, so warden nur noch die beiden Parameter 
(ch) und € tbrighbleiben. Als Lésung von (4.1) ist also, wie anfanglich vorausgesetzt, eine zwei- 
parametrige Kurvenschar zu erwarten. Es erscheint jedoch zweckmabig, die Form (4.1) beizu- 
behalten. 

Fur (4.1) lassen sich zwei spezielle Lésungen angeben, die spiter von Nutzen sein werden. | 
Setzt man die rechte Seite von (4.1) gleich einer Konstanten, so hat man den yon Hetserbes 
behandelten Fall, fiir den man die eae Lésung 


cals ss + loge (Aet*+1)+4 (4.4)) 


findet. Darin sind A und Ao Se ae Wir setzen A=0, dann wird 


x 


. x— Ww —4 2 i | 
lim ee a es : — ®& pes SaaS as 3 ) 
ae day eye FAR 4A 
foe} 
li - “+ etx 
aie Je de ae 
cI— >On = 


x 
Setzt man dies in (4.1) ein, so kommt 


1] etx 3 e-4x 
er Tq ae ) ( 4 Ae 440 ) = konst. 


1 Siehe FuBnote 1 von S. 60. 
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Daraus folgt 


: Ape 

pe leichung (4.4) hat somit die Asymptote 

a F (ko) 

3 litn: to Tx 40,151 feat Fe) 5 
4 eacogs ix + 1,151—4 log, ( E ke ) , (4.5) 


und fiir das Spektrum, das.von F(k) ~ k-5/3 mit wachsendem k in I'(k) ~ k~? tbergeht, gilt die 


exakte Form 
a ky 5/3 3 47 —4/3 
F (hk) = F (hy) (5°) [+ ere) (=) | (4.6) 
k 


v 


0 
Eine andere Lésung fiir (4.1) findet man, wenn man > 0 gehen laBt. Das Integral 


wg sy 
fe 2 dx 


x 


‘kann dann neben 1/8 vernachlassigt werden. Die Lésung fiir diesen Fall, der das letzte Stadium 
der Turbulenz kennzeichnet, bei dem die Tragheitskrafte keine Bedeutung mehr haben, lautet 
a Go 2% 

w=—(E-I) xt sig Pt to. (4.7) 


Dabei ist m wieder eine Integrationskonstante. 


Fur die vollstandige Gleichung (4.1) laBt sich keine Lésung in geschlossener Form finden. 
Eine schrittweise Integration kann wegen des breiten (x,w)-Bereiches nur mit sehr weit getrie- 
bener Rechengenauigkeit und mehrfachem Probieren zu zuverlassigen Ergebnissen fihren. 
Daher wird (4.1) durch Reihenentwicklungen gelést!. Um bei der Berechnung der in diesen 
Reihen auftretenden Koeffizienten die in (4.1) notwendigen Quadraturen zu vermeiden, wird 
(4.1) dreimal nach x differentiiert, wobei zwischen den einzelnen Differentiationen geeignete 
Umformungen vorgenommen werden, so daf eine gewéhnliche Differentialgleichung dritter 
.Ordnung entsteht (Striche bedeuten Ableitung nach x): . 


3 x—w 


e 2 [35—26w’+ 3w’?+ 4w”’] = ¢ [42—20w’+ 20? 40% — : (4262 w’ 4 (4.8) 


99 W242 we 3.0 22 10 w’ w+ ww’ + 14w”’ — 20’ w’”)| ‘ 


Ausgehend von der Lésung (4.4) wird w in eine Reihe entwickelt, die fur positive x-Werte kon- 
vergiert (Tafel 1). Ferner wird von (4.7) ausgehend ein Ausatz entwickelt, der fiir negative 
-x-Werte konvergiert (Tafel 2), wobei der Fall =2 eine Sonderstellung einnimmt. 

Der gesamte Kurvenzug wird bei x=0 aus den heiden Asten nach Tafel 1 und 2 zusammen- 
-gesetzt, wobei die beiden Kurventeile stetig und ohne Krimmungssprung ineinander ubergehen 
miissen. ZweckmaBig wird A=w gesetzt. Aus den Bedingungen 


w(0) =w(0)=0, 
w’(o)=w'(o), (4.9) 
w’"(0) =w” (0) 
bestimmen sich dann fiir gegebene Werte & und A (Reynoldszahl) die unbekannten GroBen 
c, A, Bund D; dabei werden w(o), w’(o) und w%(o0) nach Tafel 1, w(0), w’(0) und w’(o) nach Tafel 2 
ermittelt. 

Abb. 1 zeigt einige Naherungen fir £=2, A=0. Die Rechnungen konvergieren sehr gut, 
so daB nur wenige Glieder jeder Reihe benétigt werden. Schon die erste Naherung gibt den 
Verlauf im wesentlichen richtig wieder. Das Spektrum nach der von Heisenberg® angegebeneu 
Naherungsformel wurde ebenfalls mit eingezeichnet. 


1 Herrn Prof. W. Heisenberg ist der Verfasser fiir wichtige Hinweise bei der Lésung zu groBem Dank 
verpflichtet. 
2 Vel. FuBnote 1 von S. 63. 
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Die Gesamtenergie betragt 0 
E =F (kt) hy fice? dx (4.10 
und das Produkt Bs 
EL =F (kyt) f owas. (4.11) 
—o 
~3 -2 H] 0 7 2 
z 
Via 
We ‘ 
x \ —-|--- 
va \ 
; ie \ 
». me 
-— Br) ; 7 1 . “i abe a 
W=- LEAT —— \wetxra 
sg RS | 
+ }— AAAS NS \ 
A Ne" \ 


Via Sean 
3 avira My [\ 
. ‘ fe Nive 
; / r 3 ae . 3 y 
ef XK 


w—aAa 
Abb. 1. Naherungslésungen von (4,1) fir £=2, A=0. 
1) 1, Naherung: 


x0: w=—x +(1,382—1,137 x) eX 1A; 
2 x/3 
— (2,286x+1,278 e 218s as 


5 
x=0: w= «+ 1,340 e72 #3 


4 = — 1,382. 
2) 2. Naherung; 
«<0: w=—x + (0,775—1,710 x) e°* + (0,523—0,426 x) e8* +A; 


+ 0,152-e 


5 ee * 2 x/3 
eres: £20: wos x $0,891 e213 4 9,397 e-4 3 4. 0,069 em (22286x + 1,928 &® */%) 


A= — 1,298. 
3) 3, Naherung: 
x0: w=—x + (0,738—1,991 x) e°* +(0,585—0,499 x) e** +4; 


+43 5 


= we ae Sr ani. 
Baz Os ray 0,823 62 “/3 + 0,339 e74*/3 + 0,243 en 2* 
_— = 2 x/3 
0,064 e— (2,286x + 2,080 e x/ Dae. 


1 


A= — 1,323. 
4) Naherung nach Heisenberg (siehe FuBnote 1 auf Seite 63) 


fir —wo<x*< 0: w=—-x +2 log, (Uhet4!?) eas 
( A = — 1,386. 


Die Quadraturen lassen sich im allgemeinen nicht geschlossen durchfiihren. Nur bei der fir 


a | Ri— 0 angegebenen Loésung (4.7) ist dies moglich. Es ist 
seo ee If! (3) 
lim exw dx eal ? oe 
R>0 | 9 (575) | | 
a ; (4. 12) 
; ae 
lim ede = —— : 
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. Tafel 1. 
- Reihenansatz fiir ee Ot 
w(x) fae eee l oe (44ety-i9 Bae aN 
= 6) is e x e x\—1L 
v=) v i 2, e( sa ai 
+> Be en? [mxtne®! * (144 of 1/3) Sean (1 ¥ 4 x) —pi/3 oie Aet* % 
e *) A . 
ron] a (l+Ae a bral 1 Tia) f 
AaB = Integrationskonstanten. 
BF 
ie Se HO 42/3 Sep 
‘3 R el" A ee dx. 
as 
32 Brats, 
“ig es ae =) ee 9 = 77 (3+ +) c* et; 
32 1485 +41550/é PRC Riere: 
30 = Fa] 1701 ie ms Soe 
ie 32 1559,057 + 3997,143/& + 2867 848/E2 2\2 
“ 4 Se : , is 
a “40 = 44] 1701 = + Pip 
me 16 22 2 16 
Be naa (Ft) Da (G+) ee. 
ss 9 eh Veen? 
a Urea ¢ ae ij enaeg Vaal 
byo9 = 1 boog = */43 
i ; 11046812996 £4812 ay, 
az 110 = 1764 ae 


he = — Se EE mies NED se 


eh 
Baro (97,7664 — 10,6531 £ + 0,0459 £2) £2 


b bito 9+ 64/é mele 
3 A NED) 24 : 


. 
3 


i) 
Or 


_wobei I ) wie tblich die Gammafunktion bedeutet. Fiir die gesuchte dimensionslose Funktion 


f(y, Re) wird dann Hee F(kst) —" 
x: EL co é (4. 13) 
ff e—v dx 
feuds 
kL = e® —- 25 (4. 14) 
pevds 


Re = § ——___- . (4.15) 


Fir das Spektrum kénnen wir im Falle Re=0, ausgehend von (4.7) und unter Benutzung von 
(4.12), eine geschlossene Formel angeben: 


ee eG) 0. 10) 
2 
(F) 


P(E/2) 


mit 


a= 


Tafel 2. 


Reihenansatz fir x«S0: 
L)) Se 2s & % 
Dyes Eee Deen D+) © SY" diy e2#*; 
=n =e 


@, D = Integrationskonstanten. 


| (2-1) e 
A u € 1 fs — eu e —w/2 
ae D = nie =e) lim (f e 2 HER rat yaa ae, )| ° 


E>— 0 


rng Q(QUSE\E ) ay l? 
ae (§—2)(2+6)? ce ‘ 


825-38) (See) cee 


day PVE = 2) (2EREES en 
, 24358 2441008 £+ 1546 £4 86464 eo Fo 
Sue dog = (144) (€—2)7(248" 99 ae 
ee £(2442E—3 &) ecit 
} (2+) (4+) (6+6)? 
Ox 2) &=2: 


= kai—1 ” 
w(x) = w—x + e?*(D+f,x) + iy earn oe pas 
i=2 K=0 


w, D = Integrationskonstanten. 


D Sle + lim ul Saban det ee-o)] : 


&>— 0 
o fi =- i ; 
f Osean pt Dene an 
Seo 576 c2 48 ¢ x 
nl iawen 
* Ja =~ ag 
% SNA ed Ot ROT ey 0 7 2 3 4 Gg 


5. Ergebnisse der Rechnung. 
Fir eine Darstellung der Er- 
gebnisse wird die GroBe 
S(y, Re) =F (k,t)/EL tber y=kL 
im logarithmischen Mafstab auf- — 
getragen. Abb. 2 gibt zunachst | 
einige Spektren fiir unendlich | 
groBe Reynoldszahl, aber mit ver- 
schiedenem Parameter & wieder. 
Der Einflu8 der verschiedenen £ — 
wirkt sich hauptsachlich im Ge- | 
biet kleiner Wellenzahlen aus, | 


> = 43 wahrend fiir groBe k-Werte alle 
i riko Kurven in das Gesetz F(k,t)~ k-/3 | 
ria ual S 4-8 tibergehen. Das asymptotische- 
4 4 7-1 = : Verhalten bei kleinen k-Werten 
S F(k) ~ k5-! ist dabei nicht davon 
es 8 ie | abhangig, ob der benutzte Aus- 
! 


druck (3.16) fir die Turbulenz- 
zahigkeit in dieser Form hier 
noch berechtigt ist. 

Der Einflu8 der Reynoldszahl Re wird in Abb. 3 fiir den Fall E=2 (Abklinggesetz Ei ~ 1) 5! 
veranschaulicht. Bei endlicher Reynoldszahl haben alle Spektren fir groBe k-Werte die von ¢ 


Abb, 2. Gleichgewichtsspektren fiir Re= oo und verschiedene Parameter &. 


II. Band 1950. Rotta: Das Spektrum isotroper Turbulenz im statistischen Gleichgewicht. 69 
CW. Heisenberg? gefundene Asymptote F(k,t) ~k-7. Der Darstellung Abb. 3 ist zu entnehmen, 
daB die Funktion f(y, Re) bei groBen Re-Zahlen im Bereich kleinerer yx-Werte fast unabhangig 


yon der Re-Zahl ist. Dieses Verhalten trifft fur alle anderen &-Werte ebenfalls zu. Damit finden wir 
die in Abschnitt 3 getroffene Voraussetzung 


-dafir, daB die Ableitung 9 f(z, Re)/9 Re in 
(3.12) auBer acht gelassen werden darf, 
weitgehend erfillt. Nur im Gebiet groBer 
k-Werte, wo etwa das* Spektrum in 
das Gesetz F(k,t)~k~’ abbiegt, bleibt 
Of (%; Re)/d Re nicht mehr klein. Fiir groBe 
Reynoldszahlen leitet man z. B. aus (4.5) ab 
ee ttt: 0 F'(k, t) 

; lim cane grea Re 


:: Re>o, k->0o 


= 4F(k,t). 


Zum Vergleich sei der Wert fir k—>+0 
‘bei einer mittleren Re-Zahl (Re ~ 100) 
und mit =2 gegenibergestellt: 


eS ~ —0,1 F(h) . 


Tatsachlich wird das Spektrum fir 
-endliche Reynoldszahlen in diesem Gebiet 
durch die Beziehung (3.17) nicht ganz 
richtig’ wiedergegeben. Die Abweichungen 
vom wirklichen Verlauf sind jedenfalls © 
nicht erheblich; fiir die Grenzfalle Re= co 
‘und Re=0 verschwinden sie. Im ubrigen 
hat das Nichterfilltsein der genauen Be- 
ziehung (3.12) in diesem Gebiet keinen 
-EinfluB auf das tibrige Spektrum. Man 
_k6énnte nattirlich daran denken, diese als 
erste Naherung aufzufassende Lésung von 
(3.17) durch Einsetzen in (3.12) iterativ 
zu verbessern. Dies wurde bisher jedoch 
nicht fiir notwendig erachtet und deshalb 
nicht versucht. 


Abb. 3. Gleichgewichtsspektren fir §=2 (E ~ t“) 
und verschiedene Reynoldszahlen. 


Sehr wichtig ist die Frage nach dem zeitlichen Abklingen der Gesamtenergie, die bei be- 
-kanntem Spektrum beantwortet werden kann. Aus (2.3) ergibt sich 
ive) 
E 
S=sy f 2f(gRe) x? dz. (5.1) 


; ” 0 


In Abb. 4 ist das Integral Ve 2 f(y, Re) v7? dy, das aus den ermittelten Gleichgewichtsverteilungen 
0 


fiir €=2 berechnet wurde, tiber der Reynoldszahl logarithmisch dargestellt. Bei groBen Reynolds- 
: zahlen ist dieser Integralausdruck proportional zu Re, so daB hier der bekannte Zusammenhang 
dE E32 
espana» 


(5.2) 


lim 

Reo 
-besteht. Bei kleineren Reynoldszahlen ergeben sich Abweichungen von diesem Zusammenhang 
und fir den Grenzfall Re—0 wird 
lim dE E 
Reo dt aE * 


(5.3) 


Fiir andere Parameter € ist der Verlauf des Integrals ttber Re ganz ahnlich wie in Abb. 4. 


_,  1.FuBnote 1 von S. 60. 
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Wenn man das Dissipationsgesetz ermitteln will, braucht man zundachst den Zusemmenhons 
zwischen E, L und Re. Die Differentiation der in logarithmischer Form geschriebenen Gleichung 


(2.6) ergibt 


d log Re ai L dlogL 1 an : (5.4) 
dlog Hise 22 dlog E 


| 
Setzt man fur ¢ den in (3.10) ermittelten Wert ein, so kann man (5.4) geschlossen integrieren | 
Die Werte E und L lassen sich dann fir jede Reynoldszahl berechnen, wenn sie ftir irgendeine 
Re-Zahl Re= Re, bekannt sind: | 
QE 2 ;| 


2 1 
Rey \==3 E f(y, Re) le : 
Re a ee re 


Mit den so bekannten Zusammenhangen ist die Integration von (5.1) méglich, die den zeitlichen: 
Ablauf der Dissipation ergibt. Fir die Grenzfalle Re—->0 und Re—>0 besteht nach (3.4) de 
einfache Zusammenhang zwischen L und E 


1 
LE (3-9) 


(5.5) 
L(Re) = L(Re,) ( | 


| 
da hier die Funktion f(y) fir y—>0 nicht mehr von Re abhangt. Das Abklinggesetz kann durch: 
Einsetzen von (5.6) in (5.2) und (5.3) ermittelt und in geschlossener Form angegeben werden:: 


lim Re—-oo | lim Re>0 | 
5 | 
SeRUERE SATE! aa mae eee | 
‘ ‘ a . E eats At / 
i S “: a a 
go] -— Sth PRee, iP eeee eeone 
50 eS a 
: Z | 2—§ a 
40 5+ Soe 
ee: us Re| ~1t2+é ib aS 
30 a 
r | a 
th fn 
A 30 Asympt Re-~ co 
5% i a S~ a Der gesamte Abklingverlauf der 
ie “ | “ T es Energie von grofen Re-Zahlen be- 
10 We ginnend bis Re—0 kann nur durch 
a [see ce ew ewer numerische Integration gewonnen 
Hs CESS werden. Das Abklinggesetz wird 
4 . s | 
05 [ae 7 durch &€ wesentlich beeinfluBt. | 
ae loge Re Fiir den Fall £=2 gibt (3.17) die: 
7 0 7 2 3 m7 5 6 7 exakte Lésung von (3.12) an, wie man | 
leicht sieht. Nach (3.10) wird exakt | 
Abb. 4. Dissipation abhangig von der Reynoldszahl fiir £=2. C= 1/2 Flom (3.6) folgt, dal dia 


Reynoldszahl wahrend des ganzen Ab- | 
klingvorganges konstant bleibt. Es gilt das Dissipationsgesetz E ~ t~! und die mittlere Wellen- | 
lange wachst wie L~Yt. Das Spektrum bleibt fiir ¢=2 beim Abklingen vollkommen 4hnlich. 
Doch scheint kein hinreichender Grund vorzuliegen, dem Fall £=2 wegen dieser Eigenschaft 
eine bevorzugte Stellung einzuraumen. Fiir das Verhalten des Spektrums bei den gréBten Wellen- | 
langen, das ja durch die GroBe € beschrieben wird, sind namlich andere physikalische Vorgange 
verantwortlich, als die Einwirkung der kinematischen Zahigkeit. 


An Hand von (3.7) kann man den Abklingvorgang der Funktion F(k,t) in den einzelnen| 
Wellenbereichen quantitativ verfolgen. Mit C nach (3.5) laBt sich (3.7) fur eine bestimmte | 
Wellenzahl k=k tber ein begrenztes Zeitintervall integrieren, wenn man ¢~ 1/€ setzt: 

= = = 

; 1 k dF (k, t) 1 1 Re oF (k,t) 
F(k,t) ~ E% mit q=1 1 =: = - = : Har’. 
eh) a = ( | F(k,t) dk EG: =) F(k,t) Re (5:7) 


F'(k,t) nimmt also mit einer von 1 verschiedenen Potenz von E ab. Der Exponent hangt dabei 
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ae ko F(t 
hauptsachlich von —— (b, Us ab. In den verschiedenen k-Bereichen gilt fiir groBe Reynolds- 
Be BRO. ak 
zahlen: i 
k liegt im Gebiet: F (k,t) 
lim F(k) ~ ké-} wx konst 
k—+>0 } 
max F'(k,t) je. o 0| RAS s 
F(k,t) ~ k-9/8 ~ Ei+2/36 
F(k,t) = k-? ~ Bs+rs 


Man ersieht hieraus, da die kinetische Energie der kurzen Wellen schneller abstirbt als die 
der langen. Die spektrale Funktion F'(k,t) der Wellenzahl k, welche anfanglich im Gebiet F'(k) ~ k5—} 
liegen modge, andert ihre Grofe zunachst nur wenig. Nach einer gewissen Zeit, wenn die Energie 
des Bereiches k>k entsprechend abgeklungen ist, wird bei k etwa das Maximum der Spektrums- 
kurve liegen und die Energie jetzt etwa wie F(k, t) ~ E1—"'5 kleiner werden. Nach Ablauf einer 
weiteren Abklingzeit, in welcher wieder hauptsachlich die Energie des Bereiches k>k abnimmt, 
wird dann bei k das Gesetz F(k,t) ~ k—5/3 gelten und F(k,t) nun wie E!+2/3€ abnehmen. Nach 
weiterem Abklingen der in k>k enthaltenen Energie wird dann das Spektrum in der Nahe 
von k schlieBlich dem Gesetz F(k,t) — k-7 entsprechen und nun nach der Beziehung F(k.t) ~ E3+?/ 
verschwinden. Sich iiber diesen typischen Abklingvorgang der Funktion F(k,t) einer bestimmten 
Wellenzahl Klarheit zu verschaffen, ist fiir die Beurteilung des Turbulenzverhaltens und die 
zeitliche Veranderung der spektralen Verteilung wesentlich. 

Um die Kenntnisse tiber das Einstellen der Verteilung des statistischen Gleichgewichtes zu 
-vertiefen und eine Vorstellung von den hierzu erforderlichen Zeiten zu bekommen, soll die Frage 
untersucht werden, wie sich das Spektrum gegentiber Storungen der Gleichgewichtsverteilung 
verhalt. Wahrend in Abschnitt 3 das Gleichungssystem (3.11), (3.15) und (3.16) beniitzt wurde, 
um bei vorgegebenem d/f'(k,t)/dt das Spektrum zu ermitteln, kann auch umgekehrt damit unter- 


sucht werden, wie sich ein im Zeitpunkt t, vorhandenes Spektrum zeitlich andert. Setzt man 
(3.16) in (3.15) und diesen Ausdruck in (3.11) ein und differentiiert nach k, so erhalt man 


Ss ae BED Be ee ak | 2F (ht) k®— 27% 2) fon (k’,t)k2dk’. (5.8) 


Mit Hilfe dieser Beziehung kann man bei gegebenem Spektrum aie zeitliche Anderung von 
F(k,t) verfolgen. 

- Wir nehmen nun eine iiberall stetige Funktion F'(k,t) an, die beispielsweise der Gleichgewichts- 
bedingung (3.12) geniigen mége, jedoch sei in dem Bereich zwischen k und k-+dk eine unstetige 
Storung AF (k,t) vorhanden. Die Integrale in (5.8) andern ihre Betrage infolge der Storung 
AF(k,t) nicht, da die Breite der Stérung infinitesimal angenommen wurde. Der Wert OF (k,t)/dt 


bleibt im ganzen Wellenbereich unbeeinfluBt von der Stérung bei k, mit Ausnahme an der Stelle k 
‘selbst. Hier wird 


4 [F kt) + Arey] =| yn jyre 


a 2 Lei arto 


re 1) + AF, 52) feran (kyt)h2 dh. 
k3 


Nun sei ferner A F(k,t) <F(k,t). Subtrahiert man von (5.9) he Wert —0F(k,t)d/t nach (5.8) 
‘und eliminiert au®erdem den Ausdruck fiir die Zahigkeit (» $f . dk) mit Hilfe von (5.8), 


so erhalt man folgende Differentialgleichung fir die zeitliche Anderung der Stérung AF (k,t): 


dAF (kt AF (i) aF (kt) , es F (kt) - 
: ) = lle [2Paoreay (5.10) 


Ot F(k,t) ot 
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Aus dieser Beziehung erkennt man unmittelbar, da 
1 aA F (Kk, t) —1. dF¢,t) 
A F (k,t) dt F(k, 1) dt 


ist. Damit ist bereits der Nachweis erbracht, daB eine Stérung AF (kt) schneller abklingt al 
die Funktion F(k,t) im gleichen Wellenbereich. 
Die Gleichung (5.10) 1aBt sich leicht integrieren: 


ie 
t f 2F(k,t\)kedk 
Zid te a 


AF (ht) ~ F(t) exp] 2 rene 


(5. 11) 


Eine kleine Stérung AF an der Stelle k klingt zeitlich wie die ungestérte Funktion F an dieser: 


Stelle, multipliziert mit einer Exponentialfunktion, ab, wobei die Exponentialfunktion natiirlich)| 


das Entscheidende ist. Nehmen wir z. B. eine Stérung bei einem nicht zu kleinen k, bei welchem) 
das Gesetz F(k,t) ~ k-5/> gilt, und nicht zu kleiner i) Meas an, so wird dort mit (3.12)) 


Jytx {yr dk | f2rioueae 3? 


ae t) 


[2Fd.gieak sce. 


Daraus folgt mit (5.7) | 


ke 
Z Ace EMD age 3 E(t) 12/3 1 
: VFeéne Mera Lae ase }, 


ty 


| 
| 
j 
| 


wenn mit ty) der Zeitpunkt bei Beginn der Stérung gekennzeichnet wird, und % =k L(t) ist. Also | 
ergibt sich fiir eine kleine Stérung im Bereich des k—*/*-Gesetzes das Abklinggesetz 


Fash ss E(t) J1+2/36 E ([ EQ) 2/3 
AF (k,t) = AF (k, to) aa exp F@i i| E@ | || 4 (5. 12) | 


Man kann leicht abschatzen, daB z. B. innerhalb eines Zeitraumes, in welchem die kinetische 


Energie auf etwa die Halfte abklingt, eine Stérung A F(k,t) bis auf einen ganz unbedeutenden 
Rest verschwindet. 


Fur sehr kleine k-Werte, bei denen F(k) ~k§—! zeitlich konstant bleibt, ergibt sich fur pide 
kleine St6rung das Abklinggesetz 


AF (k,t) = AF (kt) exp |- sey V FH) 


k-Werte nur sehr langsam abklingt. 


Was hier fir eine kleine értliche Storung gezeigt wurde, gilt auch im wesentlichen fur solche 


Stérungen, die sich tber einen endlichen k-Bereich erstrecken. Die Betrage der Integrale in 
(5.8) werden weniger stark durch Abweichungen von der Gleichgewichtsfigur geandert als die 


78 (¢—t) | (5.13) | 


Aus einer zahlenmaBbigen Abschatzung entnimmt man, da eine Stérung im Gebiete kleiner | 


Betrage von F'(k,t) im Bereich der Stérungen. Die GréBe dF (k,t)/dt wird also hauptsachlich | 


nur in dem k-Bereich durch die Stérungen beeinfluBt, an dem die Stérungen selbst auftreten. 


Man kann aus diesen Ergebnissen folgern, daB sich in einem breiten k-Bereich ziemlich rasch | 
ein Gleichgewichtszustand fiir den Energieaustausch einstellt. Nur im Gebiet kleiner k-Werte | 


behalt ale Spektrum langere Zeit die Form bei, die es im Augenblick seiner Entstehung hat. 
Hochstens die grébsten UnregelmaBigkeiten gleichen sich hier aus. 


Es sollen nun noch einige Beziehungen der hier eingefihrten GréBen zu den GréRen an- | 


gegeben werden, die in dem in England verfaB8ten Arbeiten benutzt werden. Die als GréBe der 
kleinsten Turbulenzelemente anzusehende Lange A betragt 


10 L? 


Te Sines (5.14) 
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In Abb. 4 ist fir das Verhaltnis L/A ein OrdinatenmaBstab eingezeichnet. Die durch Integration 
on (ce) 
der Korrelationsfunktion f(r) erhaltene, als ,,integral scale’ bezeichnete Linge / f(r)dr unter- 


r 

_ scheidet sich von der hier eingefiihrten mittleren Wellenlinge L nach (2.7) rein einen kon- 

_stanten Zahlenfaktor S 

Ay 

4 [fears 22a. (5.15) 

i q 0 

Die Funktion f(r) beschreibt dabei die Korrelation zwischen zwei parallelen Geschwindigkeits- 

_komponenten an zwei Punkten, die auf einer parallel zu diesen Geschwindigkeitskomponenten 

_verlaufenden Linie im Abstand r voneinander entfernt liegen. — 

6. Beriicksichtigung der Loitsianskyschen Invarianten. Uber die GréSe des auftretenden 
Parameters € wurden aus den bisherigen Uberlegungen keine theoretischen Aufschliisse erhalten, 

_ auBer daB €>1 sein mu, um brauchbare Lésungen zu erhalten. Nun wurde eine sehr wichtige 

_ Bedingung, die sogenannte Loitsianskysche Invariante, noch nicht berticksichtigt. An Hand der 

“von Th. v. Kérmén und L. Howarth! aus den Navier-Stokesschen Bewegungsgleichungen her- 

_ geleiteten Beziehung zwischen der Doppel- und Dreifachkorrelationsfunktiou laBt sich namlich 
zeigen, daB der Ausdruck? 


ia AN 


A= = [ firyrtar (6.1) 


_wihrend des Abklingens konstant bleiben mu. Dabei ist f(r) die bereits bei (5.15) naher de- 
- finierte Korrelationsfunktion. Setzt man fiir die Korrelationsfunktion die von W. Heisenberg 
_angegebene Beziehung ein, so nimmt A folgende Form an: 


co oo 2 
A=f f FQj Aa bahar. (6.2) 
3 (pit) 


k3 


Der Betrag dieses Doppelintegrals hangt vom Verlauf von F'(k,t) in der Nahe k=0 ab. G. K. 
Batchelor hat nachgewiesen, daB das Spektrum in der Nahe k=0 die Form F(k) ~ k* haben muB. 
Es zeigt sich also, daB A nur fiir §=5 einen endlichen Wert hat. Fiir €>5 wird A=0. Dies hatte 
nach (6.1) zur Voraussetzung, dal} f(r) in einem gewissen r-Bereich negativ wird, was mit vor- 
liegenden Versuchsergebnissen im Widerspruch steht. Fir <5 ergibt sich andererseits A= © ; 
dies erscheint aber physikalisch ebenfalls unwahrscheinlich. Fur voll ausgebildetes statistisches 
Gleichgewicht im gesamten Wellenbereich 0<k<o, das als asymptotischer Zustand fiir eine 

_hbeliebig entstandene Turbulenz gelten soll, scheint also nur der Parameter €=5 in Betracht 
zu kommen, wofir noch weitere Uberlegungen sprechen. Damit liegen die Eigenschaften der im 
statistischen Gleichgewicht befindlichen Turbulenz fest. Fiir groBe Reynoldszahlen ergibt sich 
das Dissipationsgesetz E ~ t~1/’, welches bereits von A. N. Kolmogoroff? gefunden wurde; fur 
verschwindende Reynoldszahlen gilt E ~ t~*/* in Ubereinstimmung mit theoretischen und ex- 
perimentellen Ergebnissen von G. K. Batchelor und A. A. Townsend*. Der gesamte Ablauf des 
Abklingens laBt sich fir diesen ausgebildeten Gleichgewichtszustand nach den angegebenen 
Beziehungen rechnerisch genau verfolgen. 


Bei Versuchen im Windkanal scheint die Einstellung des statistischen Gleichgewichtes fir 
den gesamten Wellenbereich 0<k<o ziemlich langsam vor sich zu gehen, so dafS’ man nach 
vorliegenden Erfahrungen erst in Abstanden von einigen hundert Maschenweiten vom Turbulenz- 
gitter einen ausgebildeten Gleichgewichtszustand vorfindet. Dabei ist es aber theoretisch nicht 
begriindet, daB sich dieser Zustand erst bei kleinen Reynoldszahlen einstellt, wie es bei den 
Versuchen von G. K. Batchelor und A. A. Townsend der Fall ist. Man sollte vielmehr erwarten, 

da bei Turbulenz, die mit geniigend hohen Reynoldszahlen beginnt, das voll ausgebildete Gleich- 
gewicht schon bei entsprechend hohen Re-Zahlen praktisch erreicht wird. 


Die experimentellen Untersuchungen lassen nun weiter darauf schlieBen, daB sich bei Tur- 
bulenz, die im Luftstrom von einem Gitter erzeugt wird, schon nach verhaltnismaBig kurzer 


1 Th. v. Karman und L. Howarth, Proc. Roy. Soc. A 164 (1938) 5. 192. 

2 G. K. Batchelor, Quarterly Applied Mathematics IV (1948) S. 97. 

3 Vel, die zuletzt zitierte Arbeit von Batchelor. 

4G, K, Batchelor und A, A. Townsend, Proc. Roy. Soc. A 194 (1948) S. 527, 
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Einstellzeit ein teilweise ausgebildeter Gleichgewichtszustand einstellt, der sich mitunter eine 

Weile fast unverandert halt. In diesem Zustand laBt sich das Spektrum fiir den Bereich k>? 

mit guter Annaherung ebenfalls durch eine Gleichgewichtsverteilung nach (3.12) mit eimem- 

passend gewahlten &(& +5) rechnerisch beschreiben. Dabei ist @ jedoch so klein (8 <1/L), dab 
eo 


die im Bereich 0<k<@ enthaltene kinetische Energie / F(k) dk klein gegen die Gesamtenergie 
0 


nach (2.1) ist und deshalb keinen Einflu® auf das Gesamtverhalten der Turbulenz hat. Fir k<@ | 
gibt die Rechnung das wirkliche Spektrum nicht richtig wieder. Die allgemeinen Eigenschaften — 


der Turbulenz sind in diesem Sta- 
dium noch wesentlich von den Be- 
dingungen abhangig, unter denen 
sie erzeugt wurde, also z. B. von 
der Form des Gitters, was durch 


— zu beriicksichtigen ist. Beim Ver- 
gleich mit Messungen kann man & 


rechnerische Abklinggesetz mit dem 


stimmt. Die Messungen von G. K. 
Batchelor und A. A. Townsend! er- 
gaben mit ziemlicher Genauigkeit in 
diesem Stadium ein Abklingen, das 
dem Parameter =2(E~ 1/t) ent- 
spricht. Da jedoch diese Messungen 
alle mit geometrisch vollkommen 
ahnlichen Gittern vorgenommen 
wurden, besteht kein AnlaB, diese 
Ergebnisse zu verallgemeinern. Den 
von Th. v. Karman? veréffentlichten 
Messungen ist zu entnehmen, dah 
die Form des Gitters einen EinfluB 
auf das Abklinggesetz im Anfangs- 
stadium ausiibt. Im _ wbrigen ist 
nicht sicher, ob sich in allen Fallen 
ein teilweise ausgebildeter Gleich- 
gewichtszustand solcher Art bildet, 
daB er tiber einen langeren Zeit- 
raum durch ein konstantes & be- 
schrieben werden konnte. Es _ ist 
denkbar, daB weitere experimentelle 
Abb. 5. Gleichgewichtsspektren yon Windkanalturbulenz. Untersuchungen noch manches 

g=2 (E~r). interessante Ergebnis wtber diese 
Frage liefern werden. 


7. Vergleich mit Messungen im Luftstrom hinter einem Gitter. Durch die Funktion F (k, t) 
sind die statistischen Raummittelwerte der augenblicklichen Geschwindigkeitsschwankungen der 
im Mittel als ruhend zu denkenden Flissigkeit beschrieben. Bei den Messungen im Luftstrom 
hinter einem Gitter wird die Strémung mit ziemlich groBer Geschwindigkeit U am MeBort vor- 
beigefiihrt und dabei der Zeitmittelwert der parallel zu U wirkenden Geschwindigkeitskomponen- 
ten u, gemessen. Oder, um beim Bilde der im Mittel ruhenden Strémung zu bleiben, die MeB- 
apparatur wird mit der Geschwindigkeit U geradlinig durch das Strémungsgebiet gefiihrt, und 
man miBt dabei die Mittelwerte der Geschwindigkeitskomponenten parallel zur Bewegungs- 
richtung der Apparatur. Die spektrale Verteilung der Komponenten uj tber die Frequenz n 


1G. K. Batchelor und A. A. Townsend, Proc. Roy. Soc. A 193 (1948) S, 539. 


2 Th, v. Karmén, Some Remarks on the statistical Theory of Turbulence. V. Int ti 1 Cong 
Applied Mechanics (1938) S. 437. i 2 a eed 


entsprechende Wahl des Parameters 


beispielsweise so festlegen, da das : 


experimentell beobachteten tiberein- | 


sr nd. 


= 
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: wird durch die Funktion F(n) beschrieben, die wie folgt definiert ist: 


© 

Ss E nH 

Wy == =f Pn) dn, (7.1) 
0 

3s kU ; : : . 
dabei ist n = oa: Die Funktion F,(n) errechnet sich nach Heisenberg aus F(k,t) zu 
30 f Bk? 
UF (n) === a F(k,t)dh. (7.2) 


: ky 
Abb. 5 zeigt einige Kurven UF,(n)/227L ttber y=22Ln/U im logarithmischen Mafstab fiir 
&=—2 (E~1/t). 
Beim Vergleich mit Messungen hat man 
. auBer € noch zwei Parameter, tiber die man 
frei verfiigen kann, naimlich die GréBe L 
der Turbulenz und die Konstante x. Mit 
anderen Worten, man kann die Re-Zahl will- 
kirlich wahlen und auBerdem den Nullpunkt : 
auf einer unter 45° geneigten Geraden ver- 
schieben. Die GréBe L hat man so festzu- 
legen, dafB die gerechnete Kurve im Bereich 
kleinerer k-Werte mit den Messungen tber- 
einstimmt. Aus der Anpassung der Rechnung 
an die Versuchsergebnisse beim Ubergang in 
das Gesetz F'(k,t)~k~’ ergibt sich der Wert x. 


Aus Abb. 6 ersieht man, da die Rechnung 
in sehr gute Ubereinstimmung mit den Mes- 
sungen nach L. F.G. Simmons und C. Salter! 
zu bringen ist. Nach vorliegenden Versuchs- 

erfahrungen? liegt der MeBort (Abstand vom + U= 15 ft/sec 
Gitter =27 Maschenweiten) in einem Bereich, Mey acon 
in welchem noch kein voll ausgebildeter 
Gleichgewichtszustand besteht, das Spek- 
trum also noch in gewisser Weise von 
den Anfangsbedingungen abhangt. Da das 
Dissipationsgesetz ftir die Messungen von 
Simmons und Salter nicht genau bekannt 
ist, wird mit Ricksicht auf die Ergebnisse 
von G. K. Batchelor und A. A. Townsend, 
bei denen fiir ein Gitter mit dem Verhalt- 
nis Stabdurchmesser/Maschenweite = 3/16 das 
Gesetz E~t-1 beobachtet wurde, fiir das 
hier vorliegende Verhaltnis Stabdurchmesser/ Abb. 6. Spektrum von Windkanalturbulenz 
Maschenweite=0,29 das gleiche Abkling- pack) Recbauny wad Messong. 
gesetz, also €=2 der Vergleichsrechnung 
zugrunde gelegt. Weiter gehéren zur Rechnung die Werte x=0,28 und® 27L=0,25 ft. Die 
mittlere Wellenlange 22%L der Turbulenz hat rechnerisch die Grébe der Maschenweite und 
ist also klein im Vergleich zu den Abmessungen des Strahles (Durchmesser 4 ft), so das man 
wohl annehmen darf, daB von der Berandung des Stromungsgebietes kaum merkliche Einfliisse 
auf die Turbulenz in Kanalmitte ausgeiibt werden. Angesichts dieser Tatsache und da ferner 
die Turbulenz tatsachlich isotrop ist, wie aus den gemessenen Korrelationsfunktionen zu folgern 
ist, scheinen die Voraussetzungen der Theorie einigermaen erfiillt, zu sei. Da eine befriedigende 
Ubereinstimmung zwischen Rechnung und Messung zu erzielen ist, besteht somit einige Wahr- 


+ U=30 ” 


1 LL, F. G. Simmons und C. Salter, Proc. Roy. Soc. A 165 (1938) S. 73. 

2 Vel. FuBnote 1 von S. 74. 

3 W. Heisenberg hat auf Grund der gleichen Messungen die GréBe x auf etwa 0,5 geschatzt. Der Unterschied 
diirfte im wesentlichen auf die Abweichung zuriickzufiihren sein, die zwischen dem von Heisenberg numerisch 
berechneten Spektrum und der mathematisch exakten Form (4.5) besteht. 


76 Rotta: Das Spektrum isotroper Turbulenz im statistischen Gleichgewicht. Ingenieur-Archiy 
fe as am nrnrrinrr yea EINSTEINS 


scheinlichkeit, daB die Annahmen der Rechnung den physikalischen Gegebenheiten gerecht | 
werden und am Ort der Messung tatsachlich weitgehend ein statistischer Gleichgewichtszustand | 
in dem definierten Sinne besteht. Zur Bestatigung waren allerdings Vergleiche mit weiteren 
Messungen nitzlich, insbesondere um Unterlagen fiir die Konstante x zu bekommen. 


8. Zusammenfassung. Der von W. Heisenberg angegebene Naherungsansatz fir die Turbu- 
lenzzahigkeit erméglicht eine neue Behandlung des Turbulenzproblems. Aus der Tatsache, da _ 
das Turbulenzspektrum, d. i. die Verteilung der kinetischen Energie der turbulenten Bewegung — 
auf die Wellenzahlen (reziproke Wellenlange), im Laufe der Zeit einem statistischen Gleichgewicht _ 
zustrebt, in welchem seine Form von den Anfangsbedingungen unabhangig ist, wird eine Be- — 
ziehung zur Ermittlung des Spektrums hergeleitet. Die geringe Anderung des Spektrums mit 
der Reynoldszahl der Turbulenz erlaubt eine mathematische Vereinfachung. Die Loitsianansky- 
sche Invariante liefert die noch fehlende Aussage, um von der Vielzahl der erhaltenen Lésungen — 
eine ausgezeichnete festzulegen, die die spektrale Energieverteilung sowie das Dissipationsgesetz 
fiir isotrope Turbulenz im. ausgebildeten statistischen Gleichgewicht angibt. 

Bei geringen Abstinden vom Turbulenzgitter (20 bis 100 Maschenweiten), bei denen die meisten 
bekannt gewordenen Messungen von Spektren und Korrelationsfunktionen vorgenommen wur- 
den, besteht im allgemeinen noch kein voll ausgebildeter Gleichgewichtszustand. Spektrum und 
Dissipationsgesetz sind in diesem Stadium noch von den Bedingungen, unter denen die Turbu- 
lenz erzeugt wurde (Geometrie des Turbulenzgitters), abhangig. Doch hat sich in der Regel schon 
so weit ein Gleichgewichtszustand ausgebildet, daB der Vergleich solcher Messungen mit be- 
rechneten Spektren bei Auswahl einer passenden Lésung (passendes &) méglich ist. Ein der- 
artiger Vergleich mit von L. F.G. Simntons und C. Salter gemessenen Spektren zeigt befriedigende 
Ubereinstimmung. 


Max-Planck- Institut fiir Strémungsforschung. 


(Eingegangen am 25. Februar 1949.) 
Anschrift des Verfassers: Julius Rotta, (20b) Géttingen (Hann.), BottingerstraBe 6/8. 
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handlung des Anzahlbegriffs. Darstellung der Grundbegriffe der Mengenlehre im erweiferten 
Kalkiil. Die logischen Paradoxien. Der Stufenkalkil. Anwendungen des Stufenkalkiils. Ab- 
schlieBende Bemerkungen zum Stufenkalkiil. 
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